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Introduction 



One of the main open question in physics is the understanding of the hadrons as constituted of 
quarks and gluons. In effect, the idea that quarks are the building blocks of hadrons goes back to 
1964, when Gell-Mann and Zweig published their famous papers |96 tl207j . In 1973, the formulation 
of Quantum ChromoDynamics (QCD), as the theory of the strong interaction between quarks, was 
proposed by Fritzsch, Gell-Mann and Leutwyler [95]. But despite the fact that QCD has been 
adopted 35 years ago and even though this theory has been verified at high energy, we have not 
been able to reveal the detailed structure of hadrons. One might wonder why we have not achieved 
this aim yet. The answer is that the peculiar features of the strong interactions renders such a task 
very intricate. 

At high energy, the perturbative treatment of QCD, which is justified by asymptotic free- 
dom [1061 1169j , allows for the explanation of the hadronic phenomena, whose basic ingredients are 
the Parton Distributions. These distributions can be experimentally determined for a relatively 
high value of the momentum transfer Q^. Therefore, we are able to evolve those quantities and 
give predictions for different values of Q^. 

At low energy (typically below ~ 1 GeV), there is no justification for such a perturbative 
treatment of QCD and, therefore, we do not longer have a description of the relevant low energy 
observables from QCD. In this regime, hadronic models and effective theories come into play. In 
these schemes, the parameters are fixed by phenomenology. 

There exists no clear and definite frontier between these two energy regimes. Neither do exist 
definite connections between the relevant degrees of freedom: even though they are the elementary 
degrees of freedom of QCD, free quarks and gluons have never been observed. This fact has re- 
vealed that quarks and gluons were confined in hadrons. On the other hand, another important 
feature of the strong interactions is related to the chiral symmetry. In the limit of massless quarks, 
chiral symmetry, which is a property of the QCD lagrangian, is realized in the Goldstone mode, 
i.e. is spontaneously broken. The fundamental role played by this symmetry and its realization 
has been confirmed by the small mass of the pion as well as properties of the strong processes at 
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intermediate energies like the Partial Conservation of the Axial Current or the Soft Pion Theorems. 
Therefore, the structure of a hadron arises from an interplay between these two important properties 
of QCD, namely, confinement and the realization of the chiral symmetry. The energy regime from 
which we can gather information on the hadron structure resides in the frontier between the two 
regimes described above, where both confinement and chiral symmetry are supposed to be into play. 

A way of connecting the perturbative and non-perturbative worlds is the study of Parton Dis- 
tributions using models. The scheme that has been proposed runs as follows: we build models, 
whose characteristics somehow mimic the QCD's properties that are relevant at the energy range 
we are considering (~ 1 GeV); we then evaluate the Parton Distributions in these models; finally, 
we evolve these distributions to the scale of the experiment in order to compare with the data. 

The reactions allowing to access the different types of Parton Distributions have been receiving 
great attention from the hadronic physics community. These reactions are such that they enable us 
to look with a good resolution inside the hadron and allow us to resolve the very short distances, 
i.e. small configurations of quarks and gluons. Since at short distances the interactions between 
quarks and gluons become weak, this part of the process is described through perturbative QCD. 
A resolution of such short distances is obtained with the help of non-strongly interacting probes. 
Such a probe, typically a photon, is provided by hard reactions, like Deep Inelastic Scattering, 
Semi-Inclusive Deep Inelastic Scattering, Deeply Virtual Compton Scattering, .... In that scheme, 
the Parton Distributions refiect how the target reacts to the probe, or how the quarks and gluons 
are distributed inside the target. The insight into the structure of hadrons is reached at that stage: 
the large virtuality of the photon, Q^, involved in such processes allows for the factorization of 
the hard (perturbative) and soft (non-perturbative) contributions in their amplitudes. It is this 
non-perturbative part of the process that will interest us all along this thesis. 

The first steps towards the understanding of the structure of hadrons were taken by studying 
the total cross sections of fully inclusive processes or longitudinal asymmetries, which receive simple 
interpretations in parton models |113j . 

Prom a theoretical point of view, Parton Distributions are related to diagonal matrix elements of 
a bilocal operator on the light-cone from the initial hadron state to the same final hadron state, i.e. 
the same particle with the same momentum. They have simple probabilistic interpretations. For 
collinear quarks, the Parton Distributions q{x) receive the name of number densities because they 
reflect the probability density of finding a quark with a fraction x of the longitudinal momentum 
of the parent hadron, regardless of its spin orientation. Also, the definition of other appealing 
quantities arise, such as the helicity distribution (i.e. longitudinal polarization), as well as the 
more mysterious transversity distribution. Now, the distribution of transverse polarization is not 
observable in fully inclusive processes, thus, it requires the analyses from other hard processes, e.g. 
semi-inclusive processes. 

A first proliferation of distributions was the consequence of considering non-collinear quarks. It 
basically means that the intrinsic motion of quark could hardly be understood through the three 
previous densities. Nowadays, Transverse Momentum, kx, Dependent parton distributions are be- 
ing thoroughly examined in the context of Spin Physics. Since they embody an additional internal 
degree of freedom, they can be involved in non-trivial correlations with spin that might lead to a 
deeper understanding of the spin structure of hadrons. 
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In these days, the experimental facihties allow to access not only inclusive processes like Deep 
Inelastic Scattering but also exclusive deep processes like the Deeply Virtual Compton Scattering 
(DVCS) 7*p — )• 7p or Hard Exclusive Meson (M) Production (HEMP) ■j*p — )■ Mp. The first project 
of a facility devoted to the study of exclusive processes in the scaling regime was proposed as the 
ELFE project [181 1129| . which unfortunately did not succeed. In the last few years, experimental 
data related to the structure of hadrons, like asymmetries and DVCS or HEMP cross sections, have 
been obtained from different collaborations. We mention the HERMES, e.g. [4j, as well as the HI, 
e.g. [Hj, and the ZEUS, e.g. [55j, collaborations at DESY; the CLAS collaboration, e.g. |191j . at 
the Hall B of JLab. 

The theoretical description of these processes requires the generalization of the Parton Distribu- 
tions to Generalized Parton Distributions \in\ 11481 1173j . These distributions are matrix elements 
of non-diagonal distributions, namely, those which have a final state with different momentum 
than the initial one (for reviews, see, e.g., Refs. [30l [36l EH IIOIL I120j ). The Generalized Parton 
Distributions describe non-forward matrix elements of bifocal operators on the light-cone. They 
measure the response of the internal structure of the hadrons to the probes in the above-mentioned 
processes. These new distributions have been shown to be interesting theoretical tools for the study 
of hadrons; they connect, through sum rules, to the hadronic Form Factors. Also their forward 
limit connects to the usual Parton Distributions. The properties of such matrix elements that ensue 
from their analysis from first principles are presented in Chapter [T} Those properties are used as 
constraints on the modeling procedure that has been described above for the Parton Distributions. 

Already the prediction of some phenomenological models have been verified to be in agreement 
with the data mentioned above. This was taken as an encouraging sign that the framework of 
exclusive deep processes was adapted to our purposes. However, a detailed mapping as well as the 
extraction of GPDs from the experiments will require more data in order to constrain the models. 
In a near future, we expect new data on DVCS from the CLAS and HERMES collaborations (see 
references of proposals in [108J). The improvement up to 12 GeV of the CEBAF accelerator at JLab 
foreseen for 2012 [54j is already involving much of the hadronic physics community in proposals for 
the access of the involved distributions as a way to explore the structure of the hadrons. So is the 
COMPASS experiment at CERN. 

Also, collisions of a real photon and a highly virtual photon can be an useful tool for studying 
fundamental aspects of QCD. Inside this class of processes, the exclusive meson pair production in 
7*7 scattering has been analyzed in Ref. [165 1 introducing of a new kind of distribution amplitudes, 
called Transition Distribution Amplitude (TDA). They represent a generalization of parton distri- 
butions to the case where the initial and final states correspond to different particles. For obvious 
practical reasons, the first transitions that have been studied are the (mesonic) pion-photon TDA, 
alternatively governing processes like vr+vr^ — t- 7*7 or 7*7r^ — t- ^tt^ in the kinematical regime where 
the virtual photon is highly virtual but with small momentum transfer. Nevertheless, it has also 
been proposed in Ref. |165j baryonic TDAs that have been analyzed for the ir ^ N transition in 
Ref. il35j. 

This thesis aims to contribute to the understanding of the hadron structure by gathering the 
present knowledge on the pion-photon Transition Distribution Amplitudes. 
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For that purpose, we present the results for the pion-photon TDAs in a field theoretic scheme 
treating the pion as a bound state in a fully covariant manner using the Bethe-Salpeter equation, 
with the pion structure described by the Nambu - Jona Lasinio (NJL) model [671 EHl EH] . Actually, 
the NJL model is the most realistic model for the pion, based on a local quantum field theory 
built with quarks. It respects the realizations of chiral symmetry and gives a good description of 
the low energy physics of the pion |125| . The NJL model is a non-renormalizable field theory and 
therefore a cut-off procedure has to be defined. The Pauli-Villars regularization procedure has been 
chosen because of the symmetries it respects. Also, the NJL model together with its regularization 
procedure, is regarded as an effective theory of QCD. 

In terms of the physical processes, the NJL model is used to describe the soft (non perturbative) 
part of the deep processes, while for the hard part conventional perturbative QCD must be used. 

For the sake of clarity, we will develop this topic step by step, following the theoretical improve- 
ments described above. Namely, in order to illustrate the framework of Parton Distributions as 
well as the formalism developed to calculate such quantities, the detailed calculation of the pion 
Distribution Amplitudes in the NJL model is presented in Chapter [2j In the same Chapter, the 
link between pion Distribution Amplitudes and Parton Distributions are displayed through the 
Soft Pion Theorem. We also show the consistency of the model we use by emphasizing that, once 
QCD evolution is taken into account, a good agreement is reached between the calculated Parton 
Distributions and the experimental one. 

According to the same scheme we extend our considerations to the pion Generalized Parton 
Distributions in Chapter [3j The results for the calculation of the latter in the formalism of Chap- 
ter [2] are shown |153j . One could wonder whether the results, even if constrained by the first 
principles properties, are consistent with other formalisms. With that question in mind, we display 
the parameterization of the GPDs through the so-called Double Distributions. The last step of the 
analysis is the application of the QCD evolution equations. The effect of evolution on the pion 
GPD calculated in the NJL model has been studied in [35], whose results we will be discussed in 
detail. 

These previous developments enable a better prehension of the Transition Distribution Am- 
plitudes as non-perturbative objects matching with some well-known constraints. A peculiarity 
of these transitions consists in the connection of their structure with the one given by the pion 
radiative decay tt^ — )• je~^u process |46 1 1147j . which stresses the importance of the role played by 
the Partial Conservation of the Axial Current. A detailed analysis of the matrix element of bilocal 
current defining the pion-photon TDAs is presented in Chapter |4j Their properties are displayed 
and studied. The results of the calculation in the NJL model are subsequently given. 

Different studies of the axial and vector pion-photon TDA have been performed using different 
quark models \i4\ 1671 11271 1199] ■ A qualitative comparison as well as the parameterization through 
Double Distributions are also presented, which allow us to conclude that there is a qualitative 
agreement between all the different approaches. 

The direct application of the results obtained for the TDAs to phenomenology is obviously the 
estimation of the cross section for the processes they govern. Chapter [5] is devoted to the study of 
the cross sections for the exclusive meson production in 7*7 scattering |701I131| . where, in Ref. |70j . 
the evolution of the TDA is taken into account. 
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Up to this point, we have proposed to consider, in this thesis, the structure of the pion. Within 
the study of the hadron structure that we have carried out, we wanted to extend our analyses to 
the spin structure of the proton. This is why we briefly report some recent improvements in Spin 
Physics. From the phenomenology, it is known that Semi-Inclusive DIS off a transversely polarized 
target shows azimuthal asymmetries, e.g. OE]) that could be understood through non-trivial cor- 
relations between intrinsic motion of the quarks and transverse spin. A hint on the spin structure 
of the proton could be given through the understanding of the modulation of the number density 
of unpolarized quarks in a polarized proton. This modulation is due to the correlation between the 
intrinsic transverse momentum of the quarks and the transverse component of the proton spin. The 
results for this so-called Sivers function, calculated in two different models [66^ I72j. are exposed in 
Chapter [6j 

In the final Chapter, we conclude by gathering all the questions and all the ideas that emerged 
in this work, which represent good starting points for future research. 
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Hadron Structure from Distribution Functions 



An intuitive approach of parton distribution can be provided in the context of the parton model. 
Feynman [90] introduced the parton distributions as phenomenological quantities describing the 
properties of the nucleon in high-energy scattering. Actually a nucleon moving at the speed of light 
can be seen as made of noninteracting massless and collinear partons. In such a scheme the parton 
distributions are the density of partons as a function of the fraction of the nucleon longitudinal 
momentum they carry. 

Nevertheless the nucleon does not exactly consits of free partons. Knowledge of the internal 
structure of the nucleon, or the pion, is required in order to properly study the parton distributions 
as a nonperturbative object. 

Prom the phenomenological point of view, the study of Deep Inelastic Scattering processes, 
where individual partons are resolved, revealed itself to be an essential tool in this context. The 
unpolarized parton distributions have been extracted with good precision from various high-energy 
scattering data [98^ I130j . Those distributions have provided important clues about the hadron 
structure but yet important pieces of information are missed out in these quantities. 



Figure 1.1: DIS: Kinematics of lepton-hadron interaction dynamics. It is expressed as the com- 




X 



In order to study the properties of hadrons we will 
examine the properties of the hadronic matrix element 
defining them. The adapted framework for these pur- 
poses is any Deep Inelastic Process, characterized by 
the Bjorken limit. The cross section for those pro- 
cesses |1.1| can be separated into a leptonic l^u and a 
hadronic tensors W^'^. While the leptonic tensor is 
known completely, W^'^ , which describes the internal 
structure of the nucleon, depends on non perturbative 



scattering. 
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Hadron Structure from Distribution Functions 



mutator of two currents defined at 2 different space-time points, 

W^'' = ^ j d''ze'^-\p,S\[J^{z),rm\P,S) , (1.1) 

and can be related to the imaginary part of the forward Virtual Compton Scattering through the 
Optical Theorem. 

1.1 On the Light Cone 

Without loss of generality, any 4- vector could be expressed as, see Appendix. |Aj 

= (a+,a-^,a") , (1.2) 
with a"*" = (a^, a?) and with the components = -^{cl^ i a?). 
We define two vectors whose norm is zero 

= ^(1,0,0,1) =p+(l, 0^,0) , 

^ ^(1,0,0,-1) = ^(0,0^,1) . (1.3) 



Those vectors follow the properties ( |A.9 ) 



The Sudakov decomposition for a vector 

a^^ = apf" + a^^" + (3n^' , 

= {a.n)p^ + a-^^" + {a.p)n^' , 

with a^^' = (0,5-^^,0), can be adopted for the vectors q, z, with q constrained by the DIS kinematics. 
In terms of the Sudakov vectors, the nucleon momentum p and the momentum of the virtual photon 
q can be written as 

P^ = P^ + —71^" , 

q^ —4 {p-q)n^^-xp^ , 
with X = Oi^ /2p ■ q the Bjorken variable. 

One could inquire about the restriction imposed by the kinematical regime of DIS on the domain 



of validity of the expression of the hadronic tensor ( 1.1 ). The argument of the exponential is crucial. 



For a generic 4-vector z^ = azP^ + z'^^ + Pzn^, the product q ■ z reads, in the Bjorken limit, 

q ■ z —4 {p ■ q) az - X /3z , 



When {p ■ q) goes to infinity, the oscillations produced by the exponential term in (1.1) imply that 



Oz < l/(p • <?). Also, causality implies that the commutator [J^(z), J'^(O)] vanishes unless > 0, 



therefore z-^"^ < 2azl3z- These two restrictions together imply that the hadronic tensor (1.1) in the 



1.2 Ordering 
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Bjorken regime is dominated by light-like distances, i.e. ~ 0. 

Since the dominant contribution to DIS comes from the hght-cone, it is natural to consider 
a dynamical formulation in which the light-cone plays a special role, namely: the light-cone dy- 
namical form [331 1126j . The vector decomposition (1.2) actually corresponds to the light-cone 



decomposition, with the light-like vectors (1.3) 



1.2 Ordering 

The domain of validity of the approach has been defined. So has been the treatment of the 
dynamics on the light-cone. We can now start to investigate on the consequences /advantages of 



these definitions by analyzing the operator describing the hadronic tensor, W^^ (1-1). 

Since DIS processes are light-cone dominated, we can apply the Wilson's Operator Product 
Expansion (OPE). A product of local operators A{x) and at short distance, i.e. x — y ~ 0, 

can be expanded in a serie of well-defined local operators Oj(x) with singular c- number coefficients 
Ci[x) for (i = 0,1,2,3,...) [T^ 



Aix)Biy) = J2 C^{x - y)d,{^) . (1.4) 

i=0 



The OPE basically represents a separation between short-distance, i.e. the singular coefficients, 
and large-distance, i.e. the local operators, objects. It is a formal step towards the factorization 



property. For instance, the Operator Product Expansion Eq. (1.4) of simple bilocal currents is, see 
for example |150j . 

J-(^)J(0) = ^C»(z2)^m...^/^.o^^^^^j0) , (1.5) 



with Cn\z'^) singular c-numbers called Wilson coefficients. The terms in the expression (1.5) 



are arranged in the order of decreasing singularity, the leading contribution being the first term, 
providing for an expansion parameter. 

The strength of light-cone singularity can be found by simple dimensional analysis. If we call 
dj the dimension of the currents, d^Q {n) the dimension of the operator ©^^...^^(O) and n the spin, 
the dimension of the singular coefficients is 

C»(z2) ~ (,2)-'^. + ^-t . (1.6) 



The strength of the singularity is then governed by = d^^ (n) — n and we call the twist. The 
smaller the value of the twist, the higher the light-cone singularity: we can use the twist as an 



^Following Dirac [84] there are three independent parameterizations of space and time, namely, three different 
forms of dynamics which are completely equivalent; namely the instant, front and point forms. The most familiar 
form is the instant form, where the quantization is specified on a space-time hypersurface at a given instant time 
x° = 0. In the light-cone approach, called the front form, the theory is quantized at a particular value of the light-cone 
time x'^ . The quantization conditions are specified on a hypersurface that is a tangent plane to the light-cone defined 
at a light-cone time a::^ = 0. 

^Light-cone dominated does not mean short distance, but the light-cone domination can be led to short distances. 
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expansion parameter. In QCD the lowest twist operator has = 2. A more informal definition of 
the twist is often used |114j . The twist is the order in M/Q, M being the target mass, at which 
an operator contributes to DIS processes. In other words we have {M/QY"~'^. 

An alternative way of isolating the leading-twist contri- 



bution of the hadronic tensor W'^'^ (1-1) is obtained by 



decomposing the hard and soft parts like in Ref. |114j . 
The vector current associated with a fermion field is 



: iIj{z)j^iIj{z) 



with : O : the normal ordering. We will need the identity 

[ViiV'i,'(A2V'2] = V'llV'i, V'2}V'2 - i'2{tp2,'4'i}ipi, as weU as 
the following one, for massless fields. 




Mz),i^{o)} 



1 

2^ 



e{zo)S{z' 



The commutator of two currents appearing in (1.1) can 
be expressed in terms of bifocal operators 



Figure 1.2: Handbag diagram. Most of the 
Feynman diagrams in this thesis have been 
drawn using JaxoDraw [33| . 



[j''iz),rm 



1 



1 



{d,e{zo)Siz')) [-tr(7V7V) 7a V'(O) - V^(0) 7. V'(^)) 

-ie'^P"'' {ij{z) 7„ 75 V(0) + HO) 7„ 75 ^(^)) 



(1.7) 



The hadronic matrix element, decomposed in the following way, 

2 



1 
Air 



d^ze'''-'dp{5iz^)eiz°)) 
tr (7'^7V7") {P, Smz) 7« ^'(O)!^, S) + . . . 



(1.8) 



with the ellipses representing the other terms present in (1.7), gives rise to the handbag diagram 



depicted on Fig. 1.2 The hard part comes from the quark propagator between and z, whereas 
the soft part is parameterized through a bilocal matrix element represented by the lower blob. The 
latter object is related to the Parton Distributions. 



The expression (1.8) has been shown without writing explicitly the color indices. However, the 



latter are implicit and, in order to preserve gauge invariance, one has to take into account the 
propagation of the quarks in the gluon background by linking the two quark fields through the 
Wilson link, i.e. 

i^iz) rv^(o) ij{z) V fexp i [ de A^io] rv^(o) . 



Since the (5-function in (1.8) selects the light-cone, we can expand the gauge link and keep only the 



leading-twist contribution. In the usual parton distribution, the leading-twist contribution involves 
the plus-component of the gauge field. An appropriate choice of the gauge reduces the Wilson 
link to unity: the light-cone gauge A'^ = is chosen that allows the interpretation of the parton 



1.3 Deeply Virtual Compton Scattering 
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distributions as probability densities. 



We assume that the handbag diagram contribution to the hadronic tensor is dominated by small 
values of the transverse momentum of the quark. We therefore consider collinear quarks. Hence, 
to leading- twist, the (unpolarized) Parton Distributions read 

^>-p+-(p|V;(0,z-,0x)7+V'(0)|p> = ^qix) . (1.9) 



2tt 



1.3 Deeply Virtual Compton Scattering 



7*('/) 




Figure 1.3: Virtual Compton Scattering 



The Forward Virtual Compton Scattering (FVCS) is re- 
lated to the hadronic part of DIS through the optical the- 
orem. What generally refers to as Virtual Compton Scat- 
tering (VCS) is any process where 2 photons are involved 
and where at least one of them is virtual. Virtual Compton 
Scattering in the Bjorken region was proposed by Ji [117J 
and Radyushkin |173j as a tool to extract new structure 
functions of the nucleon. Since in this limit the momentum 
transfer is large, we will refer to VCS in the Bjorken re- 
gion as Deeply- Virtual Compton Scattering (DVCS). The 
DVCS provides for a new ground to explore the quark and 
the gluon structure of the nucleon. 



A new kind of distributions, namely non-diagonal distributions, called off-forward, skewed or 
generalized parton distributions (GPDs), were introduced by Ji [117] and Radyushkin |173j from 
the DVCS amplitude 

r^^(P + P',g,P-P') = i j (fze-^''-'{P'\Tf{z)f{<d)\P) , (1.10) 

whose sum rule are related to the Form Factors of the target. 



Following the example of the analysis of DIS in the previous Section, we would like to determine 
the dominant contributions to DVCS. 



1.3.1 Definition of the Generalized Parton Distributions 

The kinematics are similar to the DIS case. Nevertheless the off- forward nature of the new dis- 
tributions requires an additional variable for the description of the initial and final parton states, 
called the skewness variable, ^, which reflects the longitudinal momentum asymmetry. Hence the 
Generalized Parton Distributions are functions of two more variables in comparison to the usual 
Parton Distribution Functions, namely they depend on x, in this context the variable x is just the 
Fourier variable of the light-cone distance, ^ and the four-momentum transfer (squared) t. 
In this work, we choose to follow Ji's notations instead of Radyushkin'^ Hence we introduce the 

^ The light-cone gauge is not always the most adapted choice. For instance, leading-twist Final State Interactions, 
i.e. the exchange of a gluon between the active quark and one of the "spectator", after scattering by the virtual 
photon would rather be studied in a non-singular gauge [55]. This comment will make sense in Chapter [sj 

''For a comparison of the variables see jl03) . 
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average 4- momentum = {P + P')^^ /2 which has to be coUinear to in the z direction. The 
decomposition on the hght-cone is |107j 



pp. 



91 
4e 



(1.11) 



The skewness variable is defined as ^ ~ Q'^/{Ap • q) and is related to the Bjorken variable x in the 
Bjorken limit. Here the momentum transfer t = is negative. 



Similarly to the expression (1.8) obtained for the DIS process, factorization for DVCS can be 



studied |119j . Therefore the DVCS amplitude Eq. (1.10) fulfills its role in being expressed as a 
convolution of hard and soft parts |117j 



2 J_i \x-^ + ie x + E, + te 

Hix,^, t)u{P')iu{P) + E{x,(,t)uiP') !^!!5^u(P) 



2M 



1 



1 



2 ^ ' J^i \x-i + ie x + i + ie 

~ A • n 

H{x, e, t)u{P') i ^5U{P) + E{x, e, t)—u{P'hMP) 



(1.12) 



with u{P) the nucleon spinor. 

This decomposition corresponds to the handbag dia- 
gram depicted on Fig. 1.3 The dominant subpro- 



cess Fig. |1.4| consists on a virtual photon with mo- 
mentum qf^ absorbed by a single quark with momen- 
tum k^. Subsequently the quark radiates a real pho- 
ton, whose momentum is now — A^, and falls 
back with a resulting momentum + A'^ in the 
nucleon. The momentum of the recoil nucleon is 




Figure 1.4: Kinematics of the Deeply- 
Virtual Compton Scattering |117) . 



The four distributions H{x,S,,t), E{x,(,,t), H{x,S,,t) and E{x,S,,t) represent the soft part, 
namely the Generalized Parton Distributions. They are defined through the matrix elements of 



1.3 Deeply Virtual Compton Scattering 
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quark operators at a light-like separation. To leading twist, we have 





dz 


II 




1 




'h{ 


2p+ 






dz~ 


\l 




1 




Hi 


2p+ 





z+=z^=0 

' ^ ' 2M 



u{P) 



H{x, t)u{P')^+ron{P) + E{x, t)u{P') -^^P) 



(1.13) 



The two distributions H(x, ^, t) and H(x, ^, t) conserve the nucleon helicity while the two others 
are said to be (nucleon) helicity-flipping. 

1.3.2 Properties 

Let us overview the properties of these new objects that follow from the analysis from first principles. 
Support 

The support in x of the GPDs is [—1,1]. A negative momentum fraction would correspond to an 



antiquark. In Fig. 1.4, the active quark has a (+)-component momentum fraction x + From the 

A+ 



definition of the skewness variable, 



P 



P+ + P'- 



2p^ 



(1.14) 



and from the fact that (+)-momenta of physical states cannot be negative, it follows that the 
physical region for ^ is the interval [—1, 1]. However, from the light-cone decomposition, we find 
that ^ is bounded by 



< C < 



-t 



< 1 



2^M2 - I 

Given the intervals in x and in ^, we can distinguish 3 regions in x, 



(1.15) 



In the region x G 1] both momentum fractions x + ^ and x — are positive. In this region 
the GPD describes the emission and reabsorption of a quark. 

The behavior is similar in the region x £ [—1, . Both momentum fractions x + ^ and x — ^ 
are negative and therefore the GPD describes the emission and reabsorption of an antiquark, 
respectively with momentum fractions ^ — x and — ^ — x. This is shown on the extreme sides 
of Fig. fTsl 



• In the third region, i.e. x £ [— C??]; one has that + ^ is positive but x — ^ negative. This can 
be interpreted as the emission of a quark with momentum fraction x + and the emission of 
an antiquark with momentum fraction ^ — x both emitted from the initial nucleon. 
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Forward Limit 



When P' = P, the skewness variable ^ is equal to zero. As a consequence GPDs, in the regions 
corresponding to emission/reabsorption of a quark or an antiquark, i.e. [—1, and 1], have an 
equivalence with usual PDFs. This is not true for the region of emission of a quark and an antiquark 
since the interval [— ^, ^] reduces to zero in the same limit. This means that, in the forward limit, 
when the initial and final states are equal, i.e. equal helicities and P' = P, Generalized Parton 



Distributions as defined in Eq. (1.13) are expected to reduce to ordinary Parton Distributions 



Eq. (1.9), 



q[x) 
-qi-x) 



(1.16) 



for X > and a; < respectively. 



When considering the hadronic tensor in DIS (1.8) we have omitted the spin dependent distri- 
bution function as well as the gluon distribution functions. However it seems useful to notice that, 
in the forward limit, the second helicity-conserving distribution is related to the spin dependent 
parton distribution: H{x, 0, 0) = Ag(x) for x > 0. 



Symmetries 

Time reversal interchanges P with P' and its invariance ensures the distribution will not change. 
By so doing we change the sign of ^ and find that the GPDs are symmetric in ^ 

H{x,^,t) = H{x,-C,t) . (1.17) 

Quark distributions are neither even nor odd in x. 

Sum Rules 

On the other hand, forming the first moment of the GPDs by integrating over the momentum 
fraction x, we get the sum rules 

j dxH{x,i,t) = Fi (t) j dxE{x,i,t) = F2 (t) 

j dxH{x,i,t)=gA{t) j dxE{x,i,t)=gp{t) , (1.18) 



1.3 Deeply Virtual Compton Scattering 
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where the dependence on ^ drops out. The integration over x of the matrix element Eq. (1.13) 
removes all reference to the particular light-cone direction with respect to which ^ is defined. As a 
consequence the Form Factors are ^-independent. 

The Form Factors Fi(t) and F2{t) are respectively the Dirac and Pauli Form Factors 



(P'|^(0)7^V'(0)|P) 



u{P') Fi(t)7'^ + F2(t) 



2M 



u{P) 



and gAit) and gp{t) the axial and pseudoscalar ones 



{p'\i;{o)ri5m\p) 

defined for each separate flavor. 



uiP') ga{th^i5 + gpit) 



2M 



uiP) 



[1.19) 



(1.20) 



Gluon Distributions and the Ji's Sum Rule 

We can define the gluon GPDs in a similar way as the quark GPDs 





fdz- 


J 


2tt 


1 


W{x 


2p+ 






fdz- 


¥ ) 


2tt 


1 




H3{x 


2p+ 





Jxz P (^p'^Q+ 



z+=z^=0 



H^ix, t)u{P')j+u{P) + E3{x, t)u{P 



2M 



MP) 



e-^ P-(P'|G+(-^)G+(-)|P) 



Ha{x, t)u{P'h+j5u{P) + E3{x, mP 



2M 



u{P) 



(1.21) 



An interesting sum rule concerns the quark and gluon spin contribution to the nucleon spin 
as it has been one of the main motivation for GPDs. The second moment of both //(x, ^, t) and 
E{x,^,t), the helicity conserving distributions as well as the first moments of the same gluonic 
distributions are directly related to the total angular momentum contribution to the nucleon spin 

J [m] 



1 

2 



As a consequence, the total angular momentum carried by the quarks or the gluons inside a proton 
could be known if one takes the particular value of the desired Form Factors at the forward limit 
t = 0. Since the spin part of the angular momentum is found through the helicity distributions 
A(7(x), the Ji's sum rule enables to access the orbital angular momentum of the quarks and gluons. 
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1.3.3 Polynomiality 

By writing down the most general expression in terms of relevant vectors, the first moments are 
immediately related to their Form Factors. The higher-moments of the GPDs are accessed through 
the tower of twist-2 operators, generalization of the current, of the type 

^ ^^{A'ip ...ip ^ , (1.23) 

where the action of {. . .} on Lorentz indices produces the symmetric, traceless part of the tensor. In 
a general way, the matrix element of towers of twist-two operators appear in the operator product 
expansion. 



Using symmetries, we can write down the general form for the matrix element of the latter 
tower between states of unequal momenta. All possible Form Factors therefore arise, and for the 
vector operator, we obtain |120j 



[^1 



1=0 



n — 1 ] 



+ U{P') \^ % (P) Yl i?n,2.(t)A^^..A'^^«p'^^'+i...p'^-l} 

i=0 



+ Cn(i)^(l + (-ir)^n(P')^(^)A^^'A^^..A^'"-i> . (1.24) 

Time-reversal invariance imposes that only even powers in A appears. 

We project the result (1.24) on the light-front by multiplying by n^n^^ . . . n^„_;^ and define the 
off-forward parton distributions through the Form Factors 

n,n,,...n,^_APW-''"\P) = u{P')^+u{P) H^{i,t) + u{P')^^u[P) E^i^.t) (1.^5) 

with Hn{i,t) and En{£,,t) the moments of the helicity-conserving quark GPDs which are defined 
as 

= 5^ A„,2i(t)(-20'^ + C„(t)2(l + (-ir)(-20" , 

En{i,t) = j dxx"-is(x,e,t) 

jn-l j 

= ^ i?„,2i(t)(-2e)'^-a(i)^(i + (-i)")(-2e)" . (1.26) 

Basically it means that the moments of the GPDs are polynomials of the skewness variable ^ of 
order at most n. 



We have here presented the properties of the bifocal matrix elements defining the GPDs of the 
proton. The latter can be easily adapted to, e.g., the pion. 
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1.4 Transverse Momentum Dependent Distribution Functions 

All along this Chapter the quantities have been defined considering only collinear quarks. 

Let us now account for the transverse motion of quarks [27\ . The number of Parton Distributions 
at leading-twist increases from the three usual (i.e. number density q{x), helicity density Aq{x) 
and transversity ATq{x)) to eight distributions; 

fcT-dependent 



q, Aq, Atq , giT, ^il> ArA_ 

T-odd 

For the description of the quark content of the proton, the following quantity]^ is relevant, 

<^^J{x,kT) = |^^^^e^(-P^^--^^-^^-)(p,5|V^,(0)^^ , (1.27) 

depending on the light-cone fraction of the quark momentum, x = /p'^ and the transverse 
momentum component fcy. Using Lorentz invariance, hermiticity, and parity invariance one finds 
that up to leading order in 1/Q 



-if e jf^n'^ k'' S"^ 

^>(x, kr) = -< q{x, kr) +fiTix, kr) ^"'"^ ^ ^ + 9is{x, kr) 75 i+ 

+Arq{x^ kr) ia,^^,n%S^r + ^t) + h\{x, kr) ^^^^ | , 

(1.28) 

where S is the spin of the proton target. The quantities gis and h^^ are shorthand for 

gisix^kr) = XAq(x,kT)-\ 9iTix,kT) , 

hi{x,kT) = XhiLix,kT) + ^^^hiT{x,kT), (1.29) 

with M the mass of the proton, A = M S'^/p'^ the light-cone helicity, and St the transverse spin 
of the target hadron. 

The three distribution functions q{x, kx), Aq{x, kx), Axqix, kr) reduce to their analogs under 
integration over fcr, whereas the 5 others vanish if integrated over the transverse momentum. 

The two naively T-odd distribution functions /j^, are respectively the Sivers and the Boer- 
Mulders functions. Their role in the left-right asymmetries as well as their contribution to the 
understanding of the spin structure of hadrons will be highlighted in Chapters [6] & [7[ 



^It is given in the light-cone gauge = 0. 
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Pion Distribution Amplitude 



This Chapter is devoted to the pion Distribution Amphtude (DA). The DA (p is the probabihty 
amphtude for finding, in the pion, the valence quarks sufficiently near the light cone. Its behavior 
as well as its evolution in has been extensively studied in Refs. |45 |ll40j . We also study the pion 
Parton Distribution Function (PDF). Nowadays accurate experimental data allow for a comparison 
from the model calculations. 

Since there are a well-known quantities, the pion DA and PDF represent a perfect framework for 
illustrating our approach. Also, the properties of the pion DA and PDF are placed under scrutiny. 
We consider sum rules, support and soft pion theorem; the latter linking the DA to the PDF or 
to the Generalized Distribution Amplitudes. The discrete symmetries as well as isospin relation 
are used to relate the different pion DAs. We draw the way to QCD evolution of the DA and the 
PDF. We eventually give the expression of the electromagnetic Form Factor and its well-known 
asymptotic limit; namely the Brodsky-Lepage Form Factor. 



2.1 The Pion DA in the NJL Model 



In this Section, we illustrate our approach by performing the complete calculation of a simple 
quantity: the pion DA 4>{x). This will fix the conventions and allow us to give the parameters of 
the model. Also, the result of the calculation will tell us whether the model is appropriate for the 
calculation of such quantity. In particular, in the next Sections, we will see that the normalization 
condition as well as the support in x are respected. 



By definition, the pion DA is, as depicted in Fig. 2.1 



/ 



dz 



J{x-\)p+z 



(0|g 



i^^-T 9(0) \'^{p)) 



iV2f^(t) (a 



(2.1) 
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The DA has a support in x £ [0,1] and obeys the normahza- 
tion condition 

dx(j){x) = 1 . (2.2) 

' ' P 

The field theoretical approach we propose to use here has 
been first developed in Ref. |153j . The pion is described as a 
bound state in the sense of Bethe-Salpeter. In order to have 

an exact solution for the bound state, the kernel has to be ^'^ure 2.1: The pion DA. 

chosen carefully. The Nambu - Jona-Lasinio interaction being point-like, it leads to such a solution. 
Furthermore, as it is extensively explained in Appendix [Cj the NJL model respects all the required 
symmetries of the problem. The pions appear to be the three Goldstone modes coming from the 
chiral symmetry (dynamical) breaking. 

The shortcomings of the model are the absence of confinement and the non-renormalizability due 
to the point-like interaction. A regularization scheme should be chosen indeed. Also, the choice 
of this scheme defines the model per se. The NJL model completed with its regularization scheme 
is considered as an effective theory of QCD; even if built with quarks. Lorentz covariance and 
gauge invariance will assure the recovering of the properties of the parton distributions; so that the 



Pauli-Villars scheme is chosen that fulfills these conditions (see Section C.2). The regularization 
parameters, namely the cutoff A as well as the coupling strength G, are determined by calculating 
the pion decay constant /tt and the quark condensate. For f-j^ = 93 MeV and (uu) = (dd) = 
-(250MeV)^ we obtain 

A = 859 MeV , GA^ = 2.84 . (2.3) 

From the gap equation Eq. (C.9), we also fix the value of the constituent quark mass to 

m = 241 MeV . (2.4) 

This model for the pion is used in the description of the bound state. The Bethe-Salpeter (BS) 
amplitude for the bound state vr with four-momentum p is given by 

Xf}a{xi,X2;p) = {0\Tqp{xi)qa{x2)\Tf{p)) . (2.5) 

The NJL model contains a four-fermion interaction. Direct and exchange diagrams are related 
to each other by a Fierz transformation. Since they go like 1/Nc, the exchange diagrams are 
subleading and we may only consider direct diagrams. We also work in the ladder approximation, 
which consists in considering the iteration of the simplest closed loop with the kernel given by the 
model. The NJL interaction, giving rise to the kernel 

Vap^Sj{k,k';p) = 2iG{i-f5T^)s'y{i^5T^)ai3 , (2.6) 
renders the BS equation for the bound state easy to handle. 

Precisely, the integral equation generated by the chain of diagrams is the Bethe-Salpeter equation 
|112j for the bound state. Solving the Bethe-Salpeter equation in the ladder approximation with 
the kernel y(/c, A;';p) we find the expression for the Bethe-Salpeter (BS) amplitude 

X{k;p) =is(^k + ^)$ {k,p) iS (A: - , (2.7) 



^ See Eqs. ( |C.18[ |C.19| in Appendix [c] 
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with S (k) =^t —m. The quark-pion vertex function is 



(2. 



with the quark-pion couphng constant g-i^qg determined by the standard normahzation of the BS 
equation. 

In the theoretical scheme defined above, the pion DA is evaluated by developing the matrix 
element on the r.h.s. of Eq. (2.1) and rearrange it in such a way that we can identify the BS 
amplitudes xi^'^p)- Therefore the Eq. (2.1) becomes 

dz 



27r 

dz' 



(2'"2) 



We define the Fourier transform, Ref. |153j . 



Xp{xi,X2) 



-iPX 



d^k 
(2vr)4 



e-*'-x(A:,P) 



(2.9) 



(2.10) 



with the center of mass and relative coordinates X = nixi + ^2X2 and r = xi — X2 and where 
IJ'1,2 = mi,2/{mi + 7712). So that 



dz 



2tt 

{-'i'gTvqq)Nc 



d^k 



e-'^-'Tx[x{k,v)i^^T- 



dz 



Jp+{x-\)z 



^^e-^^-^ir[S{k + |)z75r-^5(fc - |) ^ tst"] 



(2.11) 



The symbol Tr at the first line is the trace over the Dirac, isospin, flavor and color spaces. Going 
one step further in the calculation, we perform the trace over isospin for a 7r+. By making the 
change of variable k ^ —k + p/2, we find 



V2/,<A"^(x) 

= -i{-4.m)p+ {-ig^gg)NcV2 



d^k 6{p+{x-l) + k+) 

(27r)4 {{-kY -m? + ie) {{-k + p^ - + ie) 



(2.12) 



The 5-function imposes a decomposition on light-cone components in order to perform the d'^k 
integral 



d'^k 



1 



(27r)4 (27r)2 J (27r)2 



d'k 



dk- 



dk^ 



The integral over dk~^ is trivial: the 5-function imposes the +-component of the momentum carried 
by the outgoing quark. The remaining integral over d^/c-*- in Eq. (2.12) is divergent. Since it 
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' m„=140MeV 

m„=0 MeV 


/ 





Figure 2.2: The pion DA in NJL. The dotted blue line represents the DA for the physical value of the pion mass 
while the plain orange line represents the chiral limit. 

has been obtained with the NJL model interaction, the integral over cPk-^ has to be consistently 
regularized using the Pauli-Villars regularization, which is part of the model. 



The result Eq. (2.12 ) is proportional to the two-propagator integral /2,p(x, ^) defined by Eq. (E.ll ) 
i\/2/^0"^(x) = -4mig.,qgNcV2i2,p{x,0) . (2.13) 



In doing so, the normalization condition of the DA Eq. (2.2) is recovered for any value of the pion 
mass. 



The result in the NJL model is shown on Fig. 2.2 for both the physical pion mass and the chiral 
limit. In the latter case, the pion DA is a constant which is equal to 1. This result is obtained 



by plugging as given in Eq. (C.45) into Eq. (2.13). In the physical case, the DA deviates by 
less than 3% from 1. The distribution amplitude does manifestly not vanish at the end points 
for none of the pion masses. This result is obviously characteristic of the model within which the 
calculation has been performed. It is not a QCD result, and in fact, QCD evolution, as will be 
shown later on, will push the DA to zero at the end points. 



2.2 Symmetries 

In the previous paragraph, we have considered the case of a vr^. Using discrete symmetries we can 
relate, in a model-independent way, the vr^ DA to the, e.g., vr^ DA. 

For a generic particle V with momentum p, spin projection Sz and charge Qi, we write the 
particle state |P(p, s^, Qj)). Parity, time reversal invariance and charge conjugation act on this 
generic state in the following way 

P\'Pip,s„Q^)) = rjp\V{-p,Sz,Qi)) , 
T\V{p,s,,Q,)) = r]T{-ly-'^\V{-p,-Sz,Q^)) , 

C\r{p,Sz,Q^)) = Vc\r{p,sz,-Qi)) , (2.14) 



^Mind the scale in Fig.p^ 



2.3 The Pion Parton Distribution Function 
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with rjp, rj'T, rjc the respective phases. For the vr"'", which is J^^ = are r]p = t]T = — !> = 1- 



In the hght-front, the T symmetry as defined by the second expression of Eq. (B.l) changes x 



to —X and vice versa. It is therefore useful to define the V symmetry has the combination P^T, 
what gives 



V\V{p,s,,Qi)) = rjv{-iy-''\V{p,-s,,Qi)) 



(2.15) 



with rjy = rjp rjT- 

In Appendix [B] are developed in details the symmetry relations for any type of distributions. 
We here relate such relations for the pion DA as given by Eq. (2.1). Using the relations Eq. (B.l), 
Time reversal invariance leads to 



charge conjugation to 



and CPT to 



(i-x) 



'{l-x) ; 



6^ (x) . 



(2.16) 



(2.17) 



(2.18) 



As perfectly illustrated on Fig. 2.2 the pion DA is symmetric around x = 1/2 so that 



(/)(1 — x) as required by isospin: it is easily seen from Eq. (2.26) that it is implied by the fact that 
the quarks carry the same constituent mass. By Eq. (2.17), it means that there is no difference 



between the DA for a vr^ or a 7r~. This also corresponds to the isospin relation: the distribution 
amplitude for a u-quark is related by the change x — >• 1 — a; to the amplitude for a d-quark. 
Those relations remain unchanged under evolution. 



2.3 The Pion Parton Distribution Function 



The nomenclature of distribution functions includes a huge variety of transitions. As overviewed 
in the Introduction, the names are mostly related to the nature of the final and initial states. For 
instance, the DA is a transition from a physical state to the vacuum. Also, by definition, the pion 
Parton Distribution Function is the probability density to find a quark carrying a fraction x of the 
parent pion's momentum, i.e. 



/ 



dz 



^i{x-\)z p+/„+ 



|vr+(p)) 



1 



q[x) 



(2.19) 



as illustrated on Fig. 2.3 In the NJL model we obtain 



qq 



1 



^/2,p(x,0) + m^xl3^^^(x,0,0) 



(2.20) 



with the 2- and 3-propagator integrals defined Eqs. (E.IO, E.ll E.13). In the chiral limit, q{x) is 



equal to 1. The NJL model has been applied to the study of pion Parton Distribution in different 
occasions |76 | [77 1 ri79j . More elaborated studies of pion Parton Distribution have been performed 
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in the Instanton Liquid Model |14| . and in lattice calculation inspired nonlocal Lagrangian mod- 
els I156j , which confirm that the result obtained in the NJL model for the Parton Distribution 
is a good approximation. 



A low-energy theorem based on PCAC]^ links the pion DA and the distributions of two physical 
pion states. 

Following the same steps as in Ref. [86], we hereafter derive this relation between the pion DA (j){x) 
and the pion PDF q{x)\pi^ which leads to the particularly interesting result 



where x means m-,^ = MeV for the pion under scrutiny. 



(2.21) 



P 




P 



Let us take the matrix element of a pion-pion transi- 
tion with momentum transfer pi —p2', we call it f""^. 
Such a matrix element obeys the following isospin 
decomposition |17Uj 



Figure 2.3: The pion Parton Distribution Function. It is related to the U— and (i— quarks distributions 
The orange blobs represent the pions. ju _^ jd^ consider the isovector contribution 

of such a matrix element 



{x)h — q{0) 



vr'^(pi) 



f2.22) 



Using the Lehmann-Symanzik-Zimmermann reduction formula we can write 

{Ap2Mx)h^qm7^'^ip^)) 

= i [ d^xe-^P^-^ (al + ml) {0\T {7r\y)qix)h ^ g(0)l K'^(pi)) , 



where tt^ (x) is the pion interpolating field and is given by PCAC, Eq. ( C.42 ). We now integrate twice 
by parts and use the relation for the derivative of the Time-ordered product dt'T{ip{x')(p^ (x)} = 
Tmy,{x'))ipHx)} + 6{t' - t) Mx'),y,Hx)]. 



Taking the soft pion limit, i.e. P2 — )• 0, the matrix element (2.22) becomes 

hm {7r\p2Mx)n - g(0)|7r'^(pi)) 
2 



5iyo-xo){OWiy),qix)]n - 9(0)|7r"(pi)) 



+6iyo)mx)h - [A'\y),qmKiPi)) 



(2.23) 



'Partial Conservation of the Axial Current; see Section C.4 in Appendix [c| 



2.4 Towards the Support Problem 
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The axial current gives rise to the axial charge = J d'^x 5{xo)Aq{x) with the equal-time com- 
mutators obeying [Qf, q{x)] = i/2 r'^" 75 q{x). 
The expression (2.23) becomes 



lim {7r\p2Mx)n - q{0)K{pi)) 

p'i^o 2 



1 



(0|g(x) ^/2^75n y g(0)|^»(pi)) 

+(0|g(x)n y V2y75g(0)|vr'^(pi)) 



I € 



3bc 



(0|g(x)nV2-75(?(0)K'^(pi)) , 

Jn ^ 



(2.24) 



where we have used the algebra [t'^/2 , t^/2] 



^ e 



r72. 



In order to go to the pion Distribution Functions, we Fourier transform Eq. (2.22) for quarks 
at a light-light separation and replace by p and make the transition diagonal in isospin and in 
momentum. In the soft pion limit, the isovector pion PDF is exactly the pion DA 



lim 



dz 



dz 



2 " V2 

7+75 r 



Ki+(p)> 



(2.25) 



what corresponds to the definition of the pion DA, Eq. (2.1). This results takes all its physical 
meaning in the chiral limit m-,^^ = niT^^ = 0, where both pions are chosen to have zero mass by 



consistency. The PDF actually has to be diagonal in state. If one of the pion is massless {p^2 ~^ 
implies that p^^ = 0), the other pion has to be massless as well. We find the expected result of 
Eq. pliri). 



2.4 Towards the Support Problem 

We now focus on the support property that is easy to illustrate in the context of 2-propagator 
integrals. Nevertheless the discussion that follows is extended to 3-propagator integrals. 



The support in x is determined by the pole structure in k of the function given in Eq. (2.13). 



Developping this integral, Eq. (E.ll), leads to 
i V2U<Pix) 



1 f d'^k± 
-i'imp^ ig^qqNcV2-^—^ J dk 

1 



X) yK 2p+{l-x) 2p+{l-x) \ X\r) K 2p+x + 2p+x 



(2.26) 
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The integral over dk~ is performed in the complex plane making use of the Cauchy's theorem. 
Would X g/ [0, 1], both the poles would be on the same side of the real axis and the integral would 
be identically zero. The definition Eq. (2.1) therefore leads to the required support in x, i.e. 
X G [0,1]. 

It is true as long as the regularization scheme that we use is covariant. This condition is 
necessary but not sufficient: the regularization procedure should preserve the symmetries of the 
light-cone instead. Regularizing Eq. (2.26) through a 3-momentum cutoff would break covariance, 
spoiling the required support. However, it is obvious that a cutoff of the type A = (A~, A_l), with 
A~ < oo prevent us from closing the contour in the complex plane of k~ . This leads to support 
corrections of the order of 1/ A~ . Therefore the loss of light-cone invariance here produces the loss 
of the correct support. 



Many calculations of PDFs show problems with the support properties. In particular, the sup- 
port problem occurs in calculations performed in non relativistic (NR) Constituent Quark Models 
(CQM), e.g. Ref. |200j . In such NR calculations, the light-cone coordinates are defined from the 
3-momentum and the energy, = m? + k'^, so that the dependence on k~ is not the genuine one. 
A slight support violation is observed. In the same Reference |200j . the inclusion of Final State 
Interactions is proposed to restore translational invariance and therefore the support, which in turn 
will restore the number of particle conservation. A trivial solution to the support problem would 
be to change the normalization in order to rescale the support to compensate for the problem. In 
the same line, we also mention that a slight violation of the support is found in the MIT bag model 
e.g. Ref. [188] . 



Let us now comment on the role of the interaction on the support problem, staying in the 
description of the hadron as a bound state in the sense of Bethe-Salpeter. 

The vertex function ^{k,p) defines the BS amplitudes, which are afterwards plugged into the 
distribution in the BS approach (2.9). The former is defined within a model with its given sym- 
metries; whereas the latter requires covariance on the light-cone. It is important to realize that, 
in order to preserve the support property, one has to go to the light-front without breaking its 
covariance |156j . a constraint that is respected by NJL model. 

The BS equation for the pion has been studied using different vertices, e.g., [571 ISSl 11951 11961 
11971 1204j . In the five first References, use has been made of a reduction of the BS amplitudes to 
light-cone wave functions by projecting on the light-cone. In Refs. [57', '58], the light-front BS vertex 
functions were replaced by wave functions obtained in a light-front CQM. In Refs. |195l 1196] 1197] . 
the scalar Wick-Cutkosky model was adopted and successfully applied to the calculation of scalar 
meson GPDs. The correct support has been found. It is however not the case in Ref. \203\ 1204] 
where the calculation is performed in a covariant CQM, i.e., the Bonn model. Once more, the 
problem is initially non covariant but the authors attempted to restore Lorentz invariance by 
applying a boost; attempt which was nonetheless either too schematic or incomplete. In Ref. [202J 
a detailed study of the problem suggests that the dependence upon the energy, k^ , of the vertex is 
responsible for the introduction of additional dynamical poles. As a counter example we mention 
the non-local Lagrangian |155| used to calculate PDFs in Ref. [1551 1156] . The non- locality of the 
interaction obviously leads to a ^"-dependent vertex, whose definition is chosen consistently with 
the definition of the mass m{p), given there by lattice QCD. The problem is totally covariant and 
the support is recovered in the physical cases. Since the calculation in non-local models is performed 
in the Euclidean space, it is the Wick rotation that, in this case, might spoil the support property. 



2.5 QCD Evolution 
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It is therefore tempting to explain the fulfillment of the support property by the conservation 
of covariance on the light-cone. 

2.5 QCD Evolution 

In a previous Section, we have proven that the pion DA and the PDF coincide in the chiral limit. 
Nevertheless this result based on PCAC must be strongly broken by QCD evolution when we move 
from low to high . We will now discuss how evolution changes the PDFs and the DAs. Our first 
objective will be to determine the scale of the model, which can be fixed from our knowledge of 
PDFs. 

The Parton Distributions here obtained are valid at a low scale where scaling still holds and our 
model is still defined. However, one would like to link the results of a model calculation with QCD. 
In other words, one would like to know the corrections when going to higher energy. The start of 
the scaling violations is given by the so-called evolution equations, which schematically describe 
the dependence of the distribution functions on 

7 

The pioneer works on the behavior of Parton Distribution Functions (PDFs) with were done 
by Dokshitzer [H^, Gribov and Lipatov jlU4[ I141| as well as by Altarelh and Parisi The QCD 
evolution of the Distribution Amplitudes has been studied by Efremov and Radyushkin [88j as well 
as by Brodsky and Lepage |14U| . This gave rise to the so-called ERBL evolution equations. 

The role of the model calculations is to provide for the initial conditions to the evolution 
equations. 

The initial scale from the NJL calculation 

The scale of validity of the model calculation is still to be found. One needs to fix the value of Qq for 
which the quark distributions obtained in the NJL model are considered to be a good approximation 
of the QCD quark distributions. For the purpose of evolving the parton distributions, we use the 
code of Freund and McDermott [94] . 

The only information we have at hand is the momentum sum rule. Obviously momentum con- 
servation holds and we know that the sum of the momentum fractions carried by each constituent, 
namely the valence quarks (f), the gluons (^g) and the sea quarks (s), i.e. the second moment of 
the PDF /(x,Q) 

rdxxf{x,Q) (2.27) 
Jo 

has to be conserved. In the models we are interested in, only valence quarks are taken into account, 
so that, 

{x),{Qo) = l , {x)s{Qo) + {x)g{Qo) = ■ (2.28) 

To leading order (LO), scaling violations are manifested in a logarithmical dependence, of the 
PDFs, on the scale through the running coupling constant 
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with 13 = 11/3 Nc — 2/3 Nf and A is the scale of QCD. In order to find the value of Qq, we need 
to fix the value of A consistently with the evolution code; i.e. we choose A = 0.174 GeV. Following 
the NLO evolution, one is to use A = 0.246 GeV. 

Knowing that the momentum fraction of each valence quark at Q = 2 GeV is 0.235 [192J, we 
fix the initial point of the evolution in such a way that the evolution of the second moment of the 
pion Parton Distribution reproduces this result. This condition is fulfilled at a rather low value, 
i.e. 



Qo = 0.29GeV 
Qo = 0.43GeV 



for the LO evolution ; 
for the NLO evolution . 



(2.30) 




The first result that we obtained is that the PDF 
resulting from the NLO evolution is basically unim- 
proved with respect to the one evolved to LO. In 
order to illustrate the latter statement, we have de- 
picted the pion PD evolved at both LO and NLO in 



Fig. 2.4, The solid (black) line of the Figure corre- 



sponds to the LO evolution of the NJL model pre- 
diction while the dashed (red) line represents the 
NLO. Both evolved results are compared to the ex- 
perimental data |65| . The agreement with the data 
is not better in one or the other case. The effect of 
the NLO evolution is compensated in the LO evolu- 
tion going to a lower value of Qoi a result that has 
already been noticed in proton Parton Distributions 
|200| . It is therefore obvious that, for numerical rea- 
sons, we will prefer to evolve the distributions to LO. 

The method used here is also applied by the au- 
thors of Ref. |103j . There exists other ways of fixing 
using different data or comparing with lattice data. They are reviewed in the latter 



Figure 2.4: Pion parton distribution. See text. 



the scale. 
Reference. 



The QCD evolution of the Pion DA 

The Distribution Amplitudes have a logarithmic dependence in Q which is completely determined 
by the QCD evolution equations derived in Ref. | 138t 1139] I140j . The QCD evolution equations for 
the distribution amplitudes can be expressed, to leading order, in terms of Gegenbauer polynomials 

cf>{x,Q) = x(l-x)^ a„C^/2(2x-l) (inj^ ) , (2-31) 



n=0 

even 



where only even n contribute since (j){x, Q) = (j){l — x, Q) is required by isospin. The anomalous 
dimensions are 
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where Cp = {N^ - l)/(2A^c) and /3 = ll/3Nc- 2/3 Nf. 

The evolution equations would be completely known if the coefficients a„ were determined. These 
coefficients can be determined from the initial conditions which are provided by the model calcu- 
lation. The evaluation of the DA in the model provides us for the 0(xi,(5o)- By using this initial 
distribution together with the orthogonality relations for the Gegenbauer polynomials, 

a, In^ =4 / dxCl/\2x-l)^{x,Qo) , n = even , (2.33) 

V A^y (n + l)(n + 2) Jq 

one might be able to find the a„ coefficients, namely, to determined the evolution equations. 




Figure 2.5: The pion DA in the NJL model evolved to the scale Q — 2 GeV (red curve) and compared to E971 di-jet 
measurement [3] after proper normalization of the data Ref. i45j. The dashed line represent the asymptotic DA. 



All the model calculations, after evolution, should agree with the theoretical asymptotic behav- 
ior of the DA, that is given by Eq. (2.31) with 
only the first term, i.e. n = 0, survives 



oo. Since 70 = and 7„ > for positive n. 



4>{x,Q) 



(2.34) 



Using the orthogonality relations for the Gegenbauer polynomials Eq. (2.33) and the renormaliza- 
tion of the charged weak current, it is found in Ref. [140j that ao = 6. We can find the other 
coefficients through the initial conditions given by our result in the NJL model. As mentioned in 



Section 2.1, in the chiral limit, the pion DA is a constant normalized to 1 



<Pix,Qo) = 1 , (2.35) 
what holds at the scale of the model, Qq. It is then easy to apply the above-described formalism 



and determine the an for our model calculation. Plugging Eq. (2.35) into the relation (2.33), we 
fincS 



4>ix,Q) 



Axil 



00 

E 

n=0 



2n + 3 



[n + l)(n + 2] 



C'J\2x-l) 



/ln(QVA2)Y 
Vln(Q2/A2); 



(2.36) 



*/g dx CfJ'^{2x — 1) = 1 for n = even. 
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At a scale Qo, the series of the previous result should give 1. However, it must be realized that 
this result is reached for n — t- oo which is numerically not feasible. At a low scale Qo < Q ^ ^ 
GeV, where changes due to evolution are more important, one should include a high number of 



terms of the series. As proven by Fig. 2.6, the desired result is reached by increasing the n value. 



It nevertheless takes an infinite number of terms to compensate for going to 1 at the end-points. 
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Figure 2.6: The truncated series (2.361 for the pion DA in the 
NJL model at Qo- 




Q=1 GeV 
Q=2 GeV 



Q=20 GeV 
Q=40 GeV 



Asymptotic 



Figure 2.7; The pion DA in the NJL model evolved to the scales 
Q = 1, 2, 20&40 GeV as well as the asymptotic form. 



On the other hand, once QCD is switched 
on, the end-point problems for the pion DA 
calculated in the NJL model in the chiral 
limit is solved, showing the need for QCD 
evolution. 



The pion DA evolved to a scale of Q 



1 GeV is shown in Fig. |2.5| for the NJL 
model with Pauli-Villars regularization. It 
is compared to the data of the E791 mea- 
surement [3j once properly normalized |45j . 
The asymptotic DA is also plotted that 
seems to be in a better agreement with the 
data. 



In our version of the NJL model, the 
scale of the model is found to be as low as 



At 



290 MeV at LO, see Eq. ( |2.30 ). 
high Q, one could take into account only 
the 4th first terms of the series. This ap- 
proximation is justified by the fact that 
the coefficients o„ decrease with increas- 
ing n and by the asymptotic behaviour in 
Q. However, for a better resolution and 
smoothness of the plots, we choose to trun- 
cate the series at n = 512. 
In order to illustrate the result obtained in 



Eq. (2.36), we plot the series in Fig. 2.7 for 



different, high values of Q. For Q higher 
than 1 GeV (pink curve), the effect of evo- 
lution starts being less important. The 
green, blue and red lines represent the LO 
order evolved at, respectively, Q = 2 GeV, 
Q = 20 GeV and Q = 40 GeV. The result changes slowly showing that no big improvement to reach 
the asymptotic value can be done by this model calculation at least at LO evolution. One can how- 
ever easily check that the analytical asymptotic behavior for the series. (2.36) is Q) = 6x(l— x), 
as required. 



2.6 End-points 
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2.6 End-points 

Although its role seems to be only that it provides for the initial conditions, a model calculation is 
a calculation that has a physical meaning. As a consequence, the results in the model are supposed 
to obey all the properties of the initial object. Turning the problem over: the results of model 
calculations are not complete until evolution is performed and most of the properties of Parton 
Distributions that were respected at the quark-model scale, e.g. isospin relations, will still be after 
evolution. On the other hand, we will show that other intrinsic properties can be restored. 



In particular, one could evoke the problem of the distribution functions at the end-points. 
Effectively, in the free parton model, PDFs are expected to vanish at the boundaries, i.e. x = 0, 1. 
For on-shell particles, at x = 1 one of the quarks carries the whole (+)-component of the momentum 
and the other quark has zero (+)-component of the momentum, i.e. on-shellness implies an infinite 
(— )-component of the momentum. Since the (— )-component of the momentum appears in the 
propagator, the PDF vanishes. Also, this end-point behavior has been analytically predicted in 
perturbative QCD, e.g. [30] . 

It is however not what happens for the results in lots of models such as the NJL model illustrated 
in, e.g.. Fig. 2.2 This unwelcome behavior comes from the misunderstanding of the Nature of the 
quarks beyond the naive parton model. It is effectively unclear how the constituent quarks of 
the quark models can be linked to the current quarks of QCD. We are willing phenomenological 
evidences to shed some light on the constituent density distributions. The scenario suggested by 
Altarelli, Cabibbo, Maiani and Petronzio (ACMP) ^ pictures the constituent quark as a complex 
system of point-like partons. This ansatz is so far one of the only for those distributions. A variation 
of this scenario has been applied in, e.g., Refs. 11561 11861 1187j . showing a huge improvement. 

As for the calculation in the NJL model, the end-point problem shows the limits of validity 
of our model calculation. The simple considerations given above do not hold for the bound-state 
quarks in the NJL model; and as a matter of fact we cannot expect the Parton Distributions to 
vanish at the boundary. This behavior has been shown to arise in a regularization-independent 
way as they appear in both the results of Refs. [35 l I153j . 

Turning our attention to the chiral limit, the question could also be formulated as follows: what 
is, in the NJL model, the probability density of finding -massless(current)- quarks with momentum 
X inside a massless pion? The quarks can, a priori, not carry a particular fraction of the massless 
pion momentum, so that the probability density does not depend on x, explaining at the same 
time the end-point behavior and the result q{x) = 1. This problem is the opposite of the zero 
bounding energy problem, where the mass of the pion is chosen to be twice that of the constituent 
quarks. In that case, the distribution is rather peaked in x = 1/2 as shown by a quick analysis of 
the kinematics. 

As for the treatment of the quarks in the calculation given here, our choice has been to maintain 
consistency with the model in itself. The relation between the current quarks of the NJL model, 
with mass mo, and the constituent quarks, carrying a mass m, is the one given by the Bogoliubov- 
Valatin transformation which is proper to the model. |^ Therefore an additional convolution of the 
type ACMP is not necessary. 

We conclude by highlighting the fact that the obtained values about q{x) = 1 in the vicinity 
of X = 0, 1 enable a good reproduction of the evolved PDF at the end-points, see Fig. 2.4 On the 



^See Appendix [c] 
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Figure 2.8: The Feynman diagram representing the elastic scattering. 



other hand, similar considerations for the pion DA, e.g. Fig. 2.5, indicate that the obtained DA is 
too high at the end-points. 



2.7 The Pion Electromagnetic Form Factor 

In this Chapter we have examined the meson Distribution Amplitudes as an illustration of the 
method we will use throughout this thesis. The DAs contribute to both elastic and deep inelastic 
processes, what made the example interesting. In particular, they can contribute to the elastic scat- 
tering process Fig. |2.8[ It is, however, not the simplest situation for the study of elastic processes: 
The simplest objects that allow for the study of hadron properties are not the meson DAs but the 
elastic Form Factors. In effect, the Form Factor is the amplitude for the, e.g., pion to absorb large 
transverse momentum while remaining intact. In this Section, we aim to link both quantities and 
to discuss the asymptotic - and model-independent - limits of the pion electromagnetic Form Factor. 

For a charged pion vr"*" the electromagnetic Form Factor i^V is defined as 
(^+(p')l-/^'"^(0)|vr+(p)) = {p + p'),eF^{q^) , 
with = {p' — p)'^', and its charge radius by 



-6- 



dq^ 



We evaluate the diagram on Fig. 2.8 in the NJL model, with both incoming and outgoing particles 



being pions as well as its crossed version, applying the technique of the Section 2.1 We find 



nqq 



■ P ■ P 



- {-ml hip) -p-p hip - p) + ml hiP,p)) 



(2.37) 



with the 2- and 3-propagator integrals defined in Appendix |E| The charge radius evaluated in the 
NJL model is found to be {r'^)-^ = 0.31 fm^, which has to be compared to the experimental value 
of 0.44 fm^. It must be realized that our model does not include the vector mesons, which also 
contribute to the charge radius of the pion, see e.g. [57] . 



Here, we have used an effective theory of QCD in order to calculate the electromagnetic quantity. 
The expression (2.37) is a useful check for the sum rules for GPDs and will be used in the next 
Chapter. 



2.7 The Pion Electromagnetic Form Factor 
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The Asymptotic Limit: the Brodsky-Lepage Form Factor 

The electromagnetic Form Factor can also be seen as the convolution of three probability ampli- 
tudes. Basically, those probability amplitudes are the probability of finding a 2-quark valence state 
in the incoming pion, i.e. the pion DA (p; a hard amplitude Th for this 2-quark-state to scatter with 
the photon and, therefore, to produce 2 quarks which can reform in similar 2-quark state |140j . 
The latter is described by a pion DA (j)*^ 



dxi dyi5{l - T.iXi)5{l - T,iyi) (j)*{yi,Q)THixi,yi,Q) (j){xi,Q) 



(2.38) 



where Th contains all the 2-particle irreducible amplitude for "y*+qq ^ qq- This approach is similar 
to the previous one in that the pion DA is well-described by the Bethe-Salpeter amplitude on the 
light-front. However, the pion DA calculated within this formalism does not explicitly depend on 
the scale unless we switch on QCD. Now, the effects of QCD evolution need to be taken into account 



into the convolution (2.40), so that we need the evolved pion DA. 



Using the light-cone perturbation theory together with the light-cone gauge, Brodsky and Lep- 



age (1551 1139j obtained the evolution equations of Section 2.5 Those equation can in turn be 
applied to the calculation of the Form Factors. The hard amplitude gives the behavior in powers 
of Q^. For instance, plugging Eq. (2.31) with Th evaluated in the same light-cone perturbation 



theory, they found, to leading logarithmic corrections 



FAQ^ 



n=0 



log 



A2 



[l + 0{as{Q^),m/Q)] 



(2.39) 



which shows the most important dynamical feature of the Form Factors, namely, their power-law 



fall-off. The asymptotic limit is given by Eq. (2.34) with known qq 

Q2^^ 4^Cpa,(Q2) 9/2 



F^iQ' 



N, 

WnasiQ )- 2 



'Q2 



(2.40) 



with Nc = 3. 
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Pion Generalized Parton Distributions 



The adjective generalized in the name Generahzed Parton Distributions (GPDs) appeared in the 
early nomenclature of distribution functions \in\ IllSj . Their introduction through the factor- 
ization theorems for deep exclusive processes demonstrated the wide applicability of the diagonal 
transitions obtained with inclusive processes. Deep exclusive processes have also been discussed in 
terms of off-forward parton distributions jl73j for they are non-diagonal in momentum, i.e. there 
is a momentum transfer between the initial and final states. The off-forward distributions are 
univocally related to the GPDs. The initial and final states we will nevertheless consider as corre- 
sponding to the same quantum numbers when referring to Generalized Parton Distributions. Such 
distributions, with states not corresponding to the same particles, are called Transition GPDs or 
Transition Distribution Amplitudes. The study of the latter quantities is actually the main topic 
of this thesis and will be fully developed in Chapter |4j 

In Chapter [l] we have introduced the Generalized Parton Distributions for the nucleons. Never- 
theless such bilocal matrix elements can also be defined for mesonic states, e.g. pions. Let us now 
focus on the study of the pion GPDs applying the rules developed in the NJL model which have 
been given in the previous Chapter about the Pion Distribution Amplitude. 

The set of distributions proposed for the proton in Refs. [79l |82] is adapted to the pion case 
including both the chiral-even and the chiral-odd GPDs. The tensorial structure for the 7*7r — )• ^tt 
process allows a single twist-2 chiral-even GPD. In the light-front, the tensorial decomposition of 
this GPD is found, using the recipe of Ref. \15\ I79| . from the vector current 

= -l^H^{x,i,t){P + P')-n=\H^{x,i,t) , (3.1) 
with q{x) the quark field and where = (P' + P)^. On the other hand, the axial current gives 



az 
~2tt' 



vrfP') 
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rise to a twist-3 GPD, involving the 7"'"75 current. This distribution function will not be studied 
in this thesis. The effect of twist-3 distributions is small but, for nucleons, their inclusion is useful 
for the restoration of electromagnetic gauge invariance in DVCS, e.g. Ref. [15]. In the pion case, 
the twist-3 GPD can be associated to a Single-Spin Asymmetry jl2j . 



Taking into account chiral-odd, namely helicity-flip, GPDs makes the number of twist-2 pion 
GPDs increase. This case will be overviewed in the Future Perspectives in Chapter [7] 

In the literature, different approaches of the non-diagonality -in momentum- have been proposed 
for the hadronic matrix element 

7T{P')\q[-'-)r^q[^)\7,{P)) . (3.2) 

One bases its analysis of the longitudinal asymmetry on the kinematics of the process. Since the 
dependence upon the skewness variable is hence explicit, this approach leads to the so-called skewed 
parton distributions [in\ 1118] \17'6\ [T7^ . Another approach is to write the spectral decomposition 
of the matrix element of the bilocal operator at a light-like distance without kinematical assump- 
tions on the skewness variable, allowing one to make statements about this very dependence. This 
second approach leads to the so-called Double Distributions (DD), which have been proposed in 
Refs. [TiHl[T73l[T7i] . 



The Nambu - Jona-Lasinio model is an appropriate choice in the context of model calculation of 
both the pion skewed PDs and DDs. As above-mentioned the pion can be treated as a bound-state 
in a fully covariant manner following the Bethe-Salpeter approach inside the NJL model. Given 
the tools exposed in Chapter [2] we are now able to extend the application of the model for the pion 
to the definition of GPDs. So, in this Chapter we will aim to analyze the twist-2 pion skewed PD 
following the lines of Ref. |153] and link it to the DD approach [45j. In particular, semi-kinematical 
approaches deriving from the original DD as well as the resonance exchange contribution will be 
investigated. The QCD evolution already described in the context of Distribution Amplitudes will 
be overviewed. 



3.1 Kinematics for Skewed Parton Distributions 

We define the kinematics of the GPDs following the conventions of Ref. |107j . As it has been 
explained in Chapter [T| the momenta in deep processes are expressed in light-cone coordinates. 

The GPDs depend on three kinematical variables: the Fourier transform variable x, the mo- 
mentum transfer t and the longitudinal momentum asymmetry, also called skewness variable, ^. 



^ The tensor structure of the twist-3 GPD is given by 



2p 



o H {x,^,t) 



where by _L we understand the component 1 or 2. 



3.2 The Chiral-Even Pion GPD 
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The latter is defined, for an incoming momentum P and outgoing one P' , by 

P+ _ p'+ 



^ P+ + P'+ ' 

The momentum P of the incoming pion and the momentum P' of the outgoing pion are therefore 



defined, see Fig. 3.1 



A-i-fJ- 1 t _ 

P^' = + —+—AM'^n^' 



with 



2 2 
A-L^ 1 + f 



(3.3) 



1 f 2 A^2x i 



and where we have used the on-shellness conditions P^ = P'^ = m^. is related to the momentum 
transfer t by = t, and defined as follows 



A^ = (p' - py 



-4 ^2 ^2 _ ^±2 



(3.4) 



{P + P 



M2 



From the relation (|3.4l) and the fact that A-'"^ is a positive quantity, we can deduce the interval of 



kinematically allowed values for the skewness variable ^ 



< C < 



2\/m2 



(3.5) 




/ 



Figure 3.1: u and d-quark contributing diagrams to the GPD in the NJL model. 



3.2 The Chiral-Even Pion GPD 



The calculation of H'^ in the NJL model has been performed in Ref. |153] . In this Section, we 
comment these results. 



38 



Pion Generalized Parton Distributions 



For the tt"*", the two contributing diagrams are depicted in Fig. 3.1 The bilocal u-quark current 
in the BS approach is given by 

{tt^{P')\u{x') 7'^u(x)|7r^'(P)) 

= J d4x2Tr{xi„(x',X2)^(l + r3)7^[xM^,^2)(^?(')-m2)]} 

+ yd^xiTr{x^p,(xi,x)(ii W-mi)x^p(xi,x')^(l + r3)7^} , (3.6) 

the j index selects the isospin: j = 7r+,7r~ or vr*^. Here the projector ^(1 + t^) selects the contri- 
butions from the li-quark present in Eq. (3.6). The two terms on the r.h.s. correspond each one to 
one of the 2 components of the bound-state, namely the u or d-quark, as active quark. |^ 

We now set the quark masses equal to the constituent quark mass, i.e. mi = m2 = m. For the 
first diagram of Fig. |3.1[ namely the active u-quark, the variables become pi = P — k and P2 = k 
while for the active d-quark the variables are pi = —P + k and p2 = —k. That is, due to isospin 
symmetry also, we can relate both diagrams 



(3.7) 



We can then concentrate on the calculation of the u-quark distribution. The GPD (3.1) becomes 



I ^4 ^[^-^+'f)^'{xpUP'-k)lil + rs)^+x'^\P-k){-^-m)} . 



(3.8) 



The BS amplitudes are replaced by their expression in terms of Feynman propagators and the 



pion-quark vertex function Eq. (C.20). 



After having performed the traces and integrating over k~^, we can analyze the poles in k 



similarly to the case of 2-propagator integrals in Section 2.4, In the present case it is the interplay 
between 3 poles that defines the support property: the integral is non-vanishing if there are poles 
on both sides. Thus, two kinematical regions are found that correspond to a GPD for the li-quark 
of supports X £ 1] and x £ From Fig. 1.5, it is easily seen that the u-quark GPD con- 

tributes to the region of emission and reabsorption of a quark as well as to the region of emission 
of a quark- antiquark pair. 



We obtain a combination of 3-propagator, as it might be expected, as well as 2-propagator 
integrals; 



qq 



+ (m2x+^(l-x))/3«-(x,e,t) 



(3.9) 



We refer to Section 2.1 for the details of the technique. 
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+ . . . 



Figure 3.3: Effective interaction in the Random Phase Approximation where only the direct terms are con- 
sidered. 

with the Hght-front integrals given in Appendix [E) 



For the pion GPD, one should take into account the isoscalar meson-exchange in the t-channel: 
in the ERBL region, the two quarks can couple with quantum numbers corresponding to the a- 



meson. Since the NJL model describes the pion as well as the cr- meson, the diagrams 3.1 may already 
contain this particular exchange. Moreover, in the region — ^ < a; < a new diagram Fig. 3.2 should 
be considered as the model's expression of this c-exchange. This meson-exchange, which is an 
explicit manifestation of chiral symmetry, can be calculated in the Random Phase Approximation; 



2iG + 2iG{-iUs{k'^))2iG + ... 
2iG 



1 - 2Gn,(A:2) 



(3.10) 



where Ils{k ) is the scalar proper polarization given in Eq. (C.35). Notice that the expression 



of the interaction is similar to the one obtained with the Bethe-Salpeter equation in the ladder 



approximation Eq. (C.31). The central part of Fig. 3.2 is given exactly by Eq. (3.10). 




The second contribution to the GPD for the u-quark 
is then 



(3.11) 



with 



Git) 



(-4m^ 



P.P'l3{m,P,P')-l2{m,-A) 
m^/2(m, m^) -|- (4m2 — t) him, A) 



Figure 3.2: u-quark diagram contributing to 
the GPD in the NJL model coming from the 
cr coupling to the two photons in the — ^ < 
X < ^ region. The corresponding diagram 
for the d-quark must be taken into account. 



We automatically find that the resonance exchange 
obeys 



(3.12) 



what corresponds to the symmetry imposed by Time- 
reversal and hermiticity. 
As it will be shown, the c-exchange term contributes to both the sum rule and the polyno- 
miality property. Such a resonance exchange should always be taken into account. It reflects the 
interpretation of GPDs as depending on the interval in x. In the ERBL region, the GPD can no 
longer be considered as a pure quark or an antiquark distribution but rather as a combination of 
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Figure 3.4: Total GPD for the u-quark in tlie NJL model i 
t = -1 GeV^ . 



m„=1 40 MeV and t=-1 .0 GeV^ 

1.0 I , , 




0.0 I ' ' ' 1 

-1.0 -0.5 0.0 0.5 1.0 



X 

(b) 

the case where, respectively, (a) t = —0.5 GeV'^ and (b) 



a q or q distribution as well as a meson Distribution Amplitude; what arises naturally from the 



diagrammatic description of the skewed parton distributions, i.e. Fig. 1.5 
The total distributions for the u- and d-quarks are then 

Hf'°\x,i,t) = Hf{x,C,t) + Hf{x,tt) , 

Hf'°\x,^,t) = Hf{x,^,t)-Hf{-x,^,t) = -HZ^'°\-x,^,t) , (3.13) 

with the -u-distribution being defined with the support x G [— ^, 1] and the d-distribution x £ [—1, S], 
as required. The support property has been brought up in Section [274} So was the discussion about 
the support violation; therefore we do not need to develop this argument here. 



Isospin Decomposition 



As it was the case for the PDFs, the total 7r+ distribution can be decomposed in an isoscalar and 
an isovector parts. The two terms are correspondingly defined as, in terms of operators, 



6abH'=^{x,^,t) 



dz 

dz^ 
in 



7 



(i) 



z^=z+=0 ' 



\z-^=z+=0 ' 



(3.14) 



with the pion field |7r+ ") = {\tt^) ± i\TT^)) and |7r°) = \tt^). 



Using the isospin relation (3.7) and Eq. (3.13) it is easily seen that the cr-meson exchange does 
only contribute to the isoscalar pion GPD. Effectively, the isoscalar decomposition is antisymmetric 
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m X 

H'=\x,i,t) = Hf{x,C,t)-Hf{-x,^,t) + 2HZ,{x,C,t) = -H'=\-x,^,t) . 

(3.15) 

Also, the isovector decomposition is hence symmetric in x and corresponds to the valence quark 
distribution 

r/=l('^ /- +\ tri"'''/'^ c +^ I zj'"''^ f „ & 4-\ ttI=1 



H'='{x,C,t) = HZ ix,^,t)+H: {-x,^,t) = H'='{-x,^,t) 



(3.16) 



We obtain the corresponding relation for the pions 

Hfix,tt) = \H'=\x,i,t) , 

Hf{x,i,t) = ^{H'=''{x,tt)+H'=\x,C,t)) , 

Hfix^tt) = \{H'=\x,^,t)-H'=\x,^,t)) . (3.17) 



Let us now comment on the numerical results for the pionic GPDs in the NJL model. The GPD 



for the u-quark in the NJL model is discontinuous at x = ^, see Fig. 3.4, So that the isoscalar 
GPD is discontinuous at x = and ^. The discontinuity comes from the formulation of the 
(7-exchange in the NJL model, as a realization of the chiral symmetry. This could be understood 
by considering the cr-exchange as a separate DA. In effect, the Distribution Amplitude in the 
NJL model calculation does not vanish at the end-points unlike its asymptotic form, leading to 
discontinuities at its end-points. This result has also been obtained in the Chiral Soliton Quark 
Model when the dependence of the dynamical quark mass is not taken into account jl71j , a situation 
that is similar to the one presented here. However, this peculiarity of the model calculation will 
disappear when switching on QCD just like for the Distribution Amplitudes as it has been discussed 
in Section [221 



Sum Rule and Polynomiality 

We can now check that the properties of GPDs are respected in this model calculation. The sum 
rule for the isovector GPD reads El 

j dx^QgHf{x,tt) = j dxH^=\x,^,t) = 2F^{t) , (3.18) 



with -F7r(t) the electromagnetic Form Factor given by Eq. (2.37). This sum rule is respected and 
has been analytically verified. Regarding the isoscalar combination, the sum rule involves the 
gravitational Form Factors of the pion 9\^2{i) in the following way, 

j dxxH^=°{x,^,t) = ei{t)-f92{t) . (3.19) 



^The association J^^ dx "^^g QqHq^ {x,^,t) = J^^^ dx H' ^{x,^,t) is valid under the integral over dx only. 



42 



Pion Generalized Parton Distributions 



1.0 



11 



0.4 
0.2 



m„=14D MeV and ^=0 









^ t=DGeV^ - - 




t=-0.5 GeV^ 




t=-1.DGeV^ 



0.0 



0.2 



0.4 



0.6 



1.0 



Figure 3.5: GPDs in the NJL model in the case m.^ = 140 MeV and C = . 

Since this sum rule refers to the momentum carried by each quark when taken at t = GeV^, we 
expect, in the chiral limit, ^i(O) = ^2(0) = 1- This condition has been fulfilled as it can be seen in 
Table 13JJU 



The polynomiality property is also respected, 



2i 



(3.20) 



i=0 



where at the second line we have explicitly expressed the i^Ar(^, t) as polynomials in ^ of order not 
higher than n. On virtue of the relation 



Hgix,-C,t) 



given by Time-reversal invariance, only even powers in the skewness variable are allowed in the 
previous expansion in polynomials. The numerical results for the coefficients in both cases uIt^ = 



MeV and = 140 MeV are given in Table 3.1 



The chiral limit, where itIty = 0, and with zero momentum transfer, is perfectly defined in 
a regularization-independent way. It makes sense to study separately the ^ dependence of the 



1 



2? 

o{i-x)e{i + x 



- x)e{^ + x) + 9{x - 0^(1 - x) 



contributions coming from the direct and crossed diagrams of Fig. 3.1 and the a-exchange of Fig. 3.2 

1 

2 

1 X 

2 2^ 

The total GPD for the n-quark in a vr^ reads 

1 



(3.21) 



d{i - x)e{i + x) + e{x - i)e{i - x) 



(3.22) 



^Notice that Table |3.l| is given 7?^^'°' which is defined between [— ^, 1] so that the results have to be mutliplied 
by 2 when considering H and we extend the integral to [—1, 1]. 
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n 


1 


2 


3 


4 




^1,0 (i) 




A3,o(t) Asait) 


^4,0 W Ai,2{t) AiA{t) 


rriTr = 




t = 

t = -1 
t = -10 


1. 

0.487 
0.091 


0.5 -0.5 
0.336 -0.240 
0.119 -0.093 


0.333 0. 
0.257 -0.049 
0.118 -0.058 


0.25 0. -0.25 
0.208 -0.030 -0.113 
0.111 -0.046 -0.044 


= 140 




t = 0[-10-'^] 
t = -1 
t = -10 


1. 

0.482 
0.090 


0.5 [-0.473] 
0.332 -0.230 
0.116 -0.088 


0.332 [0.005] 
0.253 -0.045 
0.114 -0.055 


0.247 [0.005] [-0.236] 
0.203 -0.026 -0.109 
0.108 -0.043 -0.042 



Table 3.1: CoefBcients of the polynomial expansion. The pion mass is expressed in MeV and t is expressed in GeV^. 
Values between brackets correspond to t = — 10~^ GeV^ [153| . 



Evaluated for the skewed PD, the expression of the u-resonance is model dependent. However, 
since the interpretation of the cr-exchange is unique in the limit m^r — )• 0, chiral symmetry im- 



plies that the diagram 3.2 cancels out the "cr-exchange like" contribution contained in H^. In 
other words, we recover the result i/^^°(a;, 1, 0) = predicted by the chiral symmetry jl71j . 
Also, the isoscalar and isovector combinations read 



-ei- 



00(1 + x) + e{x-oo{i 



= 9{-x - 09{1 + x) 

+ 0{^-x)9{^ + x) + e{x 



what is consistent with previous results |171| . This 



result is illustrated on Fig. 3.6 



We notice also that, in the forward limit, i.e. for 
^ = t = 0, we find the quark distribution function 



2.0 r- 

1.5 - 

1.0 

0.5 - 

0.0 - 

-0.5 - 

■1.0 — 

-1.5 - 

-2.0 — 



1=0 

1=1 



Figure 3.6: Chiral limit for ^ = 0.4. 



H^'°'ix,0,0) = qix) = l , 
where, as expected, the cr-exchange does not contribute. 



To this point, we have presented the calculation of the chiral-even twist-2 GPD of the pion in the 
NJL model using a fully covariant Bethe-Salpeter approach. The only approximation employed has 
been the ladder approximation. We have solved the BS equation exactly thanks to the point-like 
interaction of the NJL model. The usual properties of GPDs have been recovered. Nevertheless 
discontinuities are encountered, at the boundaries and are associated to the off-shellness of the con- 
stituent quark. Also the expression of the cr-exchange in the NJL model introduces discontinuities 
at the X values — ^ and ^. 



44 



Pion Generalized Parton Distributions 



3.2.1 Double Distributions and the D-term 

An alternative way to describe the hadronic matrix elements defining the GPDs is through the 
so-called Double Distribution (DD) parameterization. By Lorentz invariance the matrix element 



( 3.2 ) can only depend on the two independent variables p ■ z and A-z. Hence the pion DD \n3\ I174j 
is defined as a 2-dimensional Fourier transform in these two independent variables. The analysis 



of diagrams of the type 3.1 & 3.2 yields the behavior of these variables; 



MUp-z,A-z;t) 



for the whole matrix element. 



vr+(P)) 



z^ 

2J ' " V2. 

1 /.1-I/3I 

-1 J-{\~\fi\) 



2{p- z 

-(A • z) da e*"(^-^)/2 t) 



(3.23) 



/3p- (l+a)A/2 




p- A/2 



/3p+(l-a)A/2 



p + A/2 



Figure 3.7: Momenta associated with the quarks and hadrons in the 
Double Distribution parameterization. 



Since the variables p ■ z and A • z 
are treated as independent, the DD 
refers to fixed momentum transfer t 
but with unfixed individual hadron mo- 
menta P and P' . In other words, 
these momenta do not depend on the 
skewness variable. The interpretation 
of the latter being process dependent, 
one can say that the matrix element 
on the r.h.s. of Eq. (3.23) contains 



a process-independent information and, 
hence, has a quite general nature. 
The dependence on /3 is of the parton distribution type, while the one on the parameter a 



resembles the dependence of a distribution amplitude. This feature is illustrated in Fig. 3.7 



From the naive analysis of GPDs as pa- 
rameterized through the DD, one would 
write the whole contribution coming from 
diagrams such as 3.1 or 3.7 in the form 
fg{f3,a,t); i.e. (3.23). From our experi- 



Rn 




ence in model calculations, we know that 
diagrams of this type may contain a reso- 
nance exchange in the t-channel whose ten- 
sorial structure corresponds exactly to the 
one of the remaining non-resonant part of 
the diagram. The behavior of this meson- 
exchange contribution with respect to the two independent variables p ■ z and A • z is manifestly 
different. 

As a matter of fact, it has been shown in Ref. [171] that another twist-2 "DD structure" has to 
be considered indeed. The origin of this new structure is the subtraction from the matrix element 



Figure 3.8: The D-term as a a-resonance in the t-channel. 
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of those resonances in the t-channel (Fig. 3.8 ). Namely they are terms depending solely on the A • z 



variable, i.e. non-vanishing terms in the limit 

p - z . 



(3.24) 



The function P(a; t) is therefore defined taking the latter limit and subtracting this result from 
the p • z-dependent contribution, i.e. 

MI{t>- z,^- z;t) = A^^(p-z/0,A-z;t)-A^^(0,A-z;t)+A^^(0,A-z;t) . 



The meson-exchanges must have the characteristics of Distribution Amplitudes so that we expect 
them to be parametrized by a. In terms of DD, we call the corresponding function the L'-term jl71j . 
The Z?-term is uniquely defined for a given hadronic matrix element. 



When applying the limit (3.24) to the matrix element as expressed in the NJL model calcu- 



lation (3.6) taking into account all diagrams, this r eson ance exchange corresponds in part to the 



term going like /2.A of t he integral decomposition (3.9) as well as to the term coming from the 
second diagram Fig. 



3.2 



Th e wh ole resonance exchange also includes a x-odd part of the /2,p as 



well as of the Is integrals in (3.9). 



The Double Distributions can be reduced to the usual skewed parton distribution by taking a 
one-dimensional cut, i.e. imposing a relation between p ■ z and A • z. This reduction formula reads 



1 r-i-m 

dp 

-1 J-ii-m 



da 6{x - P - Ca)T^il3,a,t) + sgni^D^ [-,t 



(3.25) 



The Z?-term obeys the symmetries of an isoscalar resonance (3.12). Therefore it contributes solely 
to the isoscalar GPD as expected. 

The parameterization through Double Distribution applies at the level of the Fourier transform 
of the matrix element. This renders the evaluation of the DD model independent. Nonetheless 
this statement is quite misleading in the sense that the relevant diagrams depend on the model 
considered. So, from a technical point of view, one has to evaluate the diagrams that apply, and, 
with the Feynman rules resulting from the features of the model. The GPDs appear to be a 



combination of two- and three-point functions, e.g. in the NJL model Eqs. (3.9, 3.11), which, in 



turn, can be expressed in terms of Double Distributions through the reduction formula. An explicit 
example would be the expression of the 3-propagator integral 



/•I /-i-l/SI 

h{x,^,t) = dp da6ix- P-aOM/3,a,t) ; 

7_i J-{i-m) 



(3.26) 



with X3 a Double Distribution given in (F.6). The 2-propagator integrals that do not behave like 
DAs can also be included in the form I{P,a,t). This reduction is totally model independent. 
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The model's description of the a-exchange, illustrated in Fig. 3.2, also accounts for a formal 
and general description on the form of the r.h.s. of Eq. (3.23). In particular, it is of a pure D-term 
structure. The total description of the matrix element simply results from the sum of the two 
contributions; namely as given by the reduction formula (3.25) plus the corresponding formula 
for H"" 



1.5 



1.0 



0.5 



0.0 



-0.5 



-1.0 



-1.5 















- 


- 

D" - - . 




F" ■ 




H" 



-0.4 



-0.2 



0.0 



0.2 



0.4 



0.6 



0.8 



1.0 



As already mentioned in Section 2.1 the 
difference between the quark models is man- 
ifest in the regularization scheme. The cho- 
sen regularization scheme is an integrant 
part of the model since it depends on the 
model's interaction and conveniences for the 
adaptation to the problem under scrutiny. 
Concretely the Double Distributions are re- 
duced to the physical propagator integrals 
by imposing a regularization scheme. In 
the NJL model as presented in the previous 
Chapter [2| it would mean that the prop- 
agator integrals follow, in the case of 3- 



Figure 3.9: The contributions from the double distribution F'^ 
(blue dotted line), from the D-term (dashed line) and the total 
distribution (red curve). 



propagator integrals, the expression (3.26) 
to which we impose the Pauli-Villars regu- 
larization scheme, i.e. Eq. (C.15). 
In Ref. 



5] the Double Distributions of 
the pion in the NJL model are given using 
a Pauli-Villars regularization in the twice- 
subtracted version of Ref. [180], as well as in the Spectral Quark Model jl81| . Both calculations 
yield results comparable to the calculation of Ref. |153] . 



One of the main advantages of the DD parameterization is that the expression (3.25) manifestly 
leads to the correct polynomiality properties of the GPDs. Effectively, when taking the Mellin 
moments of the reduction formula, one immediately finds 



n-l 



1-1/31 
'(1-1/31) 



da5{x- /3-^a)J^{/3,a,t)+sgn{ODg [-,t 



(3.27) 



where the first term on the r.h.s. is a polynomial at a power of order at most n — 1 of the skewness 
variable ^. The D-term instead gives the highest power in ^, i.e. n. |^ 



Concerning the chiral limit, the DD approach yields the same results (3.21) as the skewed PD 
evaluation |171| . In this limit, the D-term is therefore reduced to the zeroth order term of the 
resonance exchange coming from the /2,p term. 

^The latter statement can be visualized by performing the change of variable x/(, = v; which carries a ^" outside 
the integral over, then, f. 
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Let us illustrate the previous statements. In the chiral limit, the D-term reads, for the n-quark, 



^e{x + oe{-x)-^9ix-oo{x) 



(3.28) 



The remaining part is found through the reduction formula, providing that 



Using (3.22), we obtain 



F:ix,^,t) + DZ(-,t 



F:{x,^,0) = 9{x)9{l-x) 



(3.29) 



which can obviously not satisfy the polynomiality property. By inserting the D-term (3.28), we 
obtain the expression for the total GPD (3.29). In particular, we recover the expression (3.22) 



which, in turn, fulfills the required conditions. These results are illustrated in Fig. 3.9 



There is considerable freedom in the tensorial decomposition of the matrix element Aiq- An 
alternative decomposition to the expression (3.23) is given in Ref. |198j . 



Mq {p- z,A- z;t) 



2{p-z 
-{A-z) 



dp 



-1 



-(1-1/31) 

1-1/31 



vr+(P)) 

da e-i/3(p-^)+^"(A-^)/2 J^n^p^ ^) 



dp 



(1-1/31) 



(3.30) 



with the function ICg{P,a,t) even in a and the function Q1^{P,a,t) a-odd, as a consequence of 
Time-reversal invariance. The corresponding reduction formula reads 



H^ix,^,t) 



dp 



(1-1/31) 



da 5{x-P- ia) [mp, a,t)+i g^{P, a, t)] 



(3.31) 



The sum rule for IC^ gives the electromagnetic pion Form Factor, as expected; whereas the sum 
rule for Q'^ gives since it is an odd function of a. 

One could relate the above given expression to the usual DD decomposition by writing all the 
P dependence of Gq{P, a, t) as a JCq{P, a, t), leaving a pure A • z in the second term. The decom- 
position (3.30) is however calculationally the most useful, since it does not evoke the cumbersome 
limit (3.24). On the other hand the expression of the meson-exchanges is less explicit. 

The distribution amplitude character of the D-term implies it obeys to the ERBL evolution 
equations, that has been introduced in the previous Chapter. The total GPD does however not 
have the same behavior when going the higher Q^. From Fig. 1.5, we have defined the three regions 
as DGLAP-regions for the support in x being either [C, 1] or [—1,—^] and ERBL-region for the 
X G [— ^, ^] region. Those regions were named after the evolution equations they obey. 
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3.2.2 QCD Evolution 

The form of the Leading-Order (LO) QCD evolution has been given, for the standard PDF evolu- 
tion, by Dokshitzer ^85j , Gribov and Lipatov |104l I141j as well as by Altarelli and Parisi [9j . It gave 
rise to the DGLAP equations, whose evolution will be followed by the two regions of emission and 
reabsorption of a quark or antiquark, respectively. On the other hand, the equations describing the 
evolution of the distribution amplitudes were given by Efremov and Radyushkin f88] as well as by 
Brodsky and Lepage [140j . 



The form of the integral-equation defining the evolution equations depend upon the nature of 



the distributions. In Eqs. (3.15 3.16), the singlet or sea distribution have been defined as well as 



the non-singlet or valence distribution 



(3.32) 



(3.33) 



Still, one main ingredient of QCD evolution has to appear: the gluon distribution, obeying the 
following symmetry 



xHg{x,i,t) 



-xHg{-x,i,t) 



with the gluon distribution being zero at the scale of the model, i.e. Hg{x,^,t,QQ) = 0. 

It is important to notice that the non-singlet distribution does not mix with the singlet under LO 

evolution. However, the singlet and gluon distribution do mix under evolution. 

The explicit form of the LO QCD evolution equations for the GPDs can be found in Refs. [1031 
11181 11231 11241 11481 1174] . The effect of LO evolution on the non-singlet, singlet and gluon distribu- 
tions is depicted on Fig. 3.10 The evolution code of Freund and McDermott 



has been used 

with the value for Qq and A as given in (2.30) and in the paragraph preceding the latter expression. 



On both panels the LO evolution is performed for the initial scale given by NJL model as given by 



the expression (2.30). The results presented here are in agreement with the evolution for the NJL 



model as introduced in Ref. 

The initial distributions are represented by the thin pink curves. We can see that the isovector 
GPD (on the top) differs in the ERBL-region from the isoscalar distribution (middle of the panel), 
where the D-term, or u-resonance, contributes due to its antisymmetry. 

On the left panel, the thick green lines represent the distributions evolved at Q = 4 GeV. One 
is immediately struck by the fact that evolution changes completely the end-point behavior as well 
as smoothes the discontinuities between the two regions, i.e. at x = it^. Even if the evolution 
seems strong here, it is mentioned in Ref. [45j that the desired continuity of the functions at the 
end-points is achieved already at values of Q infinitesimally larger than Qq. Once more, we refer to 
the case of the distribution amplitude where a numerical analysis of the truncation of the solution 



of the evolution equation has been driven, see Fig. 2.6 



An important feature of the evolution in of the GPDs is that the evolution equations in the 
ERBL-region depend on the value of the GPDs in both the ERBL-region and the DGLAP-region, 
while the DGLAP evolution equation do only depend on the value of the GPDs in its very own 
region. In other words, the GPDs in the DGLAP-regions feed the ERBL-region. Moreover, as for 
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Figure 3.10: Results of the LO QCD evolution from the NJL model for several values of t — —1 GeV^ and for 
5 — 1/3. From top to bottom, respectively, the non-singlet, the singlet and the gluon distributions. The sohd pink 
lines represent the initial conditions. Left panel: The thick green lines are the evolved result to Q = 4 GeV and the 
dashed lines the asymptotic forms. Right panel: The grey lines are the evolved result to different values of Q (1,8; 
4; 10; 100 and 1000 GeV) and the dotted red line the asymptotic form. 
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the DA, evolution pushes the distributions towards low x's. This can be observed in both the left 



and right panels of Fig. 3.10 



As a matter of fact the asymptotic forms for the evolved distributions are solely defined in 
the ERBL-region. The (anti)symmetry of the distributions around the point x = shown in 



Eqs. (3.32, 3.33) remains unchanged under evolution. The analytic asymptotic solutions were 



found in Ref. |174j . We adopt the forms given in Ref. [45j 



xH3{x,i) 



3 e 



2i e 



Nf 



X e 



x^ 1 



15 



15 4Cf 



8^ ACf + N 



f 



X 



[ei{t) - e 92{t)] 



(3.34) 



for — ,f < X < ^ as well as for fixed t- values. The i-behavior is given by the sum rules (3.18 



3.19); a first principles property that remains unchanged after evolution. The presence of the gluon 



distribution nevertheless affects the isoscalar sum rule (3.19) 
■1 



J dx {xH^{x,tt,Q^)+xHg{x,^,t,Q^)) = ei{t,Q'^)-^'^92{t,Q'^) . (3.35) 



The values of the gravitational Form Factors for different t values are given in Table 3.1 



The asymptotic behavior described by Eqs. (3.34) is illustrated on Fig. 3.10 by the dashed lines. 
On the right panel, it is clearly shown that the evolution is slow in reaching the solution for Q — t- oo 
of the GPDs. On the other hand, the evolution is faster at low values of Q. This result we could 
expect as it happens for the pion DA, which obeys the ERBL evolution equations. In Fig. 2.7 is 
depicted the slow evolution of the pion DA at high Q-values. 

Also, in order to improve in reaching the asymptotic behavior, one should start the evolution 



from a different value of Qq. However, as it has been explained in Section 2.5 the scale is fixed by 
the model's results for the second moment of the PDFs. Even going to the NLO could not change 
such a peculiarity of the model calculations. Effectively, the initial scale for a NLO evolution is 
found to be higher than for the LO one. However, this lower Qq for the LO evolution compensates 
for the effect of NLO evolution: the results are consistent one to the other. 



3.2.3 Prom the Experimental Side 

From the experimental point of view the GPDs of the pion are difficult to determine. Pions decay 
into muons or photons and therefore the pion Distribution Functions cannot be obtained from 
direct DIS experiments. Rather they have been inferred from Drell-Yan process \26\ \3T\ |65] and 
direct photon production [2Jj in pion-nucleon and pion-nucleus collisions. In a near future, the 
CLAS++ detector with the CEBAF accelerator at 12 GeV at JLab should allow the observation 
of pion GPDs through large angle Compton scattering |54] . 

It has been recently proposed, in Ref. [10], to access the pion GPDs by studying the DVCS on 
a virtual pion that is emitted by a proton, i.e. the process ep — )• e 7 vr"'" n. 



Transition Distribution Amplitudes 



Here comes the heart of the thesis. The study of Generahzed Parton Distribution is extended to 
distribution amphtudes for transitions. 



Colhsions of a real photon and a highly virtual photon are an useful tool for studying funda- 
mental aspects of QCD. Inside this class of processes, the exclusive meson pair production in 7*7 
scattering, 

ill ' ill P^^" > (4-1) 



i*{q.) 



has been analyzed in Ref. [131^ I165j . This process is particularly interesting since it implies, in 
different kinematical regions, different mechanisms. 

At small momentum transfer t and in the kine- 
matical regime where the photon is highly vir- 
tual, a separation between the perturbative and 
the nonperturbative regimes is assumed to be valid. 
Through the factorization theorems the amplitude 
for such reactions can be written as a convolution 
of a hard part M^, a meson Distribution Amplitude 
(j)M and another soft part, describing the photon- 
pion transition by means of the Transition Distribu- 
tion Amplitudes. 

On the other hand, at small s, the factorization 
corresponds to the GDA's mechanism. The inter- 
play between the two mechanisms in the kinematical 
region s,t <C may lead to an interesting duality 
property. This will be considered in the Future Per- 
spectives. Also, the kinematical regime where both 
t,s ^ corresponds to large angle scattering d la 
Brodsky-Lepage and will not be considered here. 




Figure 4.1: Factorization diagram defining the TDA 
for the process 7*7 — >■ ttM at smail momentum 



transfer, t 



■ Pm) , and large invariant mass, 
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The accurate definition of these objects was first presented for the mesonic TDAs, by simphcity. 
In particular, it has been first analyzed for the pion-photon transitions |165j involved in related 
processes, like hadron annihilation into two photons vr+vr" — 7*7 and backward virtual Compton 
scattering 7*7r"^ — )• 77r^. 



The concept of transition distribution was first introduced by Frankfurt, Polyakov, Strikman 
and Pobylitsa [91^ [92] as a study of the transitions 7* + p — )■ Baryon + Meson through Skewed 
Distribution Amplitudes. 

In this Chapter, we compile the results of 
our publications on mesonic Transition Distribution 
Amplitudes calculated in the N JL model [671 EH] • 

We establish the connection between the TDA 
and the vector and axial pion Form Factors, Fy and 
Fa- The definition of the vr — )• 7 TDAs are ampli- 
fied to 7 — 7- TT transitions. The calculation and the 
results in the NJL model are given in details for the 
skewed distribution approach. The way to Double 
Distributions is given that enables us to compare 
our results with both phenomenological parameter- 
izations and other model calculations. 
In the next Chapter, as an application of the numerical results, we will use the defined quan- 



Figure 4.2: The 7r+ 
plitude. 




TDA 



7 Transition Distribution Am- 



tities to estimate the cross sections for exclusive meson production in 7*7 scattering (4.1) in the 
same kinematical regime. 



4.1 The Pion Transition Form Factors 

The pion-photon TDAs are connected, through sum rules, to the vector and axial-vector pion Form 
Factors, Fy and F^- Before giving a proper definition of the TDAs let us recall the definition of 
these Form Factors. They appear in the vector and axial vector hadronic currents contributing to 



the decay amplitude of the process vr"'" — >■ 'ye'^u depicted in Figs. |4.3| fc |4.4[ The precise definition 
of these currents is |147[ | 



(7(^7)1 ^(0)7^r-g(0) {tt (p^)) = -iee'' e^.p^pPp^ , (4.2) 



(7(^7)19(0)7^751" g(0) \Tr{Pn)) = ee" {p^^p^^ - g^^py -p^] 



FA{t) 



+ e ( (p-Y - Ptt)^ P-ku 2 - a^iu I , (4.3) 



with = 93MeV, e^^"^^ = 1 and = (ti — iT2) /2 as given in Appendix [a} All the structure of 
the decaying pion is included in the Form Factors Fy and Fa- We observe that the vector current 
only contains a Lorentz structure associated with the Fy Form Factor. On the other hand, the 
axial current is composed of two terms. The first one, defining Fa, gives the structure of the pion 
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Figure 4.3: Diagrams contributing to the process n — >■ 'ye v: the structure-dependent contributions governed, 
respectively, on the l.h.s. by the vector and, on the r.h.s., the axial Form Factor. 




Figure 4.4: Diagrams contributing to the process tt^ — >■ 76^1^: the electron (l.h.s.) and the pion (r.h.s.) inner 
Bremsstrahlung. The electron Bremsstrahlung does not contribute to the tt — )■ 7 hadronic current. 



and is gauge invariant. The second one corresponds to the axial current for a point-like pion and 
is not gauge invariant. This term coming from the pion inner Bremsstrahlung has, in turn, two 
different contributions. The first one corresponds to a point-like coupling between the incoming 
pion, the outgoing photon and a virtual pion which is coupled to the axial current. It is depicted 
in the diagram of Fig. |4.5| and can be seen as a result of PCAC, because the axial current must 
be coupled to the pion. It isolates the pion pole contribution from the axial current in a model 
independent way. The second contribution of this term, proportional to fn9fiu, corresponds to a 
pion-photon-axial current contact term. With these definitions, all the structure of the pion remains 
in the Form Factor Fa_. The total gauge invariance of the amplitude for the pion radiative decay is 
ensured by the electron Bremsstrahlung contribution, whose amplitude does not contribute to the 
hadronic currents (4.2|4.3). 



4.2 From Hadronic Currents to Transition Distribution Ampli- 
tudes 

Let us go now to TDAs. 

For their definition we introduce the light-cone coordinates as defined in Appendix [A) 
The skewness variable describes the loss of plus momentum of the incident pion, i.e. 

{ = . (4.4) 

with p+ = {pt, + p^)^ /2. 
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The explicit expressions for the pion and photon momenta in terms of the light-front vectors are 

1 



Here we have 



(1 - e) + 2 



A^ = {p,-p^Y 



(4.5) 



0, A-L,0 



Given the kinematics here above, it is easy to deter- 
mine 



A 



±2 



2e(l-e)m2 - {l-e)t 



(4.7) 



which gives the expression for t 

2e(l-e)m2-A^2 



(4.6) 



-i^\/\/\/\y 



1_^2 

The value of the skewness variable is constrained to 
range between 



(4.8) Figure 4.5: Pion pole contribution between the ax- 
ial current (represented by a cross) and the photon- 
external pion vertex associated to the last contribu- 



tion of Eq. (4.31 



{2ml -t] 



< e < 1 • 



(4.9) 



Actually there is no symmetry relating the distributions for negative and positive which could 
have constrained the values of the skewness variable to be positive, like for GPDs. 



The polarization vector of the real photon, e, must satisfy the transverse condition, 



e.p^ = , 

as well as an additional gauge fixing condition. When going from the hadronic currents to the 
parton distribution amplitudes, we need e.n = /p~^ to kinematically become higher-twist, i.e. 

£.n — )• when p"*" — )• oo 

The standard gauge fixing conditions, namely e'^ = or = 0, both satisfy the previous require- 
ment. In fact, all what we need is that the components of the polarization vector remain finite 
when p'^ goes to infinity. 

A general expression for e is 

e = e^ + a^ + pp-^n . (4.10) 
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The transverse condition for the real photon polarization provides us for a first constraint, e.g. on 



1 



1 



t 



Therefore the scalar product e • A is 

£• A 



+ o 



p 



(4.11) 



where the gauge dependence inclosed in a is kinematically higher-twist, i.e. goes like l/p~^- Also, 
the decomposition e ■ A depends on the kinematics, what is explicit through the 1 — ^ in the latter 
expression. Henceforth the tensor decomposition would rather be chosen to be e-^.A-^. 



The vector and axial TDAs are the Fourier transform of the matrix element of the bilocal 
currents separated by a light-like distance 

q{-z/2)j^[j5]qiz/2) . 

Then, they are directly related to the currents defined in Eqs. ( 4.2||4.3 ) through the sum rules 

dz~ 



dx 



-1 



27r 



(7 (1^7)1 9 -o [^5] T ^ - |vr (p^)) 



2 + =2-L=0 



1 

P+ 



(7(^7)1 ^(0)7|. [75] T-q{0)\TT{p^)) 



(4.12) 



With this connection we can introduce the leading-twist decomposition of the bilocal currents. For 
that we need the light-front vectors and defined in Eq. (A.8), which leads to the expres- 
sions ( [45p^ . 

To leading-twist, the 7r+ — )• 7 TDAs are defined 



dz 



2^ e-^^^- (7 {p,)\ q {-"-) f>T~q (0 |vr+ (p.)> 



5+=Z-L=0 



(4.13) 



/ 



^e-^"- (7 {p,)\ q {-'-) frt,r-q (0 |vr+ (p.)) 



z+=z^=0 



—e £ ■ A ^ -e (e • A) — ^ e (t) ml — — 



(4.14) 



where e (^) is equal to 1 for ^ > 0, and equal to —1 for ^ < 0. Here V{x,£,,t) and A{x,S^,t) are 
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respectively the vector and axial TDAs. They are defined as dimensionless quantities. From the 



condition (4.12) we observe that they obey the following sum rules, 



ra-n 



C dx (x, e, t) = ^^^Fa (t) 
J-1 rriT, 



(4.15) 
(4.16) 



In the second term of Eq. (4.14), we have introduced the pion Distribution Amplitude (DA) 
x), whose expression in given by Eq. (2.1). The pion DA </>(x) vanishes outside the region 



X G [0, 1] and satisfies the normalization condition (2.2). 



This second term has been introduced in order to iso- 
late the pion pole contribution of the axial current in a 
model independent way, as we have done in Eq. (4.3) 
for the vr"'" — )• ^e^u process. Therefore, all the struc- 
ture of the pion remains in the TDA A [x, ^, t). It can 
be seen as a result of PC AC, because the axial current 
must be coupled to the pion. Therefore, this term is 
not a peculiarity of the pion-photon TDAs. A similar 
pion term will be present in the Lorentz decomposition 
in terms of distribution amplitudes of the axial current 
for any pair of external particles. We emphasize that it 
is a model independent definition, because we have de- 
fined the numerator of the pion pole term as the residue 
at the pole t = m^. With this definition, all the struc- 
ture dependence related to the outgoing vr^ is included 
in A (x, ^, t). Moreover, the pion pole contribution can 
be estimated in a phenomenological way, as we will see later on. 



(i + ?)p+ 




Figure 4.6: Pion pole contribution to the axial 
bilocal current corresponding to the last term of 



equation (4.141 



A pion exchange contribution was analyzed in |145| 1160] for the axial helicity-fiip GPD and, in 
|199j . a similar structure for the axial current was obtained using different arguments. This term 



we have represented in Fig. 4.6 is only non-vanishing in the ERBL region, i.e. the x S [— i^, ^\ 
region. The kinematics of this region allow the emission or absorption of a pion from the initial 



state, which is described through the pion DA. And it can be seen from Fig. 4.6 that positive values 
of ^ corresponds to an outcoming virtual pion, whereas negative values of ^ describe an incoming 
virtual pion. The latter is related to the matrix element {-k {pt^)\ g^75r~g|0) , instead of the one 



present in Eq. (2.1 ), what gives rise to the minus sign included in e (^) . 



In order to make possible the connection with other works, it must be realized that the axial 
TDA must be defined from an electromagnetic gauge independent tensorial structure. As above- 
mentioned, the product e ■ A contains gauge dependent terms. Since the TDAs are gauge indepen- 
dent quantities, the tensorial structure defining them has to be itself gauge independent. The only 



possibility is hence, see Eq. (4.14) 



(4.17) 
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which is equal to, e.g., (1 — ^) (e.A) in the = gauge we have chosen to work in. 



In many papers present in the literature, the r.h.s. of Eq. (4.14) contains only the A term. In 



a general case the A term and the pion pole terms have different tensor structures as shown by 



Eq. (4.3). However we can fix the gauge convention and the frame conventions in such a way that 



these two tensorial structures coincide. In other words, should we have chosen another gauge and 
a frame in which all the transverse component is in the pion momentum, i.e. = 0, the tensorial 
structure of the axial TDA and the pion pole would have coincided: we have the freedom to choose 
En with only transverse component. In that case, 



e- A 



Only with this particular definition, we can change our definition of the axial TDA to 

dz~ 



27r 



(7(^7*^)1 



7 '75 T 



1 



7-/1 - V 2J ' ' " V2 



z+=z^=0 



with 



ml — 



2i 



(4.18) 



(4.19) 



The latter expression shows that the pion pole contribution to the axial TDA is closely related to 
the D-term of the Generalized Parton Distributions |171) . However it must be realized that in this 
case it gives an explicit contribution to the sum rule 



dx 



(4.20) 



In what follows, the definitions derived from the tensorial decomposition (4.13, 4.14) of the hadronic 
matrix element for the pion radiative decay will mainly be used. 



4.3 A Field Theoretical Approach to the Pion-Photon TDA 



In Chapter 2.1 we have defined a method of calcu- 
lation for the pion DA in a field theoretical scheme, 
treating the pion as a bound state of quarks and an- 
tiquarks in a fully covariant manner using the Bethe- 
Salpeter equation. We then applied the same method 
for evaluating the pion GPDs in Chapter [3j The 
same is done here for the evaluation of the pion- 
photon TDAs |67j . This method has enormous ad- 
vantages because it preserves all the physical invari- 
ances of the problem. Therefore, any property as 
sum rules or polynomiality is expected to be pre- 
served. 



1 -1-^ 

Figure 4.7: 
TDAs. 



Diagram representing the tt 
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(a) (b) 



Figure 4.8: Diagrams contributing to the TDA. We have depicted diagrams in which a quark u is changed into a 
quark d by the bilocal current. There are similar diagrams in which the antiquark d is changed into a u 



As usual, we consider that the process is dominated by the handbag diagram. Each TDA has 
two related contributions, depending on which quark [u or d) of the pion is scattered off by the 



deeply virtual photon. In Fig. 4.8 we have depicted the diagrams in which the photon scatters off 
the M-quark. We observe that there are two kinds of contributing diagrams. In the first one the 
d-antiquark appears as the intermediate state, while in the second the bi-local current couples to a 
quark-antiquark pair coupled in the pion channel. The latter is present only for the axial current 
and includes the pion pole contribution. 

The details of the method of calculation are given in Ref. [153j . In the present case we obtain. 



from the first kind of diagram of Fig. 4.8, the following contributions |67j 

1 
' 2 



27r 
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+ r3 -f" iS (p-y - k) -f^b5]T'iS{p^-k) ik,p^) iS{-k) 



5 ( xp'^ 



Tr 



1 A 

2 3+"^ 



'2k^ 



pZ -p^ 



Y iS (k) cD-^ (A;, p^) iS {k-p^) [75] T- iS {k - p^) 



(4.21) 



where S {k) is the Feynman propagator of the quark and {k,pT^) is the Bethe-Salpeter amplitude 
for the pion. Here Tr() represents the trace over spinor, color, flavor and Dirac indices. 



The first contribution in Eq. (4.21) is the one depicted in the first diagram of Fig. 4.8 The 



second contribution corresponds to a similar diagram but changing quarks u and d. In the NJL 
model, ^'^^ {k,p-n) is as simple as. 



(4.22) 



where g-i^qq is the pion-quark coupling constant defined in Eq.( C.30) 
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4.3.1 Vector TDA 



The vector TDA has contribution only from this first kind of diagrams. We can express V {x, ^, t) 
as the sum of the active u-quark and the active d-quark distributions. The first contribution will 
be proportional to the d's charge, and the second contribution to the u's charge. Therefore, we can 
write 

V^Ux,^,t) = QdvZtdi^,^,t) + Quv}!^-{x,^,t) . (4.23) 
Isospin relates these two contributions, 

vll^dix,tt) = v{x,(„t) , 
v^l,^ix,tt) = v{-x,i,t) . 



A direct calculation gives 



with 



TDA 



'x,^,t) 



6{x-l + 



dnqq ^ -'Si) 



k+\ 1 

/ — + ie) {pj — k)^ — + ie [p-jr — k) 



(4.24) 



2 2 I • 



(4.25) 



given by the expression (E.15) in Appendix. 



In this integral, we first perform the integration over k~ . The pole structure of the integrand 
fixes two non- vanishing contributions to ?;(x,^,t), the first one in the region ^ < x < 1, corre- 
sponding to the quark contribution, and the second in the region — < x < ^, corresponding to 
a quark-antiquark contribution. Given the relation (4.23), the supp ort o f the entire vector TDA 



{x,^,t), is therefore x G [—1,1]. The analytical expression for (4.25) is given by Eqs. (E.15' 



E.17). 



For the vr*^ 
to V (x, ^, t) 



7*7 process, the contributions of the first type of diagram (Fig. 4.8) can be related 

-^^ (vix,C,t) + v{-x,^,t)) , 

{v{x,C,t) + v{-x,^,t)) . (4.26) 



V2 

Qd 



4.3.2 Axial TDA 

Turning our attention to the axial TDA, we find a new contribution arising from the second diagram 
of Fig. 4.8 This second contribution comes from the re-scattering of a qq pair in the pion channel. 
Therefore 



27r 

(|i:2TD + Vm'- 



7(^7'^) 



(Ptt) 



2ig 



2gllps (t) 



(4.27) 



60 



Transition Distribution Amplitudes 



where i is an isospin index, 



{ -Tr 



c/)""^ (/c'jPtt) iS {-k') ^ ( ^ + '^3 ) -i" iS{p^- k') i-fsr' iS (p^ - k') 



-Tr 



{k',P7T) iS [k' - p^) i'j^r' iS {k' - p^) 2 (^3 + ^3 J iS [k') 
I (2^"^ (^^^ -l{Py+Pt- 2A:+)) Tr [iS {p^ - k) ^^,7- iS {p^ - k) i^^r^] 



and Yips is the pseudoscalar polarization, see Fig C.3 described by Eq. (C.23), 



Hp, (A2) = -i [ ^ {Tr [i75 iS {k) i-i^ iS (A - k)] ; 
J [2tt) 



(4.28) 

i 

(4.29) 
(4.30) 



We can now evaluate the axial current in a straightforward way. Nevertheless, in order to 
extract the axial TDA we must subtract the pion pole contribution. We need for that the pion 



amplitude studied in Chapter 2.1 In the NJL model, the pion DA is given by Eq. (2.13) in the 
[-iii] region. 



1 

2^ 



2i 



V2U 



(4.32) 



The trace over the Dirac space involves a /c-dependent 3-propagator integral that renders a 
bit more tricky the calculation. After a long but direct calculation, we obtain the expression for 
A{x,^,,t) . As the vector TDA, A{x,(,,t) can be expressed as a sum of the contributions coming 
from the active li-quark and the active d-quark. The first one will be proportional to the d's charge 
and the second to the u's charge, 



(x, t) = Qd al^d (x, t) + a^_^,, (x, t) 



(4.33) 



Isospin relates these two contributions, 

al^a{x,^,t) = a{x,^,t) , 

where the minus sign is originated in the change in helicity produced by the 75 operator. The 
expression for aj^^^ (x, ^, t) depends on the sign of ^. 

^ The pion pole structure of this contribution can be exphcitly understood if we observe that 



2ig 



2gnps {t) 



4£ 1 ~ 

-8Nc Ugnqq m — l2,A{x,^,t) 

I — t, — t 



(4.31) 
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In the 1^1 < X < 1 region we have 



a {x, ^, t) = -8 Nc fn g-Kqq m | |^1 + 

1 1 + e 1 



mlix-^) + {l-x)tl + ^ 



2 



hv {x, ^, t) 



+ 



log 



(mf — z rni^ 



m? (m^ — zm. 



with the abbreviations z = {1 — x) + x) / {1 + ^)^. And in the — |^| < a; < |^| region we have 
a {X, t) = -8 AT, U 9.qq m | (l + ""'^l"^ j"^ ' ^3. (x, 



(4.34) 



+ 



+ 



1 - C 167r2 



log 



i=l 



1 + e 



log 



2^ml - t (1 + (4^2^nf - xt) {m^ - yml) 



"^m? - xt) (4^ 



t — ml (4^2^2 _ a; ^) (4^2^2 _ 



xmt 



(4.35) 



where x = (^^ — x^) and y = z for ^ > and y = for ^ < 0. 



4.3.3 The Chiral Limit 



The expressions for both the vector and axial TDAs in the chiral limit, i.e. m^r = 0, are well 
defined. In particular, for t = 0, we find the following simple expression for v{x,^,t) 



;(x,e,0) = V2/,6F^^^(0) 



x + \C\ 
2|C|(1 + |C|) 



+ e(x- 1^1)0(1 -x)^ 



(4.36) 



where Fy ^{0) is the chiral limit of the vector pion form factor at zero momentum transfer 



F^^^(O) = lim^ 



0.17 GeV"^ 



A similar expression is obtained for a{x,^,t) in the chiral limit for t = 



aix,^,0) = -V2U6Ff^{0) 



e{e-x')e{0 



4^(1 + 1^1) 



x + ^ + {x-0 



+ e(x-iei)^(i-x) 



(l-x)(x-0 

(i-n(i-0J ' 



(4.37) 



with F'^ '^(0) the axial Form Factor at t = in the chiral limit 

Fl'^{0)= ]im Ff{0)/m^ = Ff^{0) 

m,r— >0 
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Figure 4.9: The vector Form Factor for the tt'' compared to the experimental data and to the dipole parameterization 
based on it [1U5) . 



4.3.4 Sum Rules and Polynomiality 



The vector TDA of the 7r+ must obey the sum rule given in Eq. (4.15) with the expression of the 



vector pion Form Factor in the NJL model given by Eq. (D.l). We have numerically recovered the 
sum rule for different t values. In particular we obtain the value -Fy^(O) = 0.0242 for the vector 
Form Factor at i = 0, which is in agreement with the experimental value 



Fy(0) = 0.017 ±0.008 



given in |llj . 



The tt" distribution must satisfy the following sum rule |165j 

dx (g„ Vfix, t) - Qd Vf{x, t)) =V2U F^^*^{t) 



(4.38) 



A theoretical prediction of the vr*^ Form Factor value is given in Ref. [39]. In particular, at t = 0, the 
value Ft^^*^{^) = 0.272 GeV~^ is found. The neutral pion Form Factor is directly related to the vr"'" 
vector Form Factor so that the sum rule is satisfied. We obtain the value Fjr-y*-y(0) = 0.244 GeV~^. 
The difference with the well-known result from the anomaly is due to the regularization parameters, 



see (D.4). In Ref. |105j . a dipole parameterization based on experimental data is proposed for the 
t-dependence of F-,^^*^ (t), 



F, 



7r7 7 



{t) 



7r7*7 



(0) 



1 - i/A2 



obtaining for the dipole mass A = 0.77 GeV. In Fig. |4.9| are shown the experimental results of the 
CLEO experiment [l05j (green diagonal crosses) obtained for higher values of t than our calculation 
allows. 1^ On the same plot, the red crosses represent our results for the neutral pion Form Factor. 

^Notice that the recent measurements of the BABAR collaboration for > 10 GeV^ show surprising results [2]: 
the measured neutral pion Form Factor exceeds the asymptotic limit /,r/Q^ and contradicts most of the models for 
the pion DA. 
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We have found that, for small values of t, the NJL neutral pion Form Factor can be parametrized 
in a dipole form with A = 0.81 GeV. The blue line represents the fit with the latter value for the 
free parameter. 



The axial TDA obeys the sum rules given by Eq. (4.16) with the axial Form Factor given in 
the NJL model by Eq. (D.2). This sum rule is satisfied for different t values. The numerical results 
also coincide. In particular we obtain FJ^(O) = 0.0239 for the axial Form Factor at t = 0, which 
is about twice the value 

F^+(0) = 0.0115 ±0.0005 , 
given by the Particle Data Group |llj . 



We expect TDAs to obey the polynomiality condition. However no time reversal invariance 
enforces the polynomials to be even in the A-momenta like for GPDs. That means that the 
polynomials should be "complete", i.e. they should include all powers in ^, 

,1 n— 1 

dxx''~'V{x,C,t) = Y,Cn,i{t)C ■ (4.39) 



-1 



1=0 



We have numerically tested the polynomiality and obtained it. We observe that, for the vector 
TDA, the coefficients for the odd powers in ^ are of one order of magnitude smaller than those 
for the even powers in ^. In particular we have numerically proved that, in the chiral limit, the 
polynomials only contain even powers in for any value of t, 



(4.40) 



We have analytically found that the odd powers in ^ go to zero for the polynomial expansion of the 
vector TDA. A study of the polynomiality in the limit given in Eq. (4.36) leads to the following 
analytical expression for the coefficients C^2i(0 



Cl2^m = V^U 2<^(0) (-1 + 2(-l)"-i) 



n 



n + 1 



(4.41) 



Notice they do not depend on i. The coefficients in the chiral limit, Eq. (4.40), as well as the 



coefficients for m^r = 140 MeV, Eq. (4.39), we numerically obtained are given in Table 4.1 



The ^-dependence of the moments of the axial TDA also has a polynomial form. Those poly- 
nomials contain all the powers in i.e. even and odd, 

»! n— 1 

" (4.42) 

i=0 



„1 n-1 

/ dxx^-'A{x,c,t) = Y,c'r,/t)e 

•J — 1 .—n 



An analytic study of the polynomiality in the limit given in Eq. (4.37) confirms that all the 
3 be present. The analytic values for the cc 

C'n^iO) = V2U 2Ff ^(0) (l + 2(-l)"-i) 



powers in ^ have to be present. The analytic values for the coefficients in this specific limit are 

n — 2i 



C^5m(0) = 2Ff ^(0) (l + 2(-l)""i) 



n(n + l)(n + 2) 

-2(i + l) 
n(n + l)(n + 2) 



(4.43) 
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which is in agreement with those numericaUy obtained. The coefficients of the polynomial ex- 
pansions in the chiral limit are very close to those obtained for the physical values of the pion mass. 
In Table |4.2l the coefficients C'^^{t) are therefore shown only for m.,^ = 140 MeV. 



The polynomiality property of the term containing the pion DA can also be studied. The t- 
dependence only comes from the pion pole and it is divergent at the chiral limit when t = 0. We 
have numerically found 



dx X 



n-l 



4fi 1 + e 



t 2i 



2i 



1 



n-l 



(4.44) 



which is true for both m-,^ = and 140 MeV. The even moments in x disappear due to the pion 
distribution amplitude is an even function in the interval [— ■?,^]. 

Calling the l.h.s. of Eqs. ( 4.42||4.44 ), respectively, A„(^, t) and <I'n(C) ^'^^ using the alternative 
tensorial decomposition (4.18), we tend to write a general relation about the polynomiality property 
of the whole axial bilocal matrix element, in the fashion of the alternative sum rule (4.20). In the 

= gauge, it reads 



dx X 



e £ 



e £ 



n—l 



dz 



|vr(]5^)) 



p-1 



27r 

• A^)-^[^„(e,t)-(2e)a>„(e,t)^ 

1 1 



(4.45) 



i=0 



We can therefore write a general relation about the polynomiality property of the whole axial 
bilocal matrix element 



dx X 



n-l 



-1 



dz 



(7(^7)^ 



k(p7r)> 



z+=z^=0 



e [£ 



1 " 



(4.46) 



with the highest power in ^ given by the pion pole contribution. The latter result is the expected 
one from Lorentz invariance, e.g. for the GPDs (1.26). 



4.4 Symmetries 

In the previous Section we have defined the TDAs in the particular case of a transition from a 7r+ to 
a photon, parameterizing processes like hadron annihilation into two photons and backward VCS 
in the kinematical regime where the virtual photon is highly off-shell but with small momentum 
transfer t 



HH —7- 7*7 — )• e^e 7 and — t- Hj 



(4.47) 
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n 


1 


2 


3 




C[,o{t) 




Cioit) Ciiit) Cl,{t) 


Wj, = 1-10 








t = 
t = -0.5 
t = -1.0 


22.4 
16.1 
13.2 


-3.77 4.00 
-2.97 2.57 
-2.53 1.97 


6.8 -4.9 2.3 
5.7 -3.5 1.7 
5.0 -2.8 1.3 



Table 4.2: CoefEcients of the polynomial expansion for the axial TDA. The pion mass is expressed in MeV and t is 
expressed in GeV^. Notice that the coefficients have to be multiplied by 10~^. 



with H a hadron. Symmetries relate the latter distributions to TDAs involved in other processes. 
For instance, we could wish to study the 7-7r~ TDAs entering the factorized amplitude of the 
process 

M being either pi or vr; in the same kinematical regime. 

By consistency with the processes that we are considering, we choose to change the definitions 
of the kinematical variables. The momentum transfer is A = p^^ — p^, therefore p^ = 171^/2 — t/A 
and t = A^. The skewness variable describes the loss of plus momentum of the incident photon, 
i.e. 



(4.48) 



and its value ranges between 



-1<C< 



{2ml -t) ' 



(4.49) 



Since the skewness variable defined for tt- 7 transitions is minus the one defined for 7— vr transitions, 
we should hence consider the results corresponding to the results for negative ^ is a previous Section. 



Within these conventions, we define the 7-7r='= TDAs 
dz' 



I 



27r 



z+=z-^=0 



I 



dz 
~2^ 



—i eSyC ^ Pp {p-K — Pj)a ^s/2f 



Z-L=0 



= ±- 



P 



-e[s^.{p^-p^)) 



+e (e • (ptt - P7)) ^(^) 



— t 



V2U 
2C 



(4.50) 
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With the help of Appendix |B| we find the following relations. Time-reversal relates the 7r"'"-7 
TDAs to 7-7r+ TDAs in the following way 

D^^^^x,^,t) = D'r^^^ix,-^) , (4.51) 

where D = V,A. And CPT relates the presently calculated TDAs to their analog for a transition 
from a photon to a vr" 

D-^^^'< {x,i,t) = Z)T^--(-x,-e,t) . (4.52) 



4.5 Pion-Photon TDAs in the Nambu-Jona Lasinio Model 



In Figs. 4.10 and 4.11, the vector and axial vr"*" — t- 7 TDAs obtained in the NJL model are plotted 
in function of x for several values of ^ and t. In these figures, the u-quark contributes to the 
TDAs in the region of x going from — |^| to 1 and the d-antiquark going from —1 to |^| . Therefore, 
for X £ [|^| ,1] {x G [—1, — |^|]) only n-valence quarks (d-valence antiquarks) are present (DGLAP 
regions). Besides, TDAs in the — |^| < x < |^| region (ERBL region) receive contributions from 
both type of quarks. For the vector TDA, we observe in Fig. |4.10| that the position of the maxima is 
given by the ^ value, separating explicitly the ERBL region from the DGLAP regions. The vector 
TDA is positive (negative) for negative (positive) values of x with the change of sign occurring 
around x = 0. This change in the sign of the vector TDA is originated in the presence of the 



electric charge of the quarks in Eq. (4.23). 



The process involved in the calculation of the TDAs allows negative values of the skewness vari- 
able. In the chiral limit, the skewness variable goes from ^ = —1 to = 1, for any value of t. In the 
chiral limit, the vector TDA for a negative value of is equal to the vector TDA for |^|. This can 



be seen from the polynomiality expansion in this limit, Eq. (4.40), which has only even powers of ^. 



For the physical value of the pion mass, negative values of ^ are bounded by t/ (2m^ — t) < ^. For 
each allowed value of ^, we found that the numerical value of V (x, — |.^| , t) is close to V (x, |^| , t) , 



due to the smallness of the coefficients of the odd powers of ^ in the polynomial expansion (4.39) 
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m =140 MeV and r = GeV 



Analyzing the axial TDAs, plotted in 
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Figure 4.11: 
respectively, 



The axial TDA in the case m,r 
» t = and (6) t = -0.5 GeV^ 



1,0 



140 MeV and, 



Fig. 4.11, we observe two different behav- 
iors depending on the sign of ^. For posi- 
tive ^, the position of the minima is given 
by the value of the skewness variable while 
the position of the maxima is always x ~ 
in the ERBL region and ib(l + ^)/2 
in the DGLAP regions. For negative ^, the 
value of the axial TDA at x = it^ is impor- 
tant and in some cases is a maximum and, 
in the ERBL region, t) presents a 

minimum near x = 0. As we have previ- 
ously shown, the axial TDA in the ERBL 
region receives contributions from two dif- 



ferent diagrams, depicted in Fig. 4.8 In the 
second of these diagrams, a virtual quark- 
antiquark interacting pair in the pion chan- 
nel appears. The pion pole, contained in 
this diagram, has been subtracted but the 
remaining non-resonant part contributes to 
the axial TDA. We observe that the latter 
contribution is the dominant one and pro- 
duces the maxima around x = for pos- 
itives i (Fig. |4.16| ). Now, the axial TDA 
does not change the sign when we go from 
negative to positive values of x. In the ax- 
ial TDA the change of sign originated in 
the presence of the electric charge of the 



quarks in equation (4.33) is compensated 



by the change of sign between quark and 
antiquark contributions generated by the 
75 operator present in the axial current. 



By comparing the plots in Figures 4.10 



and 4.11| for different values of the momentum transfer, it is observed that the amplitudes are lower 
for higher (— t) values, as it can be inferred from the decreasing of the Form Factors with (— t), 
connected to the TDAs through the sum rules. In the case of the vector TDA, by increasing the 
(— t) value, not only the width and the curvature of the TDAs are changed but we also observe that 
higher values of ^ are preferred, i.e. the sign of the derivative of the collection of maxima changes 
passing from a zero momentum transfer to a non-zero one. 

Isospin relates the value of the vector and axial TDAs in the DGLAP regions, 



V{x,i,t) = -\v{-x,i,t), A{x,i,t) = \A{-x,i,t) 



\S,\<x<l 



(4.53) 



being the factor 1/2 the ratio between the charge of the u and d quarks. We observe in Figs. 4.10 
and 4.11 that our TDAs satisfy these relations. It must be realized that the relation (4.53) cannot 
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Figure 4.12: The functions (a) V^" 
and for t = -0.5 GeV^ 



{x, 5, t) and (b) A''^'" {x, ^, t) for different values of the skewness variable ^ 



be changed by evolution. 

Regarding the chiral limit, we observe that both the vector and axial TDAs do not significantly 
change going from a non-zero pion mass to the physical mass, except for the change in the lower 
bound of ^. 

In Fig. 4.12 are illustrated the vector and axial TDA for the transition 7 vr", whose expres- 



sions have been given in (4.52) 



4.6 Polynomiality and Double Distributions 

In the previous Section, the definition of the pion-photon transition has been given from the tensorial 
structure of the pion radiative decay. This structure for the hadronic matrix element leads to the 
identification of the vector and axial TDAs as well as the pion pole contribution. 

Generalized Parton Distributions can also be studied through their moments |120j . In this 
Section, we build up the vector and axial-vector pion-photon transition distributions by considering 
matrix elements of towers of twist-two operators. This manifestly Lorentz-covariant approach leads 
us to define the Transition Distribution Amplitudes in terms of Double Distributions. 

The Mellin moments in x lead to derivative operators between the light-like separated fields. 



Such derivative operators correspond to the local twist-2 operators defined in Eq. (1.23) that, sand- 
wiched between two hadronic states, give generalized Form Factors. 



The pion-photon matrix elements of vector twist-2 operators can be decomposed in a fully 



Lorentz covariant fashion in terms of various twist-2 Form Factors C„fc(t) similarly to Eq. (1.24) 



adapted to the tensorial structure (4.2). Namely [199] 



i e 



(n- 1)! 



k=0 



k\{k-n-l)\ 



^ \ f^n — k 

2" 



2" 



(4.54) 



70 



Transition Distribution Amplitudes 



where the action of {. . .} on Lorentz indices produces the symmetric, traceless part of the tensor. 
The polynomiahty property immediately follows. 

The generalized moments can be written as a Double Distribution W(/3, a; t) which is the 
generating function of the twist-2 Form Factors, 



Cnk{t) 



dp 



(1-1/31) 



da /3 



n~l—k „ k 



a''W{P,a;t) 



(4.55) 



It is exactly this step from the tower of twist-2 Form Factors as generated by the double distribution 
that ensures the fulfillment of the polynomiahty property when evaluating the distribution through 
this parameterization. By definition \173\ I174j . the matrix element of bilocal vector operator is 
written in terms of DD in the following way 



dp 



(1-1/31) 



dae-*^(P-^)+*"(^-^)/2yy(^,a;t) . (4.56) 



From the definition of the vector TDA (4.13) and the skewness variable (4.4), we find the reduction 
formula 



dp 



(1-1/31) 



da5{x- p-ia)W iP,a;t) . 



(4.57) 



The analysis of the axial matrix element is quite similar to the analysis of the vector one. 
Nonetheless, it is worth reminding that the structure dependent tensorial form is not the only 
contribution. There is, actually, also the pion inner Bremsstrahlung to be separated out from the 
TDA. As already mentioned in Section 4.1 we adopt the tensorial decomposition (4.3) given in 
Ref. 



, I147| . required in order to ensure electromagnetic gauge invariance is pointed out. 
Thus, an appropriate tensorial decomposition for the definition of the Double Distributions 
would be 



{7{Pj)mO) T- 7^'^75 ^V^'^ ...i P^"-i> V'(0)k+(p.)) 



e (e- A) V2U 



■A' 



t — m: 



2p 



n-l 



n-1) 

(n-l)! 



V2U k\{k-n-l)\ 



+ 



Bn{t) 



A^ ^1 
2" 



A 

2" 



Mn-l} 



A 

2" 



f^n — k 



A 

2" 



(4.58) 



where we have chosen to decompose the second term on the r.h.s. in such a way that it does 
explicitly not depend on the gauge choice. Also, since there is a large freedom in the way of 

^In Appendix |f] the polynomiality property is derived from the expression of the Double Distribution in the 
a-representation of the propagator. 
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writing down this decomposition, we have chosen to isolate the highest power of ^, in the third 
line. However, the generalized Form Factor Bn{t) could perfectly be included into C^^(t). This is 
more naturally seen if we adopt the tensorial decomposition ■ 

Time-reversal invariance does not imply any symmetry on either term. The polynomiality 
property immediately follows. 

The decomposition of the axial twist-2 matrix element (4.58) leads to a relation in agreement 
with the general relation Eq. (4.45). From the polynomiality property in the NJL calculation, 
we know that the pion pole term yields the highest power of ^ for odd moments in x (4.44). In 
particular, in the NJL model calculation, no 5-like term is present. 

The moments of the pion DA are expressed in terms of 



dx I X 



1 



(j){x) 



with n even; what corresponds to the expression (4.44) to which we have extracted the i-dependence 
to keep trace of the pion pole. 

The generalized moments can be written as Double Distributions following (3.23). These DD's 
are the generating functions of the twist-2 Form Factors. The term depending on both p and A is 
expressed as W'(/3, a; t). 



(i/3 



(1-1/31) 



da (3 



n~l—k „ k 



a" W'{/3,a;t) 



(4.59) 



On the other hand, we have extracted the pure A-dependent generalized moments whose expression 
is therefore 



Bn{t) 



daa''-^B{a,t) 



(4.60) 



-1 



After all these considerations, we can write down the light-like separated matrix element 



(7 (pt)! 



e{e-A) V2f^ 



A-z 



+ -U-.A- 



t — 
1 



a + 1 



V2U 



2p ■ z 



A-z 



dl3 



-(1-1/31) 
dae'''^^-'^/^B{a,t) 



-1 



where the pion resonance structure has been subtracted. 



(4.61) 



From the definition of the axial TDA (4.14), we find the reduction formula 

1 r^-m 
dp 

-1 J-(1-\I3\) 



A^^{x,C,t) 



da6{x-P- Ca) W'^^ {f3, a,t) + sgn{C) B { 



(4.62) 
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The Fourier transform of the first term on the r.h.s. of Eq. (4.61) manifestly reduces to the pion 



pole term parameterization of Eq. (4.14) as expected 



A D-term was introduced in the matrix elements defining the GPDs. This new term is nec- 
essary in order to modify the tensorial structure to reproduce the correct behavior with the two 
independent variables p • z, A • z. We therefore interpret the D-ievm in the context of TDAs as the 
minimal contribution that we are compelled to introduce in order to restore gauge invariance in the 
tensorial structure adapted from the axial vector hadronic current involved in the pion radiative 
decay (4.3). In other words, the equivalent of the D-term is the pion pole contribution in the 



language of Double Distribution parameterization. 



4.7 Pion-Photon TDAs in Other Approaches 



Other studies of pion-photon TDAs have already been undertaken [1991 IHl I127j . We also mention 
the phenomenological approach, coming from the GPD analysis, used in Ref. |131| . All these 



calculations have been performed using the Double Distributions described in Section 3.2.1 



In this Section, we aim to compare the (x, ^, t) dependencies as obtained by different methods 
and models in order to give a critical point of view of our results. 



4.7.1 The a-Representation of the DD for the TDAs 

In this Section, we consider the a-representation of the propagators in the kinematics of the pion- 



photon transition described in Section 4.2 This technique leads to the identification of the Double 
Distributions defined here above. Nevertheless and as it is been highlighted before, the model 
imposes the integral decomposition as well as the relevant diagrams to be considered. The forth- 
coming considerations aim to be general and applicable to the model calculations whose results will 
be given in a further Section. 

In the vector current case, the trace over the Dirac indices yields to the exact tensorial struc- 
ture ( |4.13[ ). As a consequence, the 3-point function corresponds exactly to the vector TDA for 
the M-quark up to an overall factor g coming from the trace over isospin, flavor, color and Dirac 
spaces. 1^ In terms of Double Distribution, it reads 



v{x,(,,t) = gi3{x,(,,t) 



df3 



(4vr)- Jo 

with the Double Distribution in the same kinematics being 

^(1 -a-/3) 



da 5 {x — (3 — aS,) Isil^, a, t) 



(4.63) 



9 



m2l/3(l-/3 



a) - 1(1 + a - m 



a 



(4.64) 



A detailed derivation of this result is given in Appendix |Fj The total DD for the vector TDA 
VV(/3, a, t) is obtained applying the previous result to the isospin decomposition (4.23). The support 



property follows from the /3-integration (F.ll, F.12) 



*In the very particular case of the NJL model calculation, the factor is g = 8Nc f-n 



e.g. Eq. 14.241 
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From the DD (4.64), one can easily see that the hmit — )• imphes that 

Is{(3,a,t) = X^{l3,-a,t) . 

This so-called a-symmetry reflects the appearance of only even powers of in the polynomiality 
expansion of the vector TDA; a result that has been obtained in the calculation of skewed PDs in 
the NJL model ( |440l ). 

The trace over the Dirac space for the axial current involves 2- and 3-propagator integrals. 
Besides the terms proportional to the external momenta, the trace evaluation of the diagram of the 
type (a) of Fig. 4.8 yields a factor —2e-kpn- This renders the calculation slightly more complicated 
but the /^''-dependent integral can be reduced to the same 3-propagator integral as for the vector 
one. 

The evaluation of the Double Distributions is similar to the one undertaken for the 3-propagator 
integral of the vector TDA. The subsequent steps are as described in Appendix |F} and the result 
for the /c'^-dependent integral is 



e-/3(a;,^,t) 



g- A 

(4^ 



(i/3 



-1+1/31 



da5{x- (3-aO (1 + q - /?) IsiP, a, t) 



(4.65) 



with the same double distribution X3(/3,Q,t) defined by Eq. (4.64). 



4.7.2 Phenomenological Parameterizations of TDAs 

Phenomenological parameterizations of the (x, ^, t) dependence of the parton distributions consti- 
tute a good alternative and/or complementary approach of the diagrammatic analysis of the Double 
Distributions. In what follows, we report a fully phenomenological t-independent parameterization 
and a semi-kinematical approach improved by phenomenological considerations. 

It is easier to start describing the (x,^) dependency by setting t = 0. Hence, we will, in a 



first time, consider the t-independent Double Distributions, entering the reduction formulae (4.57 



4.62). 



In the DVCS framework, the suggestion has been made by Radyushkin \n5\ 1176] to parame- 
terize the GPDs in such a way that the forward limit, i.e. 

1-1/31 

daT'^{f3,a) = q{f3) , 

■1+1/31 

is ensured. The function q{P) is hence related to the parton distribution for the flavor q. Basically, 
the phenomenological forward quark distribution as measured from DIS can be used as an input 
in the parameterization. 

The TDA do not obey any forward limit, but it has nevertheless been proposed to apply the 
DIS phenomenology as a first approach |131] . Hence, for the TDAs, we will use the same form as 
for the GPDs, i.e. 



^'?(/3,a) = h{(3,a)q{l3) 



(4.66) 
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for T'^ corresponding either to the vector or axial DD, namely W or W' respectively. The function 
a) denotes a profile function, that can be parameterized through a one-parameter ansatz, 
following Refs. [I75l [T761IT77] . 



/i(/3, a) 



r(26 + 2) [(1 



The profile function is normalized to 1. 



a 



226+ir2(?, + i) (1 



26+1 



(4.67) 



The parameter h in (4.67) characterizes 



0.025 



0.015 



0.005 



0.000 




Figure 4.13: The TDA v'^ as d escrib ed in Ref. [HI] through the 
ansatz for the profile function (4.671; for values of the ^ variable 



the strength of the ^-dependence of the 
resulting GPD. The limiting case h — )■ 
oo corresponds to the ^-independent 
ansatz for the GPD, namely Hix^ ^ = 
q{x). The power in 6 is a free parameter 
for both the valence and sea distribu- 
tions. They can be used as fit parame- 
ters in the extraction of GPDs from the 
experiments. For instance, the twist-2 
DVCS predictions of Ref. |201] corre- 
sponds to the choice 



b = by 



al 



1 



0.2,0.4, 0.6 & 0.8, respectively, the plain red, the dashed purple, the This mild dependence on ^ is fol- 
dashed-dotted orange and the short-dashed cyan curves. lowed by the authors of Ref [I3lj for 

the parameterization of the pion-photon 
TDA. 



Also, Lansberg et al. |131| assume, as a first guess, that the /3-dependence of q for the 7 — )• vr 
transitiorj^ is given by a simple linear law|^ 



?(/3) 



2(l-/3)0(/3) 



The t-dependence comes solely from the Form Factors through the sum rules. Therefore, the 
resulting distribution d(x,^), with d = v,a, must be normalized to 1 when integrated over x. This 
leads to a TDA whose dependence on the three kinematical variables is given by 



Dix,^,t) 



(4.68) 



D denoting either the vector or axial quantities. 

In the same Reference the t-dependence coming from the Form Factors is considered negligible 
as we are concerned with small t region. In this way, we can use the values given at zero momentum 



transfer, see Section 4.3.4 



^See in Section 



4.4 



how to go link this transition to the 



The parameterization is here given for the quark distribution, the antiquark would be g(/3) = —q{—P). 
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As explained in Section 3.2.1, the reduction from the DD to skewed distributions requires a 



particular dependence on the skewness variable, imposed in the ^- function of the reduction for- 
mula. This choice leads to the identification of two kinematically different regions. The result 
obtained through the above-given ansatz is shown, for the u-quark and for a vector current, on 



Fig. 4.13, The shape of the vector u-quark TDA reasonably follows what we obtained in the NJL 
model. Fig. 4.10 The same ansatz is used for the axial TDA with the simple replacement D = A 



in (4.68). Therefore the difference between the vector and axial TDA resulting from this ansatz 
resides in the ratio F^/Fy and the charge combination. The total 'j-tt~ distribution does, in turn. 



change according to the charge combinations involved in the transition, see Section 4.4 Also, this 
parameterization from the phenomenology of GPDs does not include the term compensating for 



the gauge invariance in (4.61). However, this contribution is closely related to the pion pole whose 



importance in fulfilling the sum rule for the axial current has been highly emphasized in the pre- 
vious Sections. We conclude that the axial TDA proposed by |131j considerably and conceptually 
differs from the one obtained in the NJL model. 

Another (semi-)phenomenological approach including, this time, the phenomenology of photon 
distributions has been proposed in Ref. [199J. The t-dependent double distributions defined in 
Eqs. (4.56, 4.61) are calculated in a simple quark model. In doing so, Tiburzi determines the 
t-dependence of the distributions and reproduces the Form Factors. His intermediate results, in 
agreement with the results of Ref. [S], are Fa(0)/Fv(0) = 0.98 and Fa{0) ~ 0.026. However the 



partonic content of the model has been modified to give reasonable phenomenology while evolved. 
But positivity bounds for the pion-photon Transition Distribution Amplitudes are investigated 
instead. The positivity constraints are an application of the Schwartz inequality in the Hilbert 
space [80l 11641 11661 11671 11681 1176j . Basically, in the [—1, — ^] and 1] regions, we can find a bound 
of the type, for the pion GPD, 



where the momentum fractions of the quark before and after the scattering are 



(4.69) 



i + e 



^out 



X — ^ 
1^ 



The author of Ref. |199] provides for some estimates of the TDAs on the basis of positivity 
bounds. The latters, determined using the technique described in Ref. |164j . involve the u-quark 
and d-antiquark distributions, the first one related to the pion and the second to the photon. 

The Double Distributions obtained in the simple quark model are modified in order to saturate 
those constraints. The modifications of the model are motivated by the apparition of pion and 
photon parton distributions in the new realistic double distributions. Since the evolved PDFs 
are well-known, this allows one to establish the renormalization scale of the model avoiding evo- 
lution. For the pion, a simple analytic parameterization of the valence and sea quark has been 
chosen following Ref. [99]. The real photon distribution is parameterized in a vector- meson dom- 
inance fashion |100j- From the perspective of model building, this comes out to be an attempt 
to satisfy the known constraints and generate an input distribution at a moderately low scale. 
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Vectorial TDA. (m„=2 m and t=-0.5 GeV^ ) 




Thus the values of the vector and ax- 
ial Form Factor, as given by |TT], are used 
in order to determine the empirical con- 
stants of the parametrization. In the same 
way, the constituent quark mass parame- 
ter m is obtained (fitted) requiring that the 
model pion-photon transition Form Factor 
comes close to the experimentally parame- 
terized form. The value m = 0.2 GeV is ob- 
tained. The resulting TDAs, illustrated on 



-1.00 



-0.50 



0.00 



0.50 



1.00 



Fig. 4.14, do not satisfy the isospin relation 



(4.53) due to the different choice of the u- 



Figure 4.15: 
GeV^ 



The vector TDA for = 2m and for t — —0.5 



quark and d-antiquark distributions used in 
the saturation of the positivity bounds. We 



have studied the positivity bounds (4.69) 



for GPDs in the NJL model and noticed that it is actually an upper bound that is sometimes very 
higher than the value of the GPD itself. 

Moreover the TDAs shown by Tiburzi are rather peaked at x = ^, what is not reproduced in 
the NJL model calculation. Rather, at weak binding when = 2m, the latter presents peaks at 
2(1 — £,) and X = 5(1 + ^) for the vector TDA Fig. 4.15, This is what one expects from a 



X = ^[L - ^) and X = 2(1 + for the vector TDA Fig. 4. 
free quark picture. This has already been observed for the GPDs in the NJL model [153] and in a 
previous work of the same Tiburzi |194j . Given the previous argument, the position of the peaks 
obtained by Tiburzi is not understandable in the NJL model. 



4.7.3 Model Calculations of the TDAs 

Other model calculations of the pion-photon TDAs have already been done, respectively, in the 
Spectral Quark Model (SQM) [S] and a non-local Chiral Quark Model (xQM) [HT]. The Double 
Distributions defined through the a-representation of Section [4.6| have been used in both calcula- 
tions. We here after comment and compare their results. 

The vector and axial TDAs calculated in the SQM, NJL model and non-local xQM are com- 



pared in Fig. 4.17 
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The authors of the first reference use the asymmetric notation (F.8, F.7). The comparison is 
here awkward since the authors define, by turns, 



C = {P-y -P7r)^/pt 



< = ip. 



Ptt) 



n 



We nevertheless decide to use the standard relation between their asymmetric notations and the 
symmetric ones (F.8) [81]. Their functions Vsqm and ^sqm correspond, respectively, to the v and 



a given in Eqs. (4.23 



4.33). They are obtained through the reduction formulae, e.g. (4.57) for the 



vector TDA. As for the axial TDA, the decomposition (4.61) is not exactly adopted: Broniowski 
et al. structurally isolate the pion pole by neglecting the whole A • z structure instead. In so doing 
they select from the trace only the piece proportional to {e'^ ■ A-^ )p^ of the tensorial decomposition 



dz 



(7 (Pt)!' 



|)riT5r-g(^ 



K(p7r)) 



p-\ 



V2U 



where, here, A is seen as A (4.19) to which the pion resonance has been structurally subtracted. 



Therefore the piece proportional to (e • A)A'^ leads to both L'-term like structures, one being the 



resonant pion pole (4.14), and is represented by the ellipses. 



The normalization condition is different from the one used in the NJL model calculation and 
the results quoted by these authors must be multiplied, for the vector TDA, by a factor 



48^2 V2/.^^(Q) 



10 



before comparison. From Fig. 4.17 we conclude that there is a qualitative and quantitative agree- 
ment, for the vector TDA, between the results of Ref. |44j and those obtained in the NJL model. 
Regarding the axial TDA we observe, in addition to the normalization factor 



48vr2^FA(0) 



10 



a change in the global sign due to different definitions. In the first version of |44] the axial TDA 
for positive values of ^ is giverj^ It coincides with the results in the NJL model, as observed in 
Fig. 4.17 Surprisingly, the result presented in Ref. |44j coincides with our result for negative ^. It 



is perhaps due to the change in the definition of ^ mentioned above. Rather, one might call into 
question the selection procedure of the structure dependent quantities. 



In Ref. |127j the TDAs are calculated in three different models. The first one is a local model 
whose pole structure has some similarity with the one of the NJL model. The two other models, i.e. 
semi-local and fully non-local, follow the results of the local one. The same structure decomposition 
as in the SQM is followed. The most prominent difference between the results obtained in the non- 
local and the NJL models is the appearence of important odd powers in ^ in the polynomial 
expansion of the vector TDA. We know from Ref. |155] that, for non-local models, there are 

^There is a typographic error in Eq. (23) of the first version of this reference, where a factor My/Q must be 
dropped. 
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additional contributions to those calculated in Ref. [127j . In the case of PDFs, disregarding these 
contributions can produce small isospin violations |156j . It can therefore be considered that, on 
that point, the results of Ref. |127j must be confirmed. 

For numerical comparison, the results obtained in |127j must be corrected by a factor 



2vr2^Fy(0) 



0.45 



due to the use of a different normalization condition. We observe, see Fig. 4.17 that our results in 
the NJL model coincide with the results of the latter reference. [3 



m =140 MeV and t = -0.5 GeV 




Regarding our axial TDA, it is 
worth noticing that it differs from 
the previous calculations [199^ SU 
I127j due to the effect of the non- 
resonant part of the second dia- 
gram of Fig. |4.8[ This contri- 
bution, corresponding to the last 



Figure 4.16: Contributions to the axial tt"*" — >■ 7 TDA for both positive 
= 0.25, plain line) and negative = —0.5, dashed line) values of 
the skewness variable and for — 140 MeV and t — —0.5 GeV^. In 
each case, and in the ERBL region, the contribution coming from the 
first diagram of Fig. |4.8| is represented by the dashed-dotted lines and 
the non-resonant part of the second diagram of Fig. |4.8| is represented 
by the dotted lines. 



term of Eq. ( |4.35 ), is proportional 
to {t — m^) but with zero value 
for the residue. The presence of 
this term is crucial in order to ob- 
tain the axial Form Factor using the 
sum rule as shown by the reduc- 
tion formula (4.62). Furthermore 



this term is dominant in the ERBL 
region as we can infer from Fig. 

ESI 



It can eventually be concluded 
that there is no disagreement be- 
tween the different studies concern- 
ing the 7r-7 TDAs besides the ambi- 



guity in the definition of the skewness 
variable. Calculations of Refs. [44 1 1127] are performed in the chiral limit, where ^ runs from —1 to 
1. The symmetric nature of this interval makes difficult to check the sign of ^. On the other hand, 
the NJL calculation is given for the physical pion mass. In this case, the kinematics of the process 
imposes t/{2ml -t) <(, <l. From Eqs. ( |4.34| |4.35[ ) we observe that there is a pole in the axial 

t/(2m^ — t), preventing us from going through unphysical values of ^. 



t)<C< 

TDA for the limit value ^ = 

Moreover, the sum rules (4.15, 4.16) for both the vector and axial TDAs are here satisfied for 
physical values of ^, and broken in the unphysical regions (, < t/ (2m^ — t) and ^ > 1. We therefore 
conclude that the choice of sign in our calculation in the NJL model is consistent and gives a 
guideline for comparing with other models. 

* For the axial TDA, there is a change in the definition of the skewness variable between the caption of Figs. 3, 5 
and Fig. 9 of Ref. [127) . The agreement is valid if the convention of the caption of their Fig. 9 is chosen. 
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■ - - 5=0.5 Non Local 




Figure 4.17: Comparison of the 7r+ 7 TDAs for t = -0.1 GeV^ of Refs. glim] for = MeV and Ref. [6Zl 
for tiIt^ = 140 MeV. On the left, we have: the vector TDA for ^ = ±0.5 as a single (solid) curve for the results of 
both the NJL model and SQM (these four curves are indistinguishable); the result of the non-local xQl^ calculation 
for ^ = 0.5 (dotted line) and for ^ = —0.5 (dashed-dotted line). On the right, we have the results of the axial TDA 
with the choice for the sign of ^ as discussed in the text. 



4.8 Discussion about the Results 

In this Chapter we have defined the pion-photon vector and axial Transition Distribution Ampli- 
tudes using the Bethe-Salpeter amphtude for the pion. In order to make numerical predictions we 
have used the Nambu-Jona Lasinio model. The Pauli-Villars regularization procedure is applied in 
order to preserve gauge invariance. 

We know from PCAC that the axial current couples to the pion. Therefore, in order to properly 
define the axial TDA, i.e. with all the structure of the incoming hadron being included in A (x, ^, t), 
we need to extract the pion pole contribution. In so doing, we found that the axial TDA had two 
different contributions, the first one related to a direct coupling of the axial current to a quark of 
the incoming pion and to a quark coupled to the outcoming photon and a second related to the 
non-resonant part of a quark-antiquark pair coupled with the quantum numbers of the pion. 

The use of a fully covariant and gauge invariant approach guaranties that we will recover all 
fundamental properties of the TDAs. In this way, we have the right support, x G [—1,1], and 
the sum rules and the polynomiality expansions are recovered. We want to stress that these three 
properties are not inputs, but results in our calculation. The value we found for the vector Form 
Factor in the NJL model is in agreement with the experimental result ^llj whereas the value found 
for the axial Form Factor is two times larger than in [llj. This discrepancy is a common feature 
of quark models [44J. Also the neutral pion vector Form Factor is well described. These 

results allow us to assume that the NJL model gives a reasonable description of the physics of those 
processes at this energy regime. 

Turning our attention to the polynomial expansion of the TDAs, we have seen that, in the 
chiral limit, only the coefficients of even powers in ^ were non- null for the vector TDA. This result 
is quoted as the recovery of the so-called a-symmetry of the Double Distributions. No constraint 
is obtained for the axial TDA. Nevertheless, the NJL model provides simple expressions for the 
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coefficients of the polynomial expansions in the chiral limit, Eqs. (4.41) and (4.43). 



We have obtained quite different shapes for the vector and axial TDAs. This is in part, at least 



for the DGLAP regions, imposed by the isopin relation (4.53). We have pointed out the importance 
of the non-resonant part of the qq interacting pair diagram for the axial TDA in the ERBL region. 
It is interesting to inquire about the domain of validity of the isospin relations ( |4.53 ). These 



relations are obtained from the isospin trace calculation involved in the diagram (a) of Fig. 4.8 
Due to the simplicity of the isospin wave function of the pion these results are more general than 
the NJL model and could be considered as a result of the diagrams under consideration. These 
diagrams are the simplest contribution of handbag type. 

Also, in the last Section, we have presented the different approaches proposed to gather informa- 
tion about the pion-photon transition distribution amplitudes. The model calculations are in good 
agreement one to others. On the other hand, the phenomenological approaches are either based 
on GPD's phenomenology or pion and photon PDF's. It is unfortunately the best we can do for now. 

The Transition Distribution Amplitudes proposed by the authors of [1(55] open the possibility 
of enlarging the present knowledge of hadron structure for they generalize the concept of GPDs for 
non-diagonal transitions. Calculated here, as a first step, for pion-photon transitions, these new 
observables lead to interesting estimates of cross section for exclusive meson pair production in 7*7 
scattering |131j as we will see in the next Chapter. 




Figure 4.18: Baryonic p tt" TDA. 



Transition Distribution Ampli- 
tudes have also been extended to 
baryonic transitions |133l I134[ 11351 
[I65j . For instance, the nucleon- 
to-pion TDAs describe exactly how 
a baryon can turn into a meson, 
namely the transition from a bary- 
onic to a mesonic state. The leading 



twist TDAs for the p — t- vr' 
tion, depicted in Fig. 4.18 
fined from the correlator 



transi- 
are de- 



{^\p^)\e'i^ut{zin)u^^{z2n)ul{z^n)\p{p^,si)) , 
the latter matrix element being also related to the baryonic DAs. 



There are 8 leading-twist TDAs for the p — )■ vr" transition, namely, two vector V^^'^(xj, ^, A^), 

two axial [xi,^, A^) and four tensor Tf" (xi, ^, A^) ones. The peculiarity of the baryonic TDAs 
is that, for both the meson and photon cases, 3 quarks are exchanged in the t-channel, satisfying the 
relation X1+X2 + X3 = 2^ (^^ > 0). The interplay of the variables implies 3 kinematical regions; one 
DGLAP region, Xi > 0, and two ERBL regions, xi > 0; X2 > 0; X3 < and xi > 0; X2 < 0; x^ < 0. 



A calculation in the Meson-Cloud model of Ref. |158j has been performed in the ERBL re- 
gion \159\ 11621 1163] , where the pion cloud contribution -as coming from a higher order in the Fock 
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state expansion- plays an important role and could therefore be tested. It is so far the only non- 
perturbative approach to these new distributions. 

The experimental importance of baryonic TDAs will be shortly overviewed in the next Chapter. 
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Exclusive meson production in 7*7 scattering 



This Chapter is devoted to the applications of the results for the photon-pion Transition Distribu- 
tion Amplitudes displayed in the previous Chapter. In particular, we discuss the results of Ref. [70] . 



Cross section estimates for the processes 

ill TT+vr^ , p+vr" , (5.1) 



factorized according to Fig. 4.1 have been proposed in Ref. [131] using for the TDA the t- 



independent double distributions (4.66), in a first approach, and, in a second, the i-dependent 



Double Distribution (4.64) obtained through diagrammatic analysis. 



To implement these analyses, in this Chapter, we display the cross section for the processes (5.1 ) 
with the large-distance part calculated in models. In a previous Section we have concluded that 
there is a clear agreement between the different model calculations of the pion-photon TDAs, we 
are allowed to analyze the result of a single model, e.g. the NJL model. Furthermore, this choice 
allows us to keep trace of the pion pole contribution in order to include it in our analysis. 

5.1 Processes and Kinematics 

The 7*7 — )• M+7r~ process, with Af+ = or 7r+, is a subprocess of the 

e{Pe)+l{p^)^e{p',)+M+{pM) + TT- {p^) , (5.2) 
process. We follow all the definitions of the kinematics given in Ref. |131j . with the exception 



that our vector is defined following (A.9), which is twice the vector used in Ref. [131 J. In 
particular, for massless pionsj^we have. 



^In Section III.A and Fig. 3. 

■^The pion mass is set to — MeV throughout the Chapter. 
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Figure 5.1: Kinematics for the processes ( |5.1[ ) [131) . 



^^' = (9 + ^7)' 

Se7 = (Pe +^7)^ 



{Pe-p'eY 



(5.3) 



with pe and p'^, respectively, the momen- 
tum of the incoming and outgoing elec- 
tron; 



= (1 + Op : 



Pn = il-OP + 



A 



±2 



2(1-0 



n + A-L, 



(1 + 4)P H : —n 



1+c 



(5.4) 



where At = (0, A-'-,0), and therefore A|, = — A^^. Notice that 

A^2^(-t)(l-0/(l + 6 , 

with t <0. 



It is worth noticing that, even if the skewness variable is not constrained to be positive but 



rather to belong to ^ G [—1,1], the external kinematics of the process (5.1 ) restrict ^. Effectively, the 



kinematics (5.4) have been given with the photon momentum expressed in terms of the momentum 



transfer squared and the center of mass energy squared of the 7*7 system W'^, through which 
we can determine the skewness variable, i.e. 



y 



P-y -Pe 



2i s, 
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The momentum transfer therefore reads 



-2Cp + -^n 



(5.5) 



The terms contributing to the cross section are selected by multiplying the amplitudes by the 
longitudinal polarization of the incoming photon £l^. An expression for el^ is found knowing that 



ej = 1 and that e^.g = 0. For q as given by (5.5), we find 



2e_^ Q 



(5.6) 



following Ref. |131j for the convention for the overall sign. The real photon polarization is defined 
by the condition e~ = together with the gauge condition e"*" = 0. The kinematics are illustrated 
on Fig. 5.1 
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5.2 Cross Sections 

The differential cross section for the kinematical variables Q^,t,$, and the angle ip between the 
hadronic, where takes place the subprocess — )• M+7r~ with M being either pj^ or tt, and 
leptonic planes reads 

= ^ ^ I A .e7^eM+7r-|2 /rr 

dQ^dtdid^ 2,2 {2ttY sl^ + ' ■ ^ ^ ^ 

The square amplitude for the process 67 — )■ eM+vr^ is obtained through the square amplitude for 
the subprocess 7^7 — t- M^vr" and the contribution of the fermionic line from which the highly 
virtual photon is emitted 

|A^^^-^^"-p = ^A^^r{i>'J^{q)i>Jl{q))\M''^--\' . (5.8) 



Using Eq. (5.6), we find that the trace yields 



And the amplitude for the subprocess is the factorized amplitude of Fig. 4.1, It reads 

MPi^~{Q\i,t) = - [\x [\z^f-<f>piz)M}^{z,x,OV{x,C,t) , (5.9) 



-M""""(Q',C,t) = I dx I dzMz)M^{z,x,0 

-1 Jo 



-1 JO 

12 



(5.10) 



where z is the light-cone momentum fraction carried by the quark entering the M"*" and fp = 0.216 
GeV. The hard parts are defined as in |131j 

Mn^,x,o ^ 8.w.|^.--p.(P.-P.).^(^^+^^) , 



2 Cf (e^.A^) 1 f Qu , Qa 



Mh (Z, X, = 47r aelmOis -IT, 7^ 7Z r — — + 



Nc Q z{l — z) \x — ^ + ie X + ^ — ie 



(5.11) 



The pion DA is chosen to be the usual asymptotic normalized meson DA, i.e. (pniz) = 6z(l — z), 
what cancels the z-dependence of the hard amplitude. The integral over x remains which contains 
the mo del- dependent TDA. Separating the real and imaginary parts of the amplitude, but only for 
the structure dependent part, this integral is defined 



J-l x-(, J_i x + ^ 

+Qu Die, 1 1) - ^vr) + Qd D{-i, t) (in^ + ivr) , (5.12) 
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P7 
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(b) 
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P'l 
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Figure 5.2: Pion Compton Scattering. The direct and the cross contributions, respectively (a) and (b), as well as the 
seagull contribution (c) . The dot vertex is the usual 7r7r7- vertex and the square vertex is the Brodsky-Lepage vertex. 



where D stands for V, A. This non-perturbative part will be given by the TDAs evaluated in the 
NJL model of Section lO 



As for the pion pole contribution, the expression for the amplitude (5.10) suggests that the pion 



amplitude is defined by performing the entire integral procedure just as for the TDAs. The pion 



pole term becomes proportional to the electromagnetic pion Form Factor studied in Section [277 

1 3 Q^F^iQ^) 



r 



Qz (1 — z) with z = (x + ^) /2^, we obtain the Brodsky-Lepage pion Form Factor (2.40) 



Therefore the pion pole contribution to the amplitude (5.10) becomes 



-47r aelm 2 



t — ml 



Q 



(5.13) 



as 47r t — 

In this case we would use the asymptotic form for the pion DA by consistency. If we use (pT^ {z) 



(5.14) 



(5.15) 



In an alternative way, the analyses of the pion Compton scattering yields the same result. For the 
sake of generality, we choose, for the purpose of that calculation, the general expression for e given 
by e = + ap. Since it does not depend on a, we can state that the result obtained is independent 
of the gauge's choice. Therefore, in the a = gauges, only the direct diagram (a) of Fig. 5.2 does 



contribute to the cross section. Both the crossed (6) and the seagull (c) diagrams are higher-twist. 
Contrarily to the amplitudes involving the TDAs, the pion amplitude is model-independent thanks 
to the asymptotic expression of the pion DA that will be henceforth adopted. 



An analysis of Eq. (5.11) including the isospin decomposition (4.23 4.33) shows that only the 
direct terms of the charge combinations, i.e. tx Qg, contribute to the amplitude M^l'" . A similar 
analysis for the pion production but at zero momentum transfer indicates that the cross terms, i.e. 
oc QuQd, cancel out. 



Going back to the differential cross section Eq. (5.7), we average over the real photon polariza- 
tion, which give an overall A^/2 factor. The integration over ip being straightforward, we are left 
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(5.17) 



From Eq. (5.17) it can be observed that 

^>qV(2s67-q') 



(5.18) 



In other words, there is a (positive) lower limit on the value of ^. Let us now proceed to evaluate 
these integrals. 



5.3 Results from the Phenomenological Parameterization 

The first estimates for the cross sections of meson production have been given for the phenomeno- 



logical parameterization of the Double Distribution (4.66), what is referred to as Model 1. The 



t-dependent Double Distribution (4.64) are considered for the vector TDA and is called Model 2. 
In none of these approaches the pion pole contribution to the cross section for the pion production 
is taken into account. It means that the term proportional to the pion decay constant is not present 



in the expression for the cross section (5.17). 

Also, no QCD evolution is taken into account, the effects of which are supposedly less important 
than the uncertainty of the modeling of the TDAs. 

For the vector TDA both the real and imaginary part of the amplitude contribute significantly 
to the cross section. Since we are dealing with t-independent DDs, the t-dependence comes solely 
from the overall —t factor in the trace. 

In order to numerically estimate the cross sections, the strong coupling constant is fixed to 
a value of ~ 1. This large value of is indicated from the analyses of the electromagnetic 



Form Factor [38j. On the left panel of Fig. 5.3, the cross section for the vrp production is plotted 



for = 4 GeV^ and Se7 = 40 GeV^. The high-^ behavior is similar for both model while the 
intermediate-^ region is sensitive to specific models for the TDAs. On the other hand, as it is shown 
in Ref. |131j . the behavior at fixed ^ is model independent, what constitutes a crucial test for 
the validity of the approach. 

Turning our attention to the vrvr production, we should compare the estimates with the cross 
section for the Bremsstrahlung process in the same kinematical regime. It is shown in Ref. [131J 
that the process where the tt+tt" is produced by a photon radiated from the leptonic line does not 

factor 1/4 is missing in Eq. (23) of Ref. [131) . This typo does not affect to the numerical results reported in 
that reference 11361. 
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0.2 0.4 0.6 0.8 1 

Figure 5.3: The cross section for up production (on the left) in a t-independent parameterization (Model 1) and a 
i-dependent one (Model 2) ; and for tttt on the right [131j . 



interfere much with the tttt production from 7*7 for reasonable values of ^ and with the exception 
of small values of Se7- In other words, the hard hadronic process dominates the Bremsstrahlung 
contribution in the kinematical region under scrutiny. This statement is illustrated on the right 



panel of Fig. 5.3 The axial TDA is given for the Model 1 and the dotted line represents the 



contribution of the Bremsstrahlung. 

A general observation is that both the cross sections for p^vr^ and vr+vr^ are rather small. In 
the next Section, we will see how these results are enhanced by the pion pole amplitude for the 
pion production as well as by QCD evolution for the rho production. 



5.4 Results from the TDAs in the NJL Model 

For the nonperturbative part of the process we use the TDAs evaluated in the NJL model. 



The kinematics of the process restrict the domain of relevance of the variables which, in 



turn, define the shape and the magnitude of the distributions. The restriction (5.18) on the value 



of the skewness variable is indeed particularly interesting because the value of ^ defines the shape 
of the TDAs. In particular, the shape of the axial TDA calculated in models radically changes 
according to the sign of the skewness variable; A{x, ^, t) has, at x = it^, its maximum values for 



^ > (see right panel of Fig. 4.12) while it has its minimum values for ^ < 0. However the vector 
TDA has its maximum and minimum values aX x = (see left panel of Fig. 4.12) independently 
of the sign of the skewness variable. 

On the other hand, the magnitude of the distributions is controlled by the t-dependence. This can 



be easily understood because the TDAs, that must satisfy the sum rules (4.15, 4.16), are expected 
to decrease at least as The i-behavior of the cross sections will therefore be dictated by both 
the overall (— i) factor coming from the trace and the TDAs. 

The value ~ 1 chosen in Ref. |131j is also adopted here. Using this value for we have 



evaluated the cross section for p production. In Fig. 5.6 we plot this cross section as a function of ^. 



As we observe, the cross section is largely dominated by the imaginary part of the integral (5.16). 



5.4 Results from the TDAs in the NJL Model 
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Comparing with the previous results shown on Fig. 5.3, we observe that our predictions are higher 
by a factor 2 or 3. 




e (GeV ) 



Figure 5.5: Electromagnetic Pion FF. The Brodsky-Lepage asymptotic 
FF for Qs = 1 and — 2.45 GeV^ (blue line) and the experimental 
data given in Table [5?T| The red curve represents a fit in the monopole 
form QV(1 + 0.44(370.2336). 



The cross section for pion produc- 
tion at = 4 GeV^ function 
of ^ is given in Fig. |5.7| (left). This 
cross section is dominated by the 
pion pole contribution which is deter- 
mined by the Brodsky-Lepage pion 
Form Factor. Alternatively, if the 
pion Form Factor is experimentally 
known, we can infer phenomenolog- 
ically this contribution. And hence 
the axial TDA could be extracted 
from the interference term. The lat- 
ter shows already an impressive en- 
hancement with respect to the cross 
section shown in the right panel of 



Fig. 5.3 



From Ref. |109j we know that the 
pion Form Factor at = 2.45 GeV^ 
is 0.167 lb 0.010. The pion Form Fac- 



tor is shown on Fig. 5.5 for the dif- 
ferent experimental results, see Ta- 
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Figure 5.6: 67 — > e' p'[n~ differential cross section plotted as a function of ^ for Q = 4GeV , s^-, = 40GeV , 
t = —0.5 GeV^ . In the first layer, the dotted (dashed) line is the contribution to the cross section coming from the 
imaginary (real) part of the integral given in Eq. (5.161. In the second layer we give the cross sections before and 
after evolution. 



ble 5.1, In Fig. 5.7 (right) we depict, for each contribution, the prediction using the interval 
defined by the experimental value of the Form Factor (filled areas), including also the theoretical 
prediction using the Brodsky-Lepage pion Form Factor (lines). We observe an unexpectedly huge 
enhancement of the total cross section of a factor of about 60. 



The t-dependence of the cross section for pion production includes a strong dependence on 
t coming from the pion pole. Neglecting the pion mass in (5.17), we observe that the pion pole 
contribution to /,r is proportional to t~^. Therefore, the cross section grows as for small t values. 
For large t values we expect, as in the p production decreasing as t . 





Q' [GeV] 


F"(Q") 


F''{Q'^)+ error 


F'^iQ'^)- error 


DESY [Tj 


0.351405 


0.211314 






DESY [37J 


0.695093 


0.329482 


0.353265 


0.306529 


JLab |193j 


0.599788 
0.748309 
1.00094 
1.52456 


0.259397 
0.254932 
0.313972 
0.34985 


0.274914 
0.280412 
0.345017 
0.397251 


0.246048 
0.230928 
0.284536 
0.309278 


JLab |109j 


1.59404 
2.44972 


0.385976 
0.405233 


0.408247 
0.434364 


0.368385 
0.38488 



Table 5.1: The experimental results for the pion Form Factor. The last two columns represent the upper and the 
lower limit of the error bars. 
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Figure 5.7: 67 — >■ e'n'^-K differential cross section plotted as a function of ^ for Se^ = 40GeV^ t = -0.5GeV^ 
= 4GeV^ (left) and = 2.45 GeV^ (right). The dashed (dashed-dotted) [dotted] line is the contribution to the 
cross section coming from the pion fom factor (axial TDA) [interference term between the pion Form Factor and A] . 
The pion pole contribution is calculated using the Brodsky-Lepage pion Form Factor. The filled areas in the right layer 
correspond to the same contributions but with the experimental value for the pion Form Factor = 0.167 ± 0.010 
[1091. 



5.4.1 QCD Evolution 

We have also studied the effect of the QCD evolution on our estimates, using for this purpose the 
code of Freund and McDermott [94,j • The same procedure as for the Distribution Amplitude and 



the GPDs is followed. The value of Qq for the LO evolution in the NJL model is given by Eq. ( 2.30 ). 



Turning our attention to the vector TDA evolved at LO to a scale of 4 GeV^ (Fig. 5.4), we 
observe that the value of V {x, ^, t) at x = it^ grows for small values of ^ and decreases for large ^ 
values in comparison with the TDA at the scale of the model. This implies that the cross section 
for p production, which is largely dominated by the imaginary part, will grow appreciably in the 



small 1^ region. In Fig. 5.6 we compare the cross section after evolution, calculated only through 
contribution of the imaginary part of with the one obtained before evolution. We observe that 
this cross section is multiplied by a factor about 5 in the ^ ~ 0.2 region. In the case of the axial 
TDA, the cross section is dominated by the pion Form Factor contribution, therefore the effect of 
the evolution is expected to be small. 



5.5 The Experimental Situation 

In the previous Section we have looked at the expected cross section for vr-vr and vr-p production in 
exclusive 7*7 scattering in the forward kinematical region using realistic models for the description 
of the pion. First we confirm the previous estimates for p production, even if our results for the 
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7*(g) P'{P2} 




P{v\) Av^) 

Figure 5.8: The factorization of the process 7'P — )■ P'tt" into proton DA, hard subprocess amplitude Mh and 
proton-pion TDA. Figure from [135) . 

cross section are a factor 2 larger than the one obtained in Ref. |131j . Second, in comparison with 
this previous evaluation of the cross sections, we have improved in considering the effect of evolution 
on the vector current. In doing so, an additional factor 5 appears in the small ^ region, leading 
to a cross section for p production of one order of magnitude larger than the previous calculation. 
We have also improved in including the pion pole term in the tensor decomposition of the axial 
current. Then an even larger enhancement factor, of about 60 in this case, is found in the cross 
section for pion production. The interference term becomes a factor 15 larger than the pure TDA 
contribution, making the axial TDA more accessible experimentally. 

Data have been collected at LEP and CLEO on these reactions with some phenomenological 
success, even if mainly on pp final states in a different factorization regime |16[ I17j. More data are 
obviously needed and are eagerly waited for in the TDA region, and much hope comes from the 
high luminosity electron colliders. 

5.6 Backward Pion Electro-production 

Recently Pire and collaborators have shown in a series of articles |165[ 1133) 1134) I135| that the 
TDAs represent the ideal framework through which it is possible to describe the backward pion 
electro-production, e.g., — >• e'p'ii^ , as well as the meson production via a nucleon-antinucleon 
Drell-Yan process, e.g., NN — )• 7*7r, in the forward kinematics. 

These processes are experimentally more accessible then the ones previously analyzed. For 
instance, the pp annihilation, e.g. pp — )• 77, the meson production channel pp — t- 77r^, that will be 
intensively studied at GSI-HESR antiproton program by the planned FAIR project, |^ and by the 
meson photo-production process, P7 — )• tt^p, with forthcoming data from JLab. 

In Ref. |135| , a soft-pion limit for the baryonic TDA is used to give an estimate of the differential 
cross section for 7*^—7- P'tt^. However, the lack of knowledge of the TDAs unfortunately prevents 
from giving an accurate estimate of the latter quantity. 



*A great wealth of up-to-date informations about the experiments that will be performed at GSI-HESR can be 
found at this web-page: http://www.gsi.de/zukunftspro jekt/indexe.html. 
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Data should be available, at higher energies, as for instance at HERMES and with CEBAF 
at 12 GeV |178j . The calculation of the cross section |135j . in an admittedly quite narrow range 
of the parameters, can thus serve as a reasonable input to the feasibility study of backward pion 
electroproduction at CEBAF at 12 GeV, in the hope to reach the scaling regime, in which we are 
interested. 
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Peeping through Spin Physics: the Sivers Function 




Ph 



X 



The partonic structure of transversely polar- 
ized hadrons is one of their less known fea- 
tures. A review [27| on the transverse polar- 
ization of quarks in hadrons gathered all the 
knowledge which is, nevertheless, appreciably 
growing on both the theoretical and experi- 
mental sides. Experiments for the determi- 
nation of the transversity properties are pro- 
gressing very fast and the relevant experimen- 
tal effort has motivated a strong theoretical 
activity. 

When extending our analyses to unintegrated 



Figure 6.1: Semi-Inclusive Deeply Inelastic Scattering. 



parton distributions (1.28), we have to deal 



with the intrinsic transverse motion of quarks. 
This supplementary degree of freedom implies 
the existence of new parton distribution func- 
tions, i.e. (1.28). Besides their dependence on kx, the additional parton distributions can either be 



chiral-even or odd. In the former case, the distributions can be probed in fully inclusive DIS. The 
conservation of chirality does however not allow the identification, through the optical theorem, of 
chiral-odd distributions, for helicity must be flipped twice as shown on Fig. |6.2[ The Transverse 
Momentum Dependent Parton Distributions (TMD) studied in the context of transversely polarized 
nucleons do flip helicity and cannot be accessed through DIS. 



Semi-inclusive deep inelastic scattering (SIDIS), i.e. the process depicted on Fig. 6.1 

e{l)N{V)^e'{l')h{Vh)X{V^) , 



(6.1) 



95 



96 



Peeping through Spin Physics: the Sivers Function 




with the detection in the final state of a pro- 
duced hadron h in coincidence with the scat- 
tered electron e', is one of the proposed pro- 
cesses to access the parton distributions of 
transversely polarized hadrons. For sev- 
eral years it has been known that SIDIS off a 
transversely polarized target shows azimuthal 
asymmetries, the so called single spin asym- 
metries (SSAs) IHTl MfTOO] , 
As a matter of fact, it is predicted that the 
number of produced hadrons in a given direc- 
tion or in the opposite one in the hadronic 

plane, with respect to the reaction plane, depends on the orientation of the transverse spin of a 
polarized target with respect to the direction of the unpolarized beam, 



Figure 6.2: Representation of a chirally odd distribution. 
The corresponding handbag diagram is forbidden by conser- 
vation of chiralhty. 



^UT{(ph,4>s) 



d(y{4>h, 4>s) - dai(f>h, (t>s + tt) 
d(T{(l)h, (As) + da{(ph, 4>s + vr) 



(6.2) 



It can be shown that the SSA in SIDIS off transverse polarized targets is essentially due to two dif- 
ferent physical mechanisms, whose contributions can be technically distinguished \24: \ I35 iri281ll49j . 



One of the two different mecha- 
nism has been proposed by Collins 
|61j and is therefore referred to after 
his name. It is due to parton final 
state interactions in the production 
of a spinless hadron by a transversely 
polarized quark, and will not be dis- 
cussed here. 

The other is the Sivers mecha- 
nism, producing a term in the SSA 
which is given by the product of 
the unpolarized fragmentation func- 
tion with the Sivers Parton Distribu- 
tion [1891 1190] , describing the mod- 
ulation of the number density of un- 
polarized quarks in a transversely po- 
larized target due to the correlation 
between the transverse spin of the target and the intrinsic transverse parton momentum. 




Figure 6.3: SIDIS kinematics in the 7*p frame as defined in the Trento 
conventions 1231. 



'^In this case, we can define two different transverse directions. In general, a quantity which is transverse in a 
frame where P and h have no transverse components (thus Pt = hx ~ 0), is indicated with a subscript T, so that 
T means transverse with respect to h, while those with a subscript _L are defined in a frame where q and P have no 
transverse components (standard DIS frame), so that _L means transverse with respect to q (thus q± = P± — 0) [149) . 
There is no such distinction in DIS, so that T =_L and they can be used equivalently. 

^ Here U means Unpolarized beam, T means Transversely polarized target and it is assumed that the produced 
hadron is spinless, or that its polarization is not detected. 
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The Sivers function, f^{x, kx), is the quantity of interest of this Chapter. It is formally defined 
through the $[^^1 of the decomposition (|l.28|). Namely, [23lfTn2] . 



1 f d^-(fiT 
(27r)3 



with the operator defined by jl21l | 



(6.3) 



(6.4) 



(6.5) 



The number density of unpolarized quarks in a transversely polarized target is eff'ectively mod- 
ulated by the quantity 



= V5g(0, r , iT)Cl^ (oo, r)7+'Co(oo, 0)V,(0, 0, 0) , 
where ipq{C) is the quark field and the gauge link has been defined as follows, 



-^.-•(00,^ ) = Pexp -ig I A^{r] ,CT)dr] 



<l>[^^]^(x,fcT,5, =t)-'^^^''^'ix,kT,S, =i) 



(6.6) 



with i,j being either x,y and 5 =t / i standing for S in a transverse direction, in opposition with 
the helicity S = ±. 



The Sivers function contributes with a sin(0/j — 
4>s) weighting function to the cross section of 
SIDIS with unpolarized beam off a transversely 
polarized proton (6.2), with the angles cph, 4>s de- 



fined in Fig. 6.3 The Sivers asymmetry, accord- 



ing to the "Trento convention" [23], is defined in 
terms of the experimental cross sections. 



A Sivers 



J dcpsdcph sin((/>fe - (j)s)d'^cruT 
f dcpsdcphd^auu 



(6.7) 



The cross sections themselves can be written, 
through the factorization theorems, as a convo- 
lution of fc-r-dependent distribution and fragmen- 



tation functions [149] as shown on Fig. 6.4; the 
former being here the Sivers function. 




Figure 6.4: Factorization in SIDIS. 



Recently, the first data of SIDIS off transversely polarized targets have been published, by 
the HERMES Collaboration [5] for the proton and by the COMPASS Collaboration [7] for the 
deuteron. It has been found that, while the Sivers effect is sizable for the proton, it becomes 
negligible for the deuteron, so that apparently the neutron contribution cancels the proton one, 
showing a strong fiavor dependence of the mechanism. 
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(a) (b) (c) (d) 

Figure 6.5: The expansion of the gauge link: (a) + (6) + (c) + (d) + . . .. 



Different parameterizations of the available SIDIS data have been published \19\ l63l [20] , still 
with large uncertainties. Deeper analyses are in progress, e.g., [TH [75]. New data will be soon 
available, e.g., thanks to the improvement of the CEBAF accelerator at 12 GeV at JLab [53] or the 
improvement of the experiments at HERMES [184j , that will reduce the errors on the extracted 
Sivers function and will help discriminate between different theoretical predictions. 



6.1 The Gauge Link 

The Sivers function is a Transverse Momentum Dependent PD and vanishes if integrated over the 
transverse momentum. It is a time-reversal odd object [27j and for this reason, for several years, it 
was believed to vanish due to time reversal invariance. However, this argument was invalidated by 
a calculation in a spectator model |41] . following the observation of the existence of leading-twist 
Final State Interactions (FSI) [32]. The current wisdom is that a non- vanishing Sivers function is 
generated by the gauge link in the definition of TMD parton distributions |62[ I12H [29] , whose con- 
tribution does not vanish in the light-cone gauge, as happens for the standard parton distribution 
functions. 



Technically, Final State Interactions are reflected in the expansion of the gauge link (6.5) as 



shown on Fig. 6.5 It is hence clear by the identification of the Sivers function in (6.6) that, if the 



Sivers function were defined without the gauge link, it would be identically zero. 

To the first non trivial order, i.e. (b) — (c) of Fig. 6.5 giving a contribution to the asymmetry, 
the Sivers function is obtained for the flavor a, as follows. 



flT{x,kT) 




(6. 



where, in a helicity basis, the operator reads 



+ h.c. . 



(6.9) 



^ The gauge link is usually taken in a non-singular gauge, i.e. not in the light-cone gauge, for further details we 
refer to Ref. [29]. 



6.2 Final State Interactions at the Quark Level 
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6.2 Final State Interactions at the Quark Level 

In this theoretical and experimental scenario, due to the lack of correct QCD calculations of Parton 
Distributions, it becomes important to perform model calculations of the Sivers function. Several 
estimates exist, in a quark-diquark model |41l 11211 [22l I25j : in the MIT bag model, in its sim- 
plest version |206j and introducing an instanton contribution |56) : in a light-cone model |142] : in 
a nuclear-style framework, relevant to establish the arising of the Sivers function in proton-proton 
collisions [32j : in the Constituent Quark Model (CQM) of Isgur and Karl (IK) [66j. 



Constraints on the Sivers function are obtained from first principles, following the example of, 
e.g., sum rules in GPDs. A model independent constraint on calculations of f^{x,kT) is the 
Burkardt sum rule [4:9\ I50j. It states that the total average transverse momentum of the partons in 
a hadron, {kx), which can be defined through the first moments of f^{x, /ct), for all the partons 
in the target, has to vanish. If the proton is polarized in the positive x direction, the Burkardt sum 
rule reads, 

J2 (k?) = - j dxj dkT^fi^{x, kT) = . (6.10) 

Q=u,d,s,g.. ^ 

The Burkardt sum rule remains unchanged under evolution. 

To distinguish between the model estimates, data and model independent relations, such as this 
one, can be used. 

It is also difficult to relate the Sivers Function to the target helicity-flip, impact parameter 



dependent (IPD), generalized parton distribution E (1.13). Although simple relations or ansatze 



between the two quantities are found in models [S^ I144j , a clear model independent formal relation 
is still to be proven, as shown in Ref. |146j . However some gross features, such as the sign of the 
Sivers function, can be deduced from the results for the IPD parton distributions. 



In the context of Constituent Quark Model for the proton, one expects, as a consequence of the 
Burkardt sum rule, the u and the d Sivers functions to contribute with a similar magnitude but 
opposite sign. 

In this Section, we analyze the Sivers functions in two such models; the MIT bag model and the 
CQM of Isgur and Karl. The gluon exchange is here treated in a perturbative way; both models 
moreover assume an 5f7(6) proton wave function, even if only in the first order of the IK wave 
expansion 

If an SU{6) state is assumed, any model calculation of Parton Distributions gives u and d 
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parton distributions which are proportional to each other [113^ IllSj . In the case of one is 
not investigating single particle properties but two-body ones, due to the two-body Final State 
Interaction operator generated by the gauge-link. There is therefore in principle no reason to 
have proportionality, which is a good sign that the results in such models might a priori fulfill the 
Burkardt sum rule ( |6.10[ ).p1 



The diagrammatic representation of the leading-twist Final State Interactions contributing to 



a non-zero Sivers function (6.8) in those models is depicted in Fig. 6.6 As it has been mentioned 
above, helicity has to flip twice: both at the nucleon state and at the quark level. At the quark 
level, it is the presence of a 2-body light-like gluon interaction that enables a spin-flip. By sym- 
metry among the active quarks, the quark interacting with the virtual photon (green quark line) 



can either flip its helicity or not. The first case is represented on the left panel of Fig. 6.6 while 
the second case implies that the " active-but-non-interacting- with-the-photon" quark flips in turn 
as shown on the right panel. 

Let us now concentrate on the results of these two model calculations. Historically the first 
calculation performed in the MIT bag model has been undertaken by Yuan |206] whose conclusions 
were followed by the whole Community. Our CQM calculation has led us to satisfactory Sivers 
functions [66] even if non complying with the conclusions of Ref. [206]. This is why we "peeped 
through" the first calculation in the bag model. A more physical picture of the Sivers function as 



allowing symmetry between the two active quarks 6.6 was revealed [72j. Therefore, the results of 



Refs. [66[ [72] are presented here as illustrations of 1 and 3-body approaches of the Sivers function. 



6.2.1 1-body Model 



Here we revise the first formal calculation 
of the Sivers function that was performed 
in the MIT bag model ^206j. The details of 
the calculation are given in Ref. |72j with 
the formalism developed in Ref. [206] . 



The expression (6.8) is written inserting 



the bag model wave function in momentum 
space, (/9(A;). 1^ Using Fig. 6.7 for the defini- 
tion of the quark helicity and momentum 
labels, the Sivers function in the bag model 




Figure 6.7: The quark helicity and momentum label for the bag 
calculation. 



*A proportionality of —1 would nevertheless be more than acceptable for the fulfillment of the sum rule. 
^ The expression for the wave function of the MIT bag model in momentum space is 



V5m(fc) 



\a-kti{\k\)x. 



(6.11) 



/ 3 \l/2 

with the normalization factor A'^ = ( , . , — | where uj = 2.04 for the lowest mode and Rn is the bag 
radius. The two last quantities are related through the relation RqMp = Auj. The functions ti{k) are defined as 



ti{k) ^ V? duji{ukRo)ji{uLu). 
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can be written [206j 

13 , mi ,1712 ,m3 ,1714 

^vLS^-?l-)l^l^^m2{k) j j^^L.ihh'^l^^miiks-q^) , (6.12) 

where M is the proton mass, Ep its energy, 6g ^ is the annihilation operator for a quark with flavor 
Q, helicity m, and color index i, and T?- is a Gell-Mann matrix. 



Performing the integral over k^ and then q± and calculating the appropriate flavor-color factors, 
we find that there is actually no proportionality between the u and the d quark distributions. 
Numerical results are shown in Fig. 6.8 for the first moment of the Sivers function (plain curves), 

= I cPkT^J^^^ixM) ■ (6.13) 



The obtained Burkardt sum rule is 



(fcS) = -0.78 MeV . 

Q=u,d,s,g.. 

To have an estimate of the quality of the agreement of this result with this sum rule, we consider 
the ratio 

r = TTiii 77^ - 0.05 , 



i.e., the Burkardt sum rule seems to be violated by 5 



In order to compare the results with the data, one should realize that one step of the analysis 
is still missing. The scale of the model is much lower than the one of the HERMES data, which 
is = 2.5 GeV^. For a proper comparison, the QCD evolution from the model scale to the 
experimental one would be necessary [116j. Unfortunately, the evolution of TMDs is still to be 
understood. In order to have an indication of the effect of the evolution, we perform a NLO 



evolution of the model results assuming, for the moments of the Sivers function, Eq. (6.13), the 
same anomalous dimensions of the unpolarized PDFs. |^ The parameters of the evolution have been 
fixed in order to have a fraction ~ 0.55 of the momentum carried by the valence quarks at 0.34 
GeV^, as in typical parameterizations of PDFs , starting from a scale of ;Ug ~ 0.1 GeV^ with only 
valence quarks. 

The results show an impressive improvement of the agreement with data once the results of 
Ref. |!206j are evolved in the way just described. The data, represented by the dashed bands, are 
perfectly described for both flavors. Although one should not forget that the performed evolution 
is not really correct, the analysis shows that, after evolution, model calculations can be consistent 
with the available data. 

^Important progress concerning the evolution of the first moment of the Sivers function have been made recently, 
see [1221 [205] , 
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The shght violation of the Burkardt sum rule, which is a momentum property, is probably due 
to the fact that the proton state is not a momentum eigenstate in the MIT bag model. 



6.2.2 3-body Model 



The Constituent Quark Model (CQM) has 
a long story of successful predictions in 
low energy studies of the electromagnetic 
structure of the nucleon. The Isgur-Karl 
model [TTUl llll| is chosen to perform a cal- 
culation, in order to describe the perfor- 
mances of the approach developed in the 
context of CQM in Ref. [66j. A difference 
with respect to calculations of PDFs and 
GPDs is that, in the calculation of TMDs, 
the leading-twist contribution to the one- 
gluon-exchange (OGE) Final State Inter- 
action has to be evaluated. This is done 




Figure 6.9: The quark helicity and momentum label for the 
CQM calculation. 
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through a non-relativistic reduction of the relevant operator, according to the philosophy of Con- 
stituent Quark Models [78], 



The expression (6.8) is developed by inserting proper complete sets of intermediate one-quark 



states, using translational invariance and the definition of the gluon field as explained in Ref. [66j 



plus the definition of the labels of Fig. 6.9, one is left with the following expression for the Sivers 
function, 



{PS. = MTtjTSl b^liki)bi'^\h)b''J,Ch - q)b^^fi, + q) V{h, h, q)\PS. 



with the interaction determined by 



.15) 



V{ki,k3,q) = -^Um'^{k3)j+u^,^{k3 + q)umi{ki)'y+u^,^{ki - q) 



(6.16) 



For an easy presentation, the quantity which is usually shown for the results of calculations or 



for data of the Sivers function is its first moment (6.13 ). The results of the CQM approach are given 



by the red curves in Fig. 6.10 for the u and the d flavor following the example of Fig. 6.8 in the bag 
calculation. They are compared with a parameterization of the HERMES data, corresponding to 
an experimental scale of = 2.5 GeV^ [HS]. The two dashed curves represent the best fit proposed 
in Ref. [63] plus 30 % and minus 30 %, respectively. Such an interval gives a rough estimate of 
the experimental errors and theoretical uncertainties in the extraction of the Sivers function from 
data. 



The evolved results are given by the blue curves in Fig. 6.10 for the u and the d flavors. As it is 
clearly seen, the agreement with data improves dramatically and their trend is reasonably repro- 
duced at least for x > 0.2. 



Within our scheme, at the scale of the model, the Burkardt sum rule is 



-0.40 MeV 



obtaining r ~ 0.02, so that we can say that our calculation fulfllls the Burkardt sum rule to a 
precision of a few percent. 



In this model the proton wave function is obtained in a OGE potential added to a confining harmonic oscillator; 
including contributions up to the 2&j shell, the proton state is given by the following admixture of states 



\N) = afSi/2)s + b'\^S[/2)s + c|^Si/2)m + d\*Di/2)M 



(6.14) 



where the spectroscopic notation \^^^^Xj)t, with t — A,M,S being the symmetry type, has been used. The 
coefficients were determined by spectroscopic properties to be a = 0.931, b — —0.274, c = —0.233, d = —0.067 |97) . 
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(b) 



Figure 6.10: The Sivers functions obtained in the IK model and compared with the extraction from the data |63|.The 
red curves represent the results before evolution, whereas the blue curves are after evolution to the scale of the data. 



6.3 The Way to a Conclusion 

With respect to Ref. ^206j, the results in the MIT bag model (Fig.|6.8|) show an impressive improve- 



ment of the agreement with data once the full contribution, i.e. coming from both diagrams 6.6 is 
taken into account. The data are described rather well for both flavors. Comparing this encourag- 
ing outcome with that of the IK model, one can notice that the Burkardt SR is better fulfilled in 
the CQM. 

In closing our analyses of the Sivers functions in a 1 versus 3-body model, we can say that, for 
the first time, it has been established that correct model calculations provide phenomenological 
successful interpretations of the Sivers function, which are consistent with each other. 



Let us compare our approach with other calculations mentioned above. 
The first version of the MIT bag model calculation p06 j has non-vanishing u and d-quarks con- 
tribution of opposite sign which are proportional in magnitude. The d-quark contribution is much 
smaller than the one obtained in the revised version of the bag model as well as the one of the CQM 
and therefore does not satisfy the Burkardt sum rule. The MIT bag model modified by instanton 
effects [56] has u and d-quark contributions of the same sign and therefore does not satisfy the 
Burkardt sum rule. 

The diquark model with scalar diquarks [22j has no contribution for the d-quark and therefore does 



not satisfy the Burkardt sum rule, Eq. (6.10). The diquark model with axial- vector diquarks has 



contributions to both u and d-quarks and with opposite sign, but with the magnitude of the d 10 
times smaller than that of the u. The Burkardt sum rule is not satisfied. The results are improved 
in the version of Ref. [25J . 

As a summary, we can say that the calculations presented here [66l ES] , despite the naive wave 
functions used, are in better agreement with the data with respect to the other approaches, and 
fulfill the Burkardt sum rule. The model calculations have still to be improved and their evolution 
properly undertaken. So have to be the extraction of the data as well as the experiments. Much 
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efforts in that direction is being put from the whole Transversity community. |^ 



We cite, e.g., the latest workshop "Transversity 2008: 2nd International Workshop On Transverse Polarization 
Phenomena In Hard Processes", 28-31 May 2008, Ferrara, Italy. 
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Our main source of information about the internal pion and nucleon structure is provided by Deep 
Inelastic Scattering. Unpolarized DIS has revealed the partonic distributions of the hadrons. It 
has shown that the quarks carry only about half of the total hadron's momentum. Polarized DIS 
has provided information about spin distributions, telling us that only a small fraction of the spin 
of the proton is carried by the intrinsic spin of the quarks. 

In the last decade, a theoretical treatment of semi-inclusive as well as exclusive deep reac- 
tions has been developed, increasing our knowledge on the structure of hadrons. For instance, 
Semi-Inclusive Deep Inelastic Scattering allows to obtain information also on chiral-odd Parton 
Distributions. In a transversely polarized hadron, the transversity distribution tells us what is 
the number density of quarks with polarization parallel to that of the hadron, minus the number 
density of quarks with antiparallcl polarization. Moreover, the complete study of the features of 
the produced hadron in SIDIS requires that we take into account the transverse motion of quarks. 
The theoretical tools for the description of this process are therefore the Transverse Momentum 
Dependent PDs. Among this class of distributions, the Sivers function is the object of important 
theoretical studies as it has been proposed to explain single-spin asymmetries in SIDIS. 

On the other hand, exclusive electroproduction of photons (Deep Virtual Compton Scattering) 
have been the most studied processes. Under the conditions of large virtuality of the exchanged 
photon and low momentum transfer t, i.e. near the forward direction, the amplitude of the reaction 
factorizes into a hard and a soft part. This soft part have opened new windows on the structure of 
hadrons. 

The theoretical tools for the understanding of theses processes are the Generalized Parton 
Distributions. The GPDs connect in different limits to other well known physical quantities: in the 
limit of low momentum transfer, t — > 0, they connect to the usual Parton Distributions. Integrated 
over the momentum fraction of the quark, GPDs lead to the electromagnetic hadron's Form Factors. 
Moreover, through Ji's sum rule, a particular combination of moments of GPDs can be related to 
the spin of the proton. 
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Exclusive hadron production in 7*7 scattering at small momentum transfer and large invariant 
mass is also assumed to factorize into hard and soft parts. The soft part provides for additional 
information on hadron structure through the Transition Distribution Amplitudes. 

In this thesis, we have considered some of the distributions that we have just described. 
Throughout the Chapters [2| |3j [4] & [5| the partonic structure of the pion has been scrutinized 
under the Nambu - Jona-Lasinio description of the pion. In Chapters [2] &: [3j we have reviewed the 
previous situation; whereas, in Chapters |4] & [5], we have presented our results on the Transition 
Distribution Amplitudes. 

In Chapter [2j the Distribution Amplitude and the Parton Distribution Function of the pion have 
been presented together with their relation through the Soft Pion Theorems. The study of these 
well-known quantities have provided for a way to set the scene for the subsequent and original part 
of the manuscript. To be specific, more than a simple illustration of the formalism in which we 
were to be working, the calculation of the pion DA has allowed for a first criticism of the choice 
of the NJL model. This has been done, among other things, by determining the scale of validity 
of this quark model Qq by evolving according to the QCD evolution equations. In particular, the 
Parton Distribution Functions of the pion have been found to be in a surprisingly good agreement 
with the data, once evolved to Leading-Order. We have also considered the Next-to-Leading- 
Order evolution, and we have reached the conclusion that the effect of the NLO evolution has 
been compensated at LO in the adaptation of the low model's scale Qq. On the other hand, the 
Distribution Amplitude of the pion calculated in the NJL model have been found to be in a good 
although less impressive agreement with the data. 

The conclusions that we can draw from the analysis of the pion PDF and DA also follow 
from the questions of the support property as well as the end-point behavior. First, we conclude 
that the conservation of covariance on the light-cone is a must-have feature of the model that 
one would have chosen for the analyses of distribution functions. We have actually learned that 
the support property is closely related to this symmetry, which is respected in the NJL model. 
Second, we emphasize the importance of perturbative QCD when going to higher energies in which 
the experiments are performed. We therefore advocate discarding some models only after having 
evolved their results to a perturbative scale. 

Given the symmetries it respects, and the good results it has led to, the NJL model, with all 
its merits and faults, can be considered as a realistic model for the pion not only for low energy 
observables, but also for the determination of the Parton Distributions. The emergent image of the 
pion is that it is basically determined by the chiral symmetry. 

Chapter |3] has been devoted to the Generalized Parton Distributions of the pion. In particular, 
the results of Ref. |153j have been displayed together with the parameterization of GPDs via Double 
Distributions [171j . The expected support and the polynomiality for the GPD of the pion have 
been obtained. There is no relevant discrepancy between the results obtained in the NJL model 
and the other calculations, e.g. jlU I171| . before and after QCD evolution. The agreement after 
evolution gives us confidence on the evolution code [9l] that we have used. 

In Chapter [4], we have principally presented our results for the pion-photon vector and axial 
Transition Distribution Amplitudes, respectively V{x,S,,t) and A{x,£,,t) \67\ l68l [69] . We have 



109 



used the formalism developed in Chapter [2] the Bethe-Salpeter amplitude for the pion has been 
constructed with the quark-pion vertex function coming from the Nambu-Jona Lasinio model. 
This technique allows for numerical predictions. Also, the Pauli-Villars regularization procedure 
has been applied in order to preserve gauge invariance. 



The definition of these pion-photon TDAs have been displayed in Section 4.2 The Partial 



Conservation of the Axial Current tells us that the axial current couples to the pion. Therefore, 
in order to properly define the axial TDA, namely with all the structure of the incoming hadron 
included in A (x, ^, t), the pion pole contribution must be extracted. By doing so, we have found that 
the axial TDA had, in the NJL model and the handbag approximation, two different contributions: 
the first one is related to a direct coupling of the axial current to a quark of the incoming pion and 
to a quark coupled to the outgoing photon. The second one is related to the non-resonant part of 
a quark-antiquark pair coupled with the quantum numbers of the pion. The presence of the latter 
contribution guarantees the gauge invariance of the axial TDA. 

The use of a fully covariant and gauge invariant approach guarantees that all the fundamental 
properties of the TDAs will be recovered. In this way, the right support, x £ [—1,1], has been 
obtained and the sum rules, together with the polynomiality expansions, have been recovered. We 
want to stress that these three properties are not inputs, but they are results springing from our 
calculation [67j. The value obtained in the NJL model for the vector Form Factor is in agreement 
with the experimental result provided in the Particle Data Group [TT], whereas the value found for 
the axial Form Factor is two times larger than the ones found in the PDG |llj . This discrepancy is 
a common feature of quark models as it has been obtained in different calculations |44[ I127j . Also, 
the neutral pion vector Form Factor is well described, reproducing the value given by the 

anomaly [39] . These results allow us to assume that the NJL model gives a reasonable description 
of the physics of those processes at this energy regime. 

Turning our attention to the polynomial expansion of the TDAs, it has been shown that, in 
the chiral limit, only the coefficients of even powers in ^ were non-null for the vector TDA p7] . 
No constrain has been obtained for the axial TDA. Nevertheless, the NJL model have provided for 
simple expressions for the coefficients of the polynomial expansions in the chiral limit (4.43). 

We have obtained quite different shapes for the vector and axial TDAs. At least for the DGLAP 
regions, this is in part due to the isospin relation 

V{x,^,t) = -W{-x,(,t), A{x,^,t) = ^A{-x,^,t), \(\<x<l . 

The difference in the shapes is also due to the non-resonant part of the qq interacting pair diagram 
for the axial TDA in the ERBL region. It is interesting to inquire about the domain of validity 
of the previous relations. These relations are obtained from the isospin trace calculation involved 
in the central diagram of Fig. 4.8 Due to the simplicity of the isospin wave function of the pion, 
these results are more general than the NJL model, and therefore more general than a result of the 
diagrams under consideration, which are the simplest contributions of the handbag type. 

Since there exist other model calculations as well as phenomenological approaches of the pion- 
photon TDAs [441 1127[ I13H [TUU] . a qualitative comparison and/or criticism of our approach was in 
order. This has been done in Section [4. 7[ where the Double Distributions for the TDAs have been 
defined. We have concluded that there is no disagreement between the different studies concerning 
the TT-'j TDAs, besides some ambiguities in the definition of the skewness variable. 



The TDAs, calculated, as a first step, for pion-photon transitions, have led to interesting esti- 
mates of cross section for exclusive meson pair production in 7*7 scattering. These processes involve 
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photon-pion transitions. Also, the symmetries relating the different transitions, e.g., vr"'" — )• 7 to 
7 — 7- 7r~ have been displayed. In Chapter [5] we have presented our results for a phenomenological 
application of the TDAs [70]. We have examined the expected cross section for vr-vr and vr-p produc- 
tion in exclusive 7*7 scattering in the forward kinematical region, by using realistic models for the 
description of the pion. First, we have confirmed the previous estimates for p production obtained 
in Ref. |131j . even if our results for the cross section are larger of a factor 2. Second, in comparison 
with this previous evaluation of the cross sections, we have improved our findings by considering 
the effect of QCD evolution on the vector current. By so doing, an additional factor 5 has appeared 
in the small ^ region, leading to a cross section for p production of one order of magnitude larger 
than the previous calculation. We have also ameliorated by including the pion pole term in the 
tensorial decomposition of the axial current. Then, an even larger enhancement factor, of about 
60 in this case, has been found in the cross section for pion production. The interference term has 
become larger by a factor 15 with respect to the pure TDA contribution, making the axial TDA 
more accessible experimentally. 

Besides the structure of the pion, we have also been interested in semi-inclusive processes 
involving Parton Distributions of the proton. Since the features of the proton are rather different 
than the pion's one, the description of the proton requires the use of models whose properties are 
adequate for that problem. In particular, we have used two models for the proton which do not 
include the realization of the chiral symmetry but rather confinement: the NR Constituent Quark 
Model of Isgur-Karl as well as the MIT bag model. 

We now draw some conclusions on this part of the thesis, which discusses the Sivers function. 
In Chapter [6] we have presented our results for the calculation of the Sivers function in the two 
models described above. The formalism that we have developed in Ref. [66] for the evaluation of 
the Sivers function is rather general as it can be used in any Constituent Quark Model. The crucial 
ingredient of this calculation has been the Non Relativistic reduction of the leading twist part of 
the One-Gluon Exchange diagram in the final state. It has been shown that the Isgur-Karl model, 
based also on a OGE contribution to the Hamiltonian, is a proper framework for the estimate of 
the Sivers function. The obtained results show a sizable effect, with an opposite sign for the u 
and d flavors. This is in agreement with the pattern found from an analysis of impact parameter 
dependent GPDs in the IK model [93J and allows for the fulfillment of the Burkardt sum rule. More 
precisely, this sum rule is verified with a precision of a few percents. The connection of the Sivers 
function with Impact Parameter Dependent PDs deserves a careful analysis and will be discussed 
elsewhere. In particular, some future directions of work are given in a next Section. 

The Sivers function calculated in the MIT bag model has been revised [72] , and now it fulfills the 
Burkardt sum rule to a level of 5%. Comparing this encouraging outcome with that of the CQM, 
one can say that the Burkardt sum rule is better fulfilled in the CQM. Most probably this has to 
do with the fact that the Burkardt sum rule is associated with transverse momentum conservation. 
The evolved Sivers functions, as obtained in both models, show an impressive agreement with data. 
We would like to remark that the data are described rather well for both flavors. 

In conclusion, we want to stress that we have established, for the first time, that correct model 
calculations provide phenomenological successful interpretations of the Sivers function, which are 
consistent with each other. 

In what follows, we present the possible future directions for both lines of work, trying to merge 
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Spin Physics together with pion structure in the very last Section. 



The Future of Transition Distribution Amplitudes 

In Chapters |4] and [Sj we have introduced the concept of Transition Distribution Amplitudes as 
an extension of Distribution Amplitudes as well as of Generalized Parton Distributions. We have 
presented the results in a model calculation, namely in the NJL model, and compared with different 
model calculations carried out until now. The mesonic TDAs have been applied in the estimation 
of the cross section for meson pair production in two-photon processes, awaiting experimental data. 

A glimpse of distribution parameterizations have been given including phenomenological mod- 
eling. In this short Section, we mention some other approaches and applicability frameworks for 
TDAs. 



Besides the parameterization through the Radyushkin's Double Distribution, there exists the 
so-called Dual Representation [1701 1172] . Duality comes by representing the parton distributions 
as an infinite series of distribution amplitudes in the t-channel, obtained by decomposing the GPDs 
in i-channel partial waves. This representation leads to an alternative parameterization of GPDs 
with respect to (4.66). 



A further study of the TDAs as parameterizing the processes (5.1) through this representation 
should lead to interesting, model-independent results. The latter could, in turn, be compared with 
the model-dependent estimations of the cross sections given in Section 5.4 The inclusion of the 



pion pole contribution renders more appealing the modeling. Moreover, its study is justified, from 
the experimental point of view, by the enhancement of the cross section for vr^vr" production in 
7*7 scattering due to the presence of the pion pole amplitude. 



Among the processes involving the TDAs as the large-distance part, the process 7*7—5' tttt is 
particularly interesting because different kinematical regimes lead to different mechanisms. This 
implies a description of the process through either the pions GDA or the pion-photon TDAs [13j. 
As it is depicted on Fig. [TJ there is a duality in the factorization mechanisms in describing the 
fusion of a real photon with highly-virtual and longitudinally polarized photon. According to the 
kinematical regime, i.e., 

(a) s Q"^ while t is of order Q^; 

(b) t <^ Q"^ while s is of order Q'^; 

(c) s,t<.Q'^; 

(d) s,t~g2. 

the process will factorize following the GDA[^ (regime (a)) or the TDA (regime (6)) mechanism. In 
addition there exists a kinematical domain, namely the regime (c), in which both mechanisms act 
in parallel, giving rise to duality in this region. The model-dependency of the GDA and its impact 
on duality is nevertheless pointed out. The regime (d) corresponds to the large angle scattering d 
la Brodsky-Lepage. 

This interplay between kinematical regimes leads us to think that it would be interesting to 
deepen our understanding of the description of this process through the analysis of the twist-3 pion 



'^The twist-3 GDA here is selected by the "f^jr initial state. 
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Figure 7.1: GDA vs. TDA factorization. From Ref. [S]. 



GDA, e.g. in the NJL model. 



In conclusion, the meson-photon transitions lead to many possible applications that represent 
a natural follow-up of the work displayed in this thesis. 



The Internal Spin Structure of the Pion 

Complementarily to what has been explained in the previous Chapter, it has been shown that, 
even if the target is unpolarized, the momentum distribution of its quarks ejected in SIDIS can 
already exhibit a left-right asymmetry of the intrinsic transverse quark momentum kx relative to 
their own transverse spin st |35j . The Boer-Mulders function is defined as the modulation of 
the probability density due to the correlation between these two variables, via the ^I'^+'Tsl of the 



decomposition (1.28), 



^^'^^'^^^^{x,kT,s) = n{xM)-'^^^^hi^{x,kT) (7.1) 

Like the Sivers function, the function h^'^{x,kT) requires the presence of Final State Interactions 
to allow different orbital angular momentum combinations. 



Obviously, one of the outlooks of this thesis would be to calculate the Boer-Mulders function 
in both the 1 and 3-body models in which we already have calculated the Sivers function. This 
calculation is going to be completed very soon ^73j |^ 

Another future perspective is related to the first part of this thesis, and it concerns the non 
trivial spin structure of the pion. 

Since the pion has spin zero, its longitudinal spin structure in terms of quark and gluon degrees 
of freedom is trivial. However, it would be of interest to investigate the transverse spin structure 
of the pion, which, in the case of a spinless hadron, would come from a modulation of the number 
density due to the interplay between the transverse spin of the quarks and another intrinsic variable. 
The latter is namely either the momentum kx or the Impact parameter b± , which gives the distance 



■^This work has been published after the thesis has been presented. 



113 



between the center of momentum of the hadron in the plane transverse to its motion as shown on 
Fig. 



7.2 



Following the example of the intuitive determination of the sign of the Sivers function from 



the E GPD (1.13), one expects that the sign of the Boer-Mulders function can be related to the 



distribution of transversely polarized quarks in the Impact Parameter space |5T 



Little is known about both the Impact Pa- 
rameter Dependent (IPD) distribution under 
scrutiny and the Boer-Mulders function. This 
holds both from the experimental and theo- 
retical point of views. However, what moti- 
vates our project is the existence of the first 
lattice results for the Form Factors in momen- 
tum space of the concerned pion GPD [33] . 
Our proposal idea is basically to apply the for- 
malism developed in Chapter [2] for the calcu- 
lation of pion distribution in the NJL model to 
the calculation of the pion IPD distribution of 
transversely polarized quarks. In this way, we 
hope to gain some insight about the possible 
ansatz linking the pion GPD calculated in the 
impact parameter space, i.e. correlations be- 
tween the impact parameter and the intrinsic transverse spin of the quarks, to the Boer-Mulders 
function. However, such relations involve a non trivial association of a spatial variable b± and a 
momentum k±. It is henceforth a tantalizing goal to obtain model independent relations. So that 
one should restrain itself to global considerations, such as the relative signs. The latter can be 
afterwards compared with the result in the Lattice. 



> Z 

Figure 7.2: Illustration of the GPDs of the tt^ 
parameter space. 



in the impact 



The IPD distributions have been proven to be related to the GPDs at zero skewness |47j . 
The Fourier transform of the charge density defined with respect to the impact parameter can be 
associated, without ambiguity, with the charge density as a function of the impact parameter in the 
Infinite Momentum Frame [U]. The same can be done for the GPDs, but only for ^ = 0. In other 
words, the GPDs at zero skewness allow a simultaneous measurement of the light-cone momentum 
and transverse position, i.e. impact parameter, distributions of partons in a hadron 



-lA I -b I 



F(x,C = 0-;Aj 



(7.2) 



We now consider the transverse quark polarization, which is the only one allowed since the 
helicity-flip of the state is forbidden for the pion. Quarks with transverse polarization st are 
projected out by the operator ^^[7^ ~ •S"' i'^'^'^lblQ- Their density is therefore 



;6±) 



i.e. a combination of the unpolarized (3.1) and polarized one. The latter IPD parton distribution 
is given by the chiral-odd twist-2 pion GPD, which is allowed by flipping the quarks' helicity. It is 
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defined as 

1 e+^°«A<.p^ 



\z^=Z+=0 J 



2p+ A 



H^ix,C,t) , (7.3) 



where a natural choice for some hadronic scale, say A, would be m^r in the nucleon case, [^his 
scale is rather chosen to be A = 47r /^r in Ref. [53] in order to avoid discussion within the chiral 
limit. However, the chiral limit has been studied in Ref. [83] showing that the choice A = ttItt is 
totally justified. 



The Fourier transform (7.2) of the chiral-odd distribution F^{x,£,,t) leads to a derivative of 
the IPD distribution with respect to 6^; 

4 (h^x; b^) + e*^' s'V — -^ Hf{x; b, 

The dipole term, as expected, leads to a dependence on the direction of b_\_ for fixed st- The sign 
of this contribution seems to be opposite the the one of the Boer-Mulders function (def-bm). In 
other words, the IPD describes a sideways shift in the distribution of transversely polarized 
quarks in an unpolarized proton which has the opposite sign than the Boer-Mulders function. 

In order to give such an estimation in the NJL model, the chiral-odd pion GPD defined in 



Eq. (7.3) has to be calculated using the formalism described in Section 2.1 In particular, the result 
reads 

Hf{x,i,t) = i8m^N,gl^J^P^{x,^,t) . (7.4) 

After that, the procedure is as described above: we first Fourier transform to the Impact Parameter 
space, and then match the sign conventions with the one used by, e.g., the Lattice Collaborations 
to finally yield our prediction for the sign of the Boer-Mulders function for SIDIS off pions. 

A study of the Boer-Mulders function in the quark-spectator-antiquark NJL-like model is pro- 
vided in Ref. [143J. The sign for the Boer-Mulders function is found to be negative. On the other 
hand, a naive approach of IPD pion distribution function has been given in Ref. [53], where a 
positive sign for the distribution function was found. We thus expect a negative hj-^, a result that 
our calculation should either confirm or invalidate [71j. 

Lattice calculations nevertheless yield the same sign. For example, the latest lattice calculations [43j 
give access to x-moments of the quark spin densities p'^{b'^, s±). They are obtained from the gen- 



eralized Form Factors in momentum space by applying the Fourier transform (7.2). In turn, the 



generalized Form Factors are obtained through the usual recipe of taking the Mellin moments of 



the GPDs. The results are shown on Fig. 7.3 



Note that, by analogy to the nucleon case, the standard nomenclature is iJJ or instead of H^. 




A Final Word. . . 

We conclude this extended outlook by saying that the possibilities of learning something more 
about the partonic structure of pions and protons seem to be infinite! Hopefully, we will be able to 
unveil many secrets of the hadron's spin structure by performing simple and intuitive calculations, 
which can be implemented by phenomenological insights. 




Notations and Conventions 



In this work, natural units are adopted, 

h = c = l . (A.l) 

The metric tensor is 

9^" = = diag(+l, -1,-1,-1) . (A.2) 



A.l Light- Cone Vectors 

Here we first define the vector components in terms of light-cone variables and then, the Sudakov 

parametrization [28j. 

A four-vector is generically written as 

,0 „1 2 „3 



(o^a^a^a^) . (A.3) 



We define the light-cone components 



a± = ^{a'±a') , (A.4) 

hence has the form 

fl/^ = (a+,a-^,a") , (A.5) 
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with a-*- = {a^,a^). 

The scalar product of two vectors is 

a.b = a+6~ +a~6+-a-^.6-^ , (A.6) 

while the norm is 

= 2a+a--a^^ . (A.7) 

We introduce the light-Hke vectors, i.e. vectors which norm is zero, 

= ^(1,0,0,1) , 

= p+ (1,0^,0) , 

= 0,0,-1) , 

= ^(0,0^,1) , (A.8) 
and those vectors have the following properties 

p.n = 1 , 

■n+ =p- = . (A.9) 



It is convenient to introduce the light-cone decomposition (also called Sudakov decomposition) 
for a vector 

= (a.n)p^ + a^'* + (a.p)n'' , (A.IO) 

with a-L/^ = {0,a^^',0). 

In a "light-cone basis" the metric tensor g^'^ has non-zero components = g ^ = 1 and 
g^^ = g^"^ = — 1. Moreover we can construct the perpendicular metric tensor 

^M^± = gi^^ -(p^n'' ^fn^) . (A.ll) 



A. 2 Pauli Matrices 

The Pauli matrices will often be used in the form r, where 



(A.12) 



A. 2 Pauli Matrices 
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with 

r+ = ^{r^+ ir'^) and = ^(r^ - ir^) , (A.13) 

where r* = cr*, i.e. the usual Pauh matrices. In the text, the notation r refers to 

r = (r-^r-^r--) , (A.14) 

with 

r'^^ = -^(r^ + zr^) and r" = ^jCr^ - ir^) . (A.15) 



Hence, we have the foUowing results 



on - _ / 

0/ ' ~ ll 



(A.16) 

-~-" = (o I) ' -"-' = (o o) ■ (^-1^) 



0\ + _ /I 



Therefore, the following relation can be written 



(A.18) 
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Discrete Symmetries 



In this Chapter, we treat the discrete symmetries that link the distributions for a same kind of 
particle among themselves. General relations are given and then applied to the particular cases of 
pion DAs and TDAs. 

For a generic particle V with momentum p, spin projection Sz and charge Qj, we write the 
particle state \V{p^ Sz,Qi)) ■ Parity, time reversal invariance and charge conjugation act on this 
generic state in the following way 

P\V{p,Sz,Qi)) = Vp\'P{-p,s,,Qi)) , 

T\rip,Sz,Q^)) = VT{-ir-'^\V{-p,-Sz,Qi)) , 

C\V{p,s,,Q,)) = vc\r{p,s,,-Qi)) , (B.l) 
with rjp, rjT, rjc the respective phases. For the tt"*", which is J^^ = ^, are rjp = rjT = — l, ??c = 1- 



In the light-front, the T symmetry as defined by the second expression of Eq. (B.l ) changes x"*" 
to —x~ and vice versa. It is therefore useful to define the V symmetry has the combination P^T, 
what gives 

v\r{p,sz,Qi)) = wi-ir~''\r{p,-sz,Qi)) , (B.2) 

with r]v = tIpVt- 

The distribution functions for a generic transition between the state V{p, Sz,Qi) to the state 
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V'{j? , s'^, Qi) are associated to 



(B.3) 



where V and can be different particles as well as the vacuum. 

Let us now apply the definitions of the discrete symmetries to the generic distri- 



bution (B.3) : 



• Time reversal , as given in its hght-cone version by Eq. (B.2), transforms our generic distri- 
bution in the following way 

with related to the nature of the current, i.e. Cv = 1 for T = 1, ^75,7^^ while Cv = — 1 for 
r = 7/,75,o-^!/- 



Similarly, applying charge conjugation, Eq. (B.l), we find 



with = 1 for r = 1, Z75, 7^75 while Cc = -1 for E = 7^, cj^i,. 

By applying both the above results, we obtain, for the CPT symmetry 

gq_^^'{x,p,Sz,0,s'^) = CcPTGg'^^{-X,p,-s'^,p,-S^) , 

where (cpt = CcCv- 



Let us finally give the hermitian of the generic transition (B.3) 



We can now go to particular distribution functions: 



Pion Distribution Amplitude 

The pion DA would be 



e::7°((^-l/2),P,0) = /^e^(--^)^"-(O|,-(-07+75r-g(0|vr+(p)) 



z+=z^=0 

iV2Ur^{x) . (B.4) 



The hermitian conjugate of Q'^ ^^((x — 0) is given by 



Using the relations Eq. (B.l), the time reversal invariance leads to 
charge conjugation to 

r~{x) = r\l-x) ■ 

and CPT to 



Vector Pion Transition Distribution Amplitudes 

The ir^-j and 7-7r~ vector TDAs are defined 



ixP+z 

2tt 



2J ' " V2 



—iee f+^P"" P (v -n) ^'^ ^^(^'^'^) 



27r 



2J ' ^ V2 



The skewness variable is defined as 

i = {Pn - P^)-^ /2P-^ 



The hermitian conjugate of the vr^ — )• 7 vector TDA is given by 

[Su^Ti^^p^p ^^)) = Sj^f(x,p,p,e) ; 
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Time reversal invariance leads to 

iee^e'^^'^'^n, p /^T-P^)^ V-^^^x,C,t) 



charge conjugation to 

gll^^'^ix,p,p ,e) = C,cVcGl_^''{-x,p,p ,e) , 



v2/7r 



v2/7r 



V-^^'^{x,C,t) = V---^^{-x,^,t) 

and CPT to 

Ql'_^^{x,p,p ,e) = CcPTVCVvGj^^ {-x,^,p,e) 
iee^e'^-P'^n^ P,^^^^^^V-^^^{x,^,t) 

v2/7r 



V-^^^{x,^,t) = V^^--{-x,-^,t) 



Axial Pion Transition Distribution Amplitudes 

The 7r+-7 and 7-7r~ axial TDAs are defined So that 

A-^^Ux,tt) 



e(8^.(pt-pt)) 



— t 

A'y-*--{x,-^,t) 

V2U 
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with e(^) = 1 for ^ > and e(0 = -1 for ^ < 0. 
The skewness variable is defined as 



Combining the results for Q^^'^ {x^p^p ,£) and for Q^^^^ (^'f^''^) expressions 
for the axial TDAs. 

The hermitian conjugate of the 7r+ ^ 7 axial TDA is given by 



I e [e -{p^ -pjj + e {e ■ {p^ - p^)) ^^^sTT^ e(0 

= e [e^-{p^-p^)j - 



x + C 
2^ 



V2U 



+ e{e-{p^- p^)) ^^^2^ e(-0 (^^^ 



Time reversal invariance leads to 



charge conjugation to 



and CPT to 



Those relations remain unchanged under evolution. 
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The Nambu - Jona-Lasinio Model 



A pre-QCD model for the pion was introduced in 1961 by Y. Nambu and G. Jona-Lasinio with 
nucleon degrees of freedom [1511 1152) . This study was motivated by the theory of superconductivity 
of Bardeen, Cooper and Schrieffer (BCS) and so built in analogy with it. The basic idea of the 
theory is that the mass of the "Dirac particles" , in the case of the original paper of Nambu and 
Jona-Lasinio a nucleon, could be due to some interaction between massless fermions that have the 
same quantum numbers excepted chirality. In the same fashion as in the theory of superconductiv- 
ity, the gap is caused by nucleon-nucleon interaction. Therefore, starting from a massless nucleon, 
i.e. an eigenstate of chirality, we arrive to a massive nucleon. This is the "solid state physics" 
formulation of the chiral symmetry breaking. 

The lagrangian was contructed to conserve chirality and "nucleon number" . In quark language, 
more precisely in a two-flavor quark model, it has to be understood as SUv{2) x SUa(2) x [/^/(l) in 
the second version of the model which exclude in accordance with experiment. Transforma- 

tions under the isospin SUvC^) symmetry and the chiral S'J7^(2) symmetry are respectively given 



by 



■0 — )• e 




(C.l) 



with CO and 6 the respective generators. 

The NJL lagrangian density in its 2-flavor version is 



^NJL 



ill) il) -\- G (i^ip)'^ + (i/A^^Tip) 



2 



(C.2) 



with 
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C.l The Chiral Symmetry Breaking 

Nowadays it makes sense to speak about quark degrees of freedom instead of nucleons. The quark 
mass topic is quite puzzhng: quark models give us high values while chiral symmetry seems to 
impose small quark mass. This can be understood by considering the average interaction of a 
quark with the other quarks in the mean-field approximation. The concept of constituent quark 
mass can then be introduced in opposition to current quark mass. The quark self-energy in a mean- 



field approximation is given by the Hartree-Fock (HF) approximation, Fig. C.l This consists in 
considering perturbations around the free " 

4 = Co + C^ , (C.3) 

where the self-energy lagrangian Cy, = — S ^^/j is chosen in such a way that the interaction lagrangian 
does not contain any self-energy effect. 




Figure C.l: The Hartree and Fock, respectively direct and exchange, contributions to the quark self-energy. 

The total self-energy S is the sum of the scalar and pseudoscalar ones and is given by, in the 
HF approximation, 

S = T^' + TP' , 

= 2G |Tr ^^(x, x) — j;) + (i75r)Tr j;)(i75T) — (i75T)«5'(x, x)(i75T)| , (C.4) 

where S'(x, x) is the quark propagator and satisfies 

(i^^ -mo-S)5(x,x') = 6^{x-x') . (C.5) 
The current quark mass mo comes from Cq. 

Given the "observed" constituent quark mass m, the propagator S{x, x) should also satisfy the 
following equation of motion 

{i^^-m)S{x,x') = 6^{x-x') , (C.6) 
so that m — mo = S. The solution of Eq. ( C.5| ) in momentum space is then. 



Sip) = -i 2 • C.7 



If we insert this result into the expression for the self-energy Eq. (C.4), performing the traces and 
noting that the integral odd in p does not contribute, we find 

E = ™-™„ = 2«(iV.A',+ l) / Ifj^ . (C.8) 



C.l The Chiral Symmetry Breaking 
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At this point, it is a reasonable assumption to set the smah current quark mass to mo = 0. Even 
if the exact chiral symmetry forbids mass terms in the Lagrangian, there exists a solution of the 
type m 7^ when the current quark mass is zero. Actually the equation Eq. (C.8) admits two 
solutions. The trivial one is m = and is in agreement with symmetry requirements. The second 
and nontrivial one is 



1 



2iG NrN 



1 



(27r)4 p2 



(C.9) 



This equation is called the self-consistency or gap equation in analogy with the nontrivial solution 
giving rise to the gap in the BCS theory. The non-trivial solution, as solution of the gap equation 



Eq. (C.9), exists if the value of the coupling strength G exceeds a critical value |151l I152j . 



Regarding the coefficient in front of the integral in (C.8), we observe that the first term, prO' 



portional to G Nc, comes from the direct diagram of Fig. C.l| (Hartree contribution), and that the 
second one, proportional to G, is the contribution of the exchange diagram (Fock contribution). 
We can understand the NJL model as an effective theory connected to QCD in the large Nc limit. 
In that case, G Nc ^ C'(l) and G ~ Nc^- Therefore, since we can consider that we work in a 
A''~^-expansion, it is consistent to ignore the Fock contribution. 

It would make sense to work out which of the two solutions is the ground-state solution. By 
analysing the vacua of both the trivial and non-trivial solutions, we can related them in a BCS 
variational vacuum fashion. The vacuum for the non-trivial solution |0)m is expressed as a su- 
perposition of pairs of zero total momentum and zero total helicity over the trivial vacuum |0)o 
|15H 1152] . We consider the Dirac equation in both cases, and with similar initial conditions. 



and decompose the respective fields into their Fourier's components 

r{x) = ^{6"(p,5)tx-(p,s)e*?'- + dt-(^?,s)^;"(p,5)e-^P-} 



(C.IO) 



(C.ll) 



for a = or m. The spinors for particles and antiparticles are respectively u"{p,s) and v°'{p,s). 
The annihilation operators for particles and antiparticles are respectively b°'{p, s) and d°'{p, s). 
The annihilation operators for the fields and are related to each other according to they 
have the same initial conditions. Therefore the vacua defined with respect to both fields, i.e. 



b°{p,s)\0)o = d\p,s)\0)o - 

b"'{p,s)\OU = d^{p,s)\0)m -- 

can also be related one to each other. It can be found |125) 








|0)r 



n 

p,s 



COS 



(f){p) 



+ s sm 



0(p) 



6tO(p,s)dtO(-p,s)}|0)o 



(C.12) 



(C.13) 



This is a BCS-like variational vacuum: a superposition of pairs of zero total momentum and zero 
total helicity. This is also called the Bogoliubov-Valatin approach because the operators that annihi- 
late such a vacuum are found through a Bogoliubov-Valatin transformation |125j . The angle (p{p) is 
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determined by minimizing the vacuum energy VF[(/)(p)]rre =m {0\H\0)m, obtaining tan [(j){p)] = m/p. 



The energy difference of the two vacua Wm — Wo is found to be negative |15H 1152] . So \0)m is 
the true vacuum and the nontrivial solution is the solution that minimizes the energy. In particular, 
minimizing the vacuum expectation value in the nontrivial solution case leads to the gap equation 



Eq. (C.9). 



The true vacuum \0)m. is not invariant under SUaC^)- The consequence is a (dynamical) breaking 



of the chiral symmetry SUaC^). In other words the interaction given by Eq. (C.2) between two 



massless quark with the same quantum numbers excepted chirality gives rise to the gap equation 



Eq. (C.9) and generates the constituent quark mass m. 



C.2 Regular izat ion 



The gap equation Eq. (C.9) 



1 



2G NcNf h 



(C.14) 



has not been resolved yet. In fact the integral /i, given by Eq. (E.l), is divergent and we need 



to define a regularization scheme. The NJL interaction is a contact interaction that renders the 
theory nonrenormalizable. Therefore the parameters introduced by the regularization scheme are 
integrally part of the model in itself. The NJL model plus its regularization scheme is regarded as 
an effective theory of QCD. 



Several regularization schemes have been studied. We mention the noncovariant 3-momentum 
cutoff scheme, the covariant 4-momentum cutoff scheme or the proper-time regularization |125j . 

In the study of PDFs, GPDs and TDAs, it is important to preserve Lorentz covariance and 
gauge invariance. The Pauli-Villars scheme fulfills these conditions. The idea of this scheme is to 
introduce the minimum number of regulating masses rrij and constants Cj in such a way that the 
result becomes finite. The following replacement is made 

with m| = m^+jA^. We use the standard values given in |112j for the Cj coefficients, i.e. cq = C2 = 1 
and ci = —2. A and m are the regularization parameters. They are determined by calculating 
the pion decay constant fn and the quark condensate |125j . For fj^ = 93MeV and {uu) = (dd) = 
-(250MeV)3, we obtain 

m = 241 MeV , 

A = 859 MeV . (C.16) 



C.3 The Pion and the Sigma 



In its two-flavor version, the NJL model Eq. (C.2) gives rise to 4 mesonic modes. 



C.3 The Pion and the Sigma 
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The scalar interaction is responsible for a isoscalar (/ = 0) scalar J^^ = 0^^ mode. The 
scalar mode is associated with the o"- meson, also known as /o(600) [11]. The O"*"*" sector is both 
theoretically and experimentally complex. The scalar modes are difficult to solve because of their 
large decay widths. 

The pseudoscalar channel is responsible for the isovector (/ = 1) pseudoscalar J^^ = ^ 
modes identified as the pions which are then the three Goldstone bosons coming from the chiral 
symmetry breaking. 

There are two ways of expressing the pion. In the original paper |15H 1152] , Nambu and Jona- 
Lasinio consider the pion as a bound-state of two nucleons, what is nowadays reinterpreted as a 
bound-state of two quarks, and study it through the use of the Bethe-Salpeter (BS) equation in 
the ladder approximation. Another approach is to consider the effective interaction resulting from 
the exchange of a pion in the Random Phase Approximation (RPA). The two methods are identical. 

Even though both approaches will be used in the next chapters, according to we are interested 
in pions as bound-states or in pion exchange, we choose to introduce here the pion as a bound-state 
through the Bethe-Salpeter equation. 



The Bethe-Salpeter (BS) amplitude for the bound-state vr with four-momentum P is given by 

Xp{xi,X2) = (0|Tg(xi)g(x2)|vf(P)) . (C.17) 

The NJL model contains a four-fermion interaction. In the ladder approximation, only direct 
diagrams are considered. In our case, this approximation is equivalent to considering the iteration 



of the simplest closed loop Fig. C.2 with the kernel V{p,p']P) = 2iG{i'y^T) . The equation of 



Bethe-Salpeter is easy to handle due to the point-like interaction of the NJL model. 




Figure C.2: Graphs corresponding to the BS equation in the ladder approximation. The thick line represents 
the bound-state. 



The integral equation generated by the chain of diagrams Fig. C.2 is the Bethe-Salpeter equation 

m 

iS-Hp) Xp{p)iS-\p-P) + 1 {V{p,p";P)xp{p")) = , (C.18) 

here in momentum space. S~^{k) =^ —m is the inverse of the Feynman quark propagator. 
Solving the Bethe-Salpeter equation in the ladder approximation 

^S-\p)xpiP)^S-Hp-P) = Gj,T I ^2Tt {ij,T.xp{p")) , (C.19) 
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we find the expression for the Bethe-Salpeter amphtude 

XP{P) = -gnqqiS{p)-i^TiS{p- P) 



(C.20) 



Inserting the Bethe-Salpeter amplitude Eq. (C.20) back into the BS equation Eq. (C.18) gives 
a self-consistency condition 



1 



2G{-i) 



Tr 75T^"S(p)75r" Sip-P) 



(C.21) 



The integral represents the pseudoscalar proper polarization Hps shown in Fig. C.3 The mass of 



the pseudoscalar mode is obtained for P^ = m^. That is, the mass of the pion is found by solving 
the following equation 



l-2GIlps(ml) = Q 



(C.22) 



We have completly defined the pion. By resolving Eq. (C.22) at k"^ = m^, the mass m.,^ is found. 
Now we would like to show that the pion described in such a way is the researched Goldstone boson, 
i.e. that it has zero mass in absence of the current mass quarks. In order to solve such a problem. 



we allow a small current quark mass in the Lagrangian Eq. (C.2) of the type —mQipip and solve 



Eq. (C.22). 




Figure C.3: The pseudoscalar proper polarization Ilps{k^). 



The proper polarization Fig. C.3 is explicitly given by 



d'^p 



-i ANcNf 



Tr|i7V''" 



p + + m 



{p+ - m2 



+ m 



{p - 1)2 - m? 



- p2 _|_ ^ 



d^p 

(27r)4 §)2 _m2][(p- 1)2 -m2] 



(C.23) 



where the change of variables p = p + /c/2, with P = k, has been performed. We bring Eq. (C.23) 
into elementary integrals so that we can easily extract the pion mass myr. The numerator can be 



reexpressed like (m — p 



-) + -2 and the denominator 



[{p + 1)2 - m2][(p - 1)2 - m2] 2[p2 + ^ - m2] I [{p + |^2 _ ^2" 



+ 



1)2 -m2] 



(C.24) 
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So the pseudoscalar polarization Eq. (C.23) takes the forn|^ 

Upsik^) = 2NcNf[2h-k^l2{k^)] 



(C.25) 



with Ii and l2{k'^) respectively one and two-propagator integrals, see Appendix |Ej Moreover, the 
self-consistency condition Eq. (C.21) gives us the gap equation Eq. (C.14) when the mass of the 
bound-state is zero, i.e. 



The first term on the right-hand side of Eq. (C.25) can be replaced by its expression in the gap 
equation Eq. (C.9). Adding a small current quark mass term to the gap equation, we find for Ii 



SGNcNfh 
In doing so we obtain, using Eq. ( |C.22 ), 



1 



mo 
m 



1 - 2GUps(.k^) 



mo 



m 







For fc^ = m^, we may deduce the mass of the pion 

2 mo 1 



m^ 



m AGNcNfhiml) ' 



(C.26) 

(C.27) 
(C.28) 

(C.29) 



In the chiral limit, i.e. in the absence of current quark mass mo, the mass of the pseudoscalar mode 
vanishes, as expected. 

The quark-pion coupling constant g-j^qg has been determined by the standard normalization of 
the Bethe-Salpeter equation jll2j 



2 

Sirqq 



12(/2(m2) + m2 (^) |p.^^. 



(C.30) 



with l2{p) given in Eq. (E.2). 



Alternatively, we can determine the interaction properties of the bound-states in the NJL model 
through the study of gg-scattering amplitudes |15H 1152] . In the ladder approximation the bound- 
states are expected to be real stable particles. So should be the poles corresponding to the virtual 
exchange of, in the pseudoscalar channel example, pions. The bound-state is therefore the iteration 
of the bubble diagrams Ilps{k'^) in the pseudoscalar channel Eq. (C.3). 



The scattering matrix Fig. C.4 generated by the exchange of a pion is then 



iUijik^ 



{i-f5T') 2iG + 2iG{-iUps{k^))2iG + . . . {i-f^r^ 



2iG 



1 - 2GUpsik^) 



(C.31) 



shift has been performed in order to obtain the 2/i term. This still holds if one uses the Pauli-Villars 
regularization scheme. What is not obvious in other schemes. 
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+ . . . 



Figure C.4: Scattering matrix generated by the bubble diagram Fig. |C.3 



vice versa 



The r* selects the appropriate channel. For exchanging a Tr"*", r* = t'^ and = i 
for a TT~ and r* = r-' = t'^ for a vr''. This scattering matrix can be compared to the one coming 
from the minimal local interaction Lagrangian that describes the coupling of a pion field to quark 
fields. Namely, the effective pion exchange is 



(i75r*) 



mt 



(C.32) 



Hence, in the light of Eq. (C.32), we can immediately determines the mass of the pion by finding 
the pole of Eq. (C.31) at = rn^ . In other words, the self-consistency equation Eq. (C.21) of the 



Bethe-Salpeter equation comes to be the pole of Eq. ( |C.31 ) 



1 - IGIipsik^ = m^) = 



(C.33) 



In the same fashion, we can relate the quark-pion coupling constant to the residue at the pole 



Eq. (C.33) 



2 

9-Kqq 



With lipsik"^) given by Eq. 



dk^ 
(|C.23l), the 



expression Eq. 



(C.34) 

(C.30) is recovered. 




Figure C.5: 



The mass of the scalar mode a associated with the term 
{'iptp)'^ can be obtained in the same fashion. The zero of the 
function l—2GIls{k'^), with Ils{k'^) the scalar bubble diagram 



Fig. C.5, defines its mass nia- Also its coupling strength can 



be related to the residue of the pole at k 



iUs{k^) 



The scalar loop is given as 



(2vr)4 
-ANcNf 



m 



} . 



(27r)4 [(p + 1)2 - m2] [{p -l)2-rn^ 



(C.35) 



Using the same trick as in the pseudoscalar case, i.e. writing Eq. (C.35) in terms of elementary 
integrals, we find 



n,(A;2) = 2NcNf [2Ii - {k^ - Am^) h{k^ 



(C.36) 



C.4 Partial Conservation of the Axial Current 
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As before the one-propagator integral is extracted from the gap equation and we find 

l-2GUs{k^) = — +4GN,Nf{k^ -Am^)l2{k'^) , (C.37) 
m 

which is exactly zero. On the pole k'^ = rn^, we find 

2 ^ I , . 2 

m AGNcNfhiml) ^ ' 

~ 4m2 + m2 . (C.38) 

This result linking the pion and the sigma masses is independent of the regularization scheme 
employed. In the chiral limit, the mass of the c-meson is non-zero and its value is 

ml = 4m^ . (C.39) 



We have described a model for the pion. In the NJL model the pion arises through the dynami- 
cal breaking of the chiral symmetry. As a consequence its mass should be zero in the chiral limit. 

The NJL Lagrangian contains a point-like interaction that renders the theory nonrenormaliz- 
able. The model should always be completed by a regularization scheme fulfilling the conditions 
imposed by the theory: the existence of a "gap" and, therefore, zero-modes or Goldstone bosons. 
Another shortcoming is that the model is not confining. 

The Nambu - Jona-Lasinio model together with its regularization scheme is treated like an effective 
theory and should therefore be used in a characteristic scale imposed by the same regularization 
scheme. 



C.4 Partial Conservation of the Axial Current 



We showed in Section C.l that the true vacuum is not invariant under SU One might wonder 
what happens for the axial transformations. Familiar conservation laws are derived from the 
conservation of the currents associated with manifest symmetries of the, e.g., strong interactions. So 
the vector SC/y(2) x \Jy{\) transformations, which are manifest symmetries, imply baryon number 
and isospin conservation. 

The 5C/a(2) symmetry being broken, we restore the chiral symmetry by assuming that the axial 
currents 

Al{x) = g(x)7^75yg(x) , (C.40) 

are partially / approximately conserved. The symmetry is implemented in the Goldstone mode 
with the pion as massless particles. Besides the fact that it is spontaneously broken, the chiral 
symmetry is explicitly broken by the mass term in the Lagrangian. That is, when the current 
quark mass are non zero, the pion which is effectively massless, acquires a mass, just as explained 
in Section fC. 31 

Consider that the axial current matrix element between a pion state and the vacuum can be 
parameterized by 

{0\AUx)\7T'ip)) = zp^P'^Ue-^P-^ , (C.41) 
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with the pion decay constant. It is experimentally estimated to f-^ = 93 MeV pTT] through the 
decay 7r~ — >• //"z^ 



r 



G"^ mlfl {ml - mlf 

cos (7(7 , 



47r m| 

with 6c the Cabibbo angle and G the Fermi constant. 
The conservation of the axial current implies that 

iS)\d^^A^^{x)\7:\p)) = ml5=^Ue-'^-^ = Q . 

It means that m^ f^^ = 0, so that, is chosen to be zero since the Goldstone realization of the 
symmetry is in agreement with nearly massless pions. 





Figure C.6: Axial current in the nucleon. On the right panel, the pion pole contribution. 



However, when restoring the quark (current) mass term, the axial current is no longer exactly 



conserved. The pions acquire a mass as given by Eq. (C.29). Hence, it is worth finding an expression 
for the current conservation when extrapolating from pi^ = to = m^, with 7^ 0. This can 
be done by identifying the divergence of the current with a smooth interpolating pion field tt^ , 



ml vr^ (x) 



(C.42) 



This behavior of this interpolating field can be guessed from the analysis of the matrix element 
of the axial current between nucleon state^ Fig. C.6, Current conservation in the exact chiral limit 
effectively implies the exchange of a massless particle, identified here, to be the pion. It therefore 
means that the matrix elements of the divergence of the axial current are dominated by a pion pole 
(on the right panel of Fig. C.6) for small momentum transfers \q'^\ < ml 



{Aip + qWAi\B{q)) 



Res 



mt 



+ ... 



(C.43) 



By consistency with the model calculation, it 
is important to evaluate the pion decay with the 
technology developed in the previous Sections of 
this Appendix. In the NJL model, the expression 
Eq. (C.41) becomes, see Fig. C.7 



7m 75 



^See, for instance, classical notebooks Refs. [1611 1112] . 




Figure C.7: 
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2 

(C.44) 



which gives us the expression for in terms of 2-propagator integrals 

U = -2NcNfg^qqml2im,p) . (C.45) 
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Transition Form Factors 



In the NJL model, the vector form factor is 



mm^ h{p,P ) 



(D.l) 



3 V2 

In order to calculate the form factors we need the expression for the three-propagator integral. In 



the particular case where = and p'"^ = 0, the expression for I^{p.,p') given by Eq. (E.4). 



The axial form factor involves the three-propagator integral as well, but also the two-propagator 



integral given by Eq. E.2 



V2 ihip,p) + 



mi — t 



(D.2) 



The TT^ —7- 7*7 form factor can be obtained from the vector form factor through a isospin 
rotation 



7r7 7 



it) 



(D.3) 



In the chiral limit and for t = all the considered form factors have the simple expression: 



(0) 



V2 



Ff (0) 



n (0) 



V2 4^V- 



2m^ 



+ 



m 



m 



0.192 X (1 - 0.108) GeV~^ = 0.171 GeV"^ 



(D.4) 
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Transition Form Factors 



The first coefficient of the right hand side is what is expected from the axial anomaly contribution 
to TT ^ 7*7 decay. The term between brackets has a small correction to the expected value of 1 
due to the finiteness of the regularization masses. In the NJL model not only the quarks, but also 
the counter-terms run in the triangle diagram of the axial anomaly. In a proper renormalizable 
theory this correction disappears in the limit A oo. 




Elementary Integrals 



Regularization of the following integrals is insured by the Pauli-Villars regularization scheme, i.e. 
see Section IC. 21 

E.l One, two and three-propagator Integrals 

The one-propagator integral Ii is given by 

J _ . [ d'^k 1 

1 2 



The two-propagator integral /2(p) is given by, in a general case, 

d^k 1 



^c,m|ln-| . (E.l) 



(27r)4 - m?- ^ ie){{p + ky - m? + ie) 

1 



In particular, for < p < Am , Eq. (E.2) becomes 



'^(f') = (S)5E'=.{>";;^+2^^;^"=«S^^} . (E,3) 
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Elementary Integrals 



The three-propagator integral l3{pi,P2) is given by 
h{Pi,P2) 



(27r)4 - m2 + i€){{pi + fc)2 - m2 + ie)((p2 + A;)^ - + ie) 

= 7A2 r^^^ln^^ , (E.4) 

with, in the particular case where Pi=P2 = fn^-, 

D = t^z'^ -4.zH{m] + mlz{z-l)) ; 

E = zH- 2m^j - 2mlz{z - 1) ; (E.5) 

and with, in the particular case where pf = and p| = 0, we obtain 

D = t^z'^ + z^{l-zfml + 2tz^{mlz{l- z)-2m]) , (E.6) 

E = zH + z{l - z)ml - 2m] . (E.7) 

E.2 Light-front Integrals 

E.2.1 Two- propagator light-front integrals 

We can now calculate the three 2-propagator integrals, l2,A{x,^,t), l2,p' and l2,p{x,$,) defined as 

d'^k S {(x - l)p+ + k+) 



+ 



l2p, = ip' 



(27r)i ((P - A;)2 - m2 + ie){{P' - kf - + ie) 

d^k 5 {{x - l)p+ + k+) 

(27r)4 (A;2 - w? + ie){{P> - k^ - m? + ie) 



~ , ^, , [ d'^k 6 ((.X - l)p+ + ,^ ^, 

/2,P(X,0 = ^P J (fe2-m2 + .e)((P-fc)2-m2 + ze) " ^^"^^ 

The first step of this integration consists in decomposing the momenta in their hght-cone compo- 
nents. The delta function selects the +-momenta carried by the quarks so that the integral over 
dk'^ is straightforward. The integral over dk~ is performed with the help of the Residue Theorem 
so that we are left with an integration over dk^. The latter is straightforward. 

GPDs' kinematics 

The result is 

1 1 2 ^2 I t x^-e 

4A(x,C,i) = — ^- ^ c,- In ' \ e{x + O0{^ - x) . (E.9) 



E.2 Light-front Integrals 
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The same steps can be done with the /2,P' and l2,p{x,(,) integrals: 



l2,P' 



1 1 ^A , mj + m^ ^ ^^I'^^y 



(47r)2 1 - ^ 



EC,- In- 



i=o 



(x-i)B(\-x) 



(E.IO) 



The l2,p{x,^,) integral is defined on — ^ < x < 1. 



m, + m: 



j -r ' "-TT (l+g)2 



0(x + 0^(1 - a:) 



(E.ll) 



TDAs' kinematics 



In the TDA's kinematics, the Eqs. (E.9, E.ll) do not change but l2,p' becomes 



11 TTl ■ 



(E.12) 



i=o 



E.2. 2 Three-propagator light-front integrals 

The three-propagator light-front integral is given by 

GPDs' kinematics 



with 



ip 



d^k 5 {{x - l)p+ + k+) 



(2^)4 - m2 + ie){{P - kf - m2 + ie)((P' - kf - tv? + ie) 

2 



2 (47r)2 ^ ''^^ 



r 1 -m2x(l-x)+m2(l-^2)_ t(^_^)2^(;L_^)y^ 

|e(l-x)0(x-e)^ln ,v w 2V V V ; 



L> -mlx{l -x) + m|(l - _ i (1 _ ^)2 _ (1 _ x)VD 



+e{C-x)9{x + C)—ln 



1 , -ml^{^ + x) + 2m]^{l + ^) - ^{1-x){x + + {x + 0VD 



D -ml^{^ + x) + 2m]i{\ + ^ - |(1 - x)(x + - (x + i)^^ 

(E.13) 



2 



D = ml(^ + -(1 - x)^ - Am^jCmi + tm^x(l - x) - m^t(l - . (E.14) 



4 
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TDAs' kinematics 

In the TDA's kinematics, the skewness variable is bounded by t/(2rri^ — t) < ^ < 1. We therefore 
have to consider both the cases with positive and negative values of ^. 

For positive values of ^, the three-propagator light-front integral becomes 



d^k (){(,)■- l)p+ + k+) 



(27r)4 (fe2 - ^2 + ie){{P - kf - + ie){{P' - kf - + ie) 



1 ^ 



(47r) 



leix + \mm -x)^\n — — - — ^ — ^ — — 

'^'^ '^^D ml{x-0 + t{x-l) + Ae^^m]-VD 

+e{x-\i\)e{i-x)^hi^ — — ^ '- — f^^^ — —\ , (E.15) 

ml{x-i)-t{x-l) + 2^-^m] + ^^ 



where we have introduced the Pauli-Villars regularization parameters and with 

D = {ml{x-i) + t{l-x)f + A m) {2mli{l - - t{l - e)) ■ (E.16) 



On the other hand, for negative values of the skewness variable, the three-propagator light-front 
integral becomes 



d^k 5 {{x - l)p+ + k+) 



ip 



(27r)4 - m2 + ie){{P - kf - m? + ie){{P' - kf - m? + ie) 



1 ^ 

^ ' j=o 



^(^+ICI)^(IC|-x)^ln 



1 ^ -ml{x^i)^t{x-\)-Ai\S(m]^^ 



^ -ml{x + i)^t{x-\)-Ai^^m]-^fD 

1 ml{x-^)-t(x-\) + 2^m\-^^ 

+d{x-\i\)e{\-x)^\T,— ^ —\ . (E.17) 

^ '^'^ ^ ml{x-^)-t(x-V) + 2'-^m] + ^^ ^ ' 



Double Distributions and the a-representation 



This Appendix is based on the work by Broniowski, Ruiz Arriola and Golec-Bernat [45] . 

In this Appendix we aim to give an exphcit calculation of the Double Distribution of the form 
given by reduction formula (3.25 ). This is conveniently done through the Feynman a-representation 
of the scalar propagators. 
So the scalar propagator reads 

1 f°° 

Sk = ^ = / dae-"(^'-"'+^^) . 

— + le Jo 

The 3-propagator Integral 

We write the 3-propagator integral in the asymmetric notation for it is more appropriate for our 
purposes, omitting to write the +ie. It reads 



d^k 6{X -k- n) 



(2^)4 (A;2 - m2) {{k + A)2 - m^) {{k - p^Y - m^) 



with A = — Ptt. In terms of exponential representation of the propagator, the vector TDA 
becomes 

da\ d/3 / dTe-^C^'— ')-^((^'+^)'— ')-^((^-^'-)'--') . 
Jo 

We perform the following change of variables 

s = Q + /3 + 7 , y = - , z = - , 

s s 
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Double Distributions and the a-representation 



with a, 6, c > so that we can constrain y, z to have support between [0, 1] and also y + z < 1. The 
constrains on the limits for y, z are automatically imposed by the a-representation of the scalar 
propagator. 

Concentrating on the argument of the exponential, we find 



(a + /3 + 7)(/c^ - m^) - /3A - jp^ - 2/3 A • + ■ k + iXk ■ n - iXX 
-s Uk^ - rr?) + y ll^ + zpl + 2y /\ ■ k - 2zp-^ ■ k - ^A^— + i^— j ■ 



(F.l) 



The following change of variable is performed, 



k = k' + zp^-yA + i^ , (F.2) 

k"^ = k'^ + z^ pi + y'^ - 2zPt, ■ k + 2y A ■ k - iA^-^ - 2zyp^ ■ A 

s 

, .. Pn-n .. A-n 

+i\z iXy 



The argument of the exponential (F.l) becomes 



-s (^{k"^ - m^) + y{l - y) + pi z{l - z) + 2zyp^ ■ A 



^ Pt^ ■ n , A ■ n X\ ,^ . 

-iXz— +iXy + iX—] . F.3 

s s s 



The integral over dk is a Gaussian-integral when going to the Euclidean space. Rotating the 
contour of the k^ integration by performing a Wick rotation, i.e. A;^ = ifeg while k = k^;, we are 
left with 

h{X,C,t) = -—J - / dss^^ / dy dz0{l-y-z) 



(27r)4 J 2tt Jq s2 



Jo 



(-m2+A2 y{l-y)+p2 z{l- z)+2zy ■ A-iXz^^ +iXy ^ 



We can now perform the integral over dX. The dependence upon A comes only from the expo- 
nential 

^-iX{X +y A-n- z pTT-n) 
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The expression for the 3-propagator integral hence becomes 

1 



(4^)2 



'1 rl 

ds I dy dzO [1 — y — z) 



10 JO JO 

5iX - zp^-n-yA-n) e-«'»'e-ri^(i-^)^e-^'2y(i-y)«e-22/^sp.-A _ 



In turn, the integral over ds gives 



1 



, dy dze{l-y- z) 

6 {X — zPtj ■ n — yA ■ n) 
m? — z{l — z) — y(l — y) — lyzp-j^ ■ A 



The latter expression can be reduced to a Double Distribution following Eq. (3.25) 

1 



dy dz5{X - zpT,-n-yA-n)l3{z,y,t) ; 



(4vr)2 7o 'Jo 

where the Double Distribution J^^{z,y,t) has be determined to be 

e{l-y-z) 



— p2 2(1 — z) — A'^ y{l — y) — 2yz Pt^ ■ A 



(F.4) 



(F.5) 



(F.6) 



By definition, the DD is not constrained by the kinematics. When linking the 2-dimensional 
Fourier transform that defines the DD to the 1-dimensional Fourier transform defining the skewed 
parton distributions, one has to use the kinematics imposed by the process and to relate the two 



variables p ■ n and A • n. This is done through the J-function in Eq. (F.5). In other words, the 
function I3 depends also on the skewness variable. 

So far, we have hardly used the typical kinematics of the transition and the result is somehow 
general. We can now apply the previous developments to the TDA's kinematics. We will use the 
asymmetric notation introduced by Radyushkin. The kinematics are slightly changed with respect 



to the one described in Section 4.2 In Ref. [H], the following conventions are used 

0, p^-n = l, A-n = 



n 



2 2 



A' 



2pn-A 



(F.7) 



with the asymmetric skewness variable a priori defined between C G [0, 1]. The couple of symmetric 
variables (a;,^) is related the the asymmetric ones {X,C,) by 



X 



1 + e 



c 



(F.e 



Also, in this asymmetric notation, one should change the couple of variables 

(/3,a) (z,(2y + z-l)) , 
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what is equivalent to changing (x,^) for {XX)- 



The expression (F.4) can be written in the form of a reduction formula for the TDA's kinematics 

h{X,C,t) = [\y [\z6{X-z-yOMz,y,t) , (F.9) 

(4vr)^ 7o Jo 

with the double distribution in the same kinematics 

Mz,y,t) = 2 n ^ • (F.IO) 

— z{l — z — y) — t y{l — y — z) 



Performing the z-integration, its value is set at z = X — yQ. In the TDA's kinematics, the value 



of the skewness variable is not constrained to be positive, see (4.9). This leads to two support 
decompositions according to the sign of C,. For positive C, values 

x/c r{i-X)/{i-0\ 



h{x,c,t) = [9{x)e{c-x) +e{i-x)e{x-c) 



(OO) ; 



^2 -ral{X- yC){l - {X - yQ + y) - ty{l - y + {X - yQ) 

(F.ll) 

leading to a support x G [0,C] and x G [Cil]- On the other hand, for negative values of the 
support is decomposed in the following way 



1 / [■{i-X)/{i-0 rii-x)/{i-(r 

/3(x,c,i) = j—^ (eix)e{c-x) / +9ix)9{i-x) / 

(4vr)^ \ J x/c Jo 

dy 



m^-ml{X-yC){l-{X-yC)+y)-ty{l-y + {X-yO) ^ ' 



fF.12) 



which corresponds to a support in x between [(, 0] and [0, 1]. 

The /^''-dependent 3-propagator Integral 

In the case of an integral proportional to the loop variable, i.e. 

d'^k 6{X -k- n) k" 



(2^)4 (fc2 - m2) {{k + A)2 - m2) ((A; - p^)'^ - m?) ' 
the shifting operation changes to y, z dependence of the DD. 

In the particular case of the axial TDA, the Dirac trace taken over the axial matrix element 



brings a factor of —2k^ into the integral. The change of variable (F.2) produces a shift of 

-2F — > -{2k + 2zp^-2yA)'' . 
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Using the fact that the photon polarization impHes e ■ p-f = 0, we obtain e ■ p-,^ = —e ■ A; 



e-h{X,C,t) = i J Jo '^^ Jo 

^_s((fc2-m2)+A2 y{l-y)+pl z{l-z)+2zyp^-ll-i\zl^+i\y^+i\^) 

(F.13) 

The Gaussian integration over dk carries out a ^y^vr/s^ . The term proportional to k cancels out. 
The following steps are as described in above for /3, we obtain 

e-/3(X,C,t) = j\y j\z6{X-z-yO{2y + 2z)D{z,y,t) , (F.14) 



with the double distribution D{z,y,t) defined by Eq. (F.IO) 



Polynomiality Property 



Turning our attention to the polynomiality condition (4.39) resulting from Lorentz invariance, we 
find that it is an inherent property of the Double Distribution parameterization, 



dxx^-^h{x,c,t) 



1 



dy I dze{z + yC) e{l-z- yC) {Cy + zf'^ X^{z, y, t) 



2 7r2 

n-l 



i=0 



The latter expression implies the expression for the general form factors as follows 

Cn,i(t) oc ^ dy dzy' z""^^-' l3{z,y,t) . 



(F.15) 



Since the kinematical range for the skewness variable is — 1 < C ^ Ij the condition on the integra- 
tion limits imposed by the ^-function are verified. 



The expression (F.15 ) corresponds to the Double Distribution form of the General Moments (4.55), 
as expected. 
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